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Abstract. Going beyond the linearized study has been a longstanding problem 
in the theory of Landau damping. In this paper we establish exponential Lan- 
dau damping in analytic regularity. The damping phenomenon is reinterpreted in 
terms of transfer of regularity between kinetic and spatial variables, rather than 
exchanges of energy; phase mixing is the driving mechanism. The analysis involves 
new families of analytic norms, measuring regularity by comparison with solutions 
of the free transport equation; new functional inequalities; a control of nonlinear 
echoes; sharp scattering estimates; and a Newton approximation scheme. Our 
results hold for any potential no more singular than Coulomb or Newton interac- 
tion; the limit cases are included with specific technical effort. As a side result, the 
stability of homogeneous equilibria of the nonlinear Vlasov equation is established 
under sharp assumptions. We point out the strong analogy with the KAM theory, 
and discuss physical implications. 
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Landau damping may be the single most famous mystery of classical plasma 
physics. For the past sixty years it has been treated in the linear setting at various 
degrees of rigor; but its nonlinear version has remained elusive, since the only avail- 
able results [131 ESj prove the existence of some damped solutions, without telling 
anything about their genericity. 

In the present work we close this gap by treating the nonlinear version of Landau 
damping in arbitrarily large times, under assumptions which cover both attractive 
and repulsive interactions, of any regularity down to Coulomb/Newton. 

This will lead us to discover a distinctive mathematical theory of Landau damping, 
complete with its own functional spaces and functional inequalities. Let us make 
it clear that this study is not just for the sake of mathematical rigor: indeed, we 
shall get new insights in the physics of the problem, and identify new mathematical 
phenomena. 

The plan of the paper is as follows. 

In Section [1] we provide an introduction to Landau damping, including historical 
comments and a review of the existing literature. Then in Section [2], we state and 
comment on our main result about "nonlinear Landau damping" (Theorem 12. 6p . 

In Section [3] we provide a rather complete treatment of linear Landau damping, 
slightly improving on the existing results both in generality and simplicity. This 
section can be read independently of the rest. 

In Section |4] we define the spaces of analytic functions which are used in the 
remainder of the paper. The careful choice of norms is one of the keys of our 
analysis; the complexity of the problem will naturally lead us to work with norms 
having up to 5 parameters. As a first application, we shall revisit linear Landau 
damping within this framework. 

In Sections [5] to [7] we establish four types of new estimates (scattering estimates, 
short-term and long-term regularity extortion, echo control); these are the key sec- 
tions containing in particular the physically relevant new material. 

In Section [S] we adapt the Newton algorithm to the setting of the nonlinear Vlasov 
equation. Then in Sections [9] to [11] we establish some iterative estimates along this 
scheme. (Section [11] is devoted specifically to a technical refinement allowing to 
handle Coulomb/Newton interaction.) 

From these estimates our main theorem is easily deduced in Section [12] 

An extension to non-analytic perturbations is presented in Section [13] 
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Some counterexamples and asymptotic expansions are studied in Section [TH 
Final comments about the scope and range of applicability of these results are 
provided in Section [T5l 

Even though it basically proves one main result, this paper is very long. This is 
due partly to the intrinsic complexity and richness of the problem, partly to the need 
to develop an adequate functional theory from scratch, and partly to the inclusion 
of remarks, explanations and comments intended to help the reader understand the 
proof and the scope of the results. The whole process culminates in the extremely 
technical iteration performed in Sections [TD] and [TTJ A short summary of our results 
and methods of proofs can be found in the expository paper [61] . 

This project started from an unlikely conjunction of discussions of the authors 
with various people, most notably Yan Guo, Dong Li, Freddy Bouchet and Etienne 
Ghys. We also got crucial inspiration from the books [H [10] by James Binney and 
Scott Tremaine; and [2] by Serge Alinhac and Patrick Gerard. Warm thanks to 
Julien Barre, Jean Dolbeault, Thierry Gallay, Stephen Gustafson, Gregory Ham- 
mett, Donald Lynden-Bell, Michael Sigal, Eric Sere and especially Michael Kiessling 
for useful exchanges and references; and to Francis Filbet and Irene Gamba for 
providing numerical simulations. We are also grateful to Patrick Bernard, Freddy 
Bouchet, Emanuele Caglioti, Yves Elskens, Yan Guo, Zhiwu Lin, Michael Loss, Peter 
Markowich, Govind Menon, Yann OUivier, Mario Pulvirenti, Jeff Ranch, Igor Rod- 
nianski, Peter Smereka, Yoshio Sone, Tom Spencer, and the team of the Princeton 
Plasma Physics Laboratory for further constructive discussions about our results. 
Finally, we acknowledge the generous hospitality of several institutions: Brown Uni- 
versity, where the first author was introduced to Landau damping by Yan Guo in 
early 2005; the Institute for Advanced Study in Princeton, who offered the second 
author a serene atmosphere of work and concentration during the best part of the 
preparation of this work; Cambridge University, who provided repeated hospitality 
to the first author thanks to the Award No. KUK-Il-007-43, funded by the King 
Abdullah University of Science and Technology (KAUST); and the University of 
Michigan, where conversations with Jeff Ranch and others triggered a significant 
improvement of our results. 



1. Introduction to Landau damping 



1.1. Discovery. Under adequate assumptions (collisionless regime, nonrelativistic 
motion, heavy ions, no magnetic field), a dilute plasma is well described by the 
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nonlinear Vlasov-Poisson equation 




V./ = 0, 



where f = f{t,x,v) > is the density of electrons in phase space {x= position, 
v= velocity), m is the mass of an electron, and F = F[t,x) is the mean-field (self- 
consistent) electrostatic force: 



Here e > is the absolute electron charge, E = E{t, x) is the electric field, and 
p = p{t, x) is the density of charges 



Pi being the density of charges due to ions. This model and its many variants are of 
tantamount importance in plasma physics [H [5l HSl HQ] . 

In contrast to models incorporating collisions [86], the Vlasov-Poisson equation is 
time-reversible. However, in 1946 Landau [H] stunned the physical community by 
predicting an irreversible behavior on the basis of this equation. This "astonishing 
result" (as it was called in [SU]) relied on the solution of the Cauchy problem for 
the linearized Vlasov-Poisson equation around a spatially homogeneous Maxwellian 
(Gaussian) equilibrium. Landau formally solved the equation by means of Fourier 
and Laplace transforms, and after a study of singularities in the complex plane, 
concluded that the electric field decays exponentially fast; he further studied the rate 
of decay as a function of the wave vector k. Landau's computations are reproduced 
in [m Section 34] or [H Section 4.2]. 

An alternative argument appears in ^91 Section 30]: there the thermodynamical 
formalism is used to compute the amount of heat Q which is dissipated when a 
(small) oscillating electric field E(t,x) = E e'^^^'^~''^^'^ {k a wave vector, u; > a 
frequency) is applied to a plasma whose distribution /° is homogeneous in space 
and isotropic in velocity space; the result is 



(1.2) 



F = -cE, 



E = VA^\4np). 



(1.3) 




(1.4) 




where 0(fi) = / P{vi,V2,v^) dv2dvz. In particular, (11. 4p is always positive (see the 
last remark in [IHl Section 30]), which means that the system reacts against the 
perturbation, and thus possesses some "active" stabilization mechanism. 



ON LANDAU DAMPING 



5 



A third argument [IHl Section 32] consists in studying the dispersion relation, 
or equivalently searching for the (generahzed) eigenmodes of the hnearized Vlasov- 
Poisson equation, now with complex frequency uj. After appropriate selection, these 
eigenmodes are all decaying (9a; < 0) as t — *• cxo. This again suggests stability, 
although in a somewhat weaker sense than the computation of heat release. 

The first and third arguments also apply to the gravitational Vlasov-Poisson equa- 
tion, which is the main model for nonrelativistic galactic dynamics. This equation 
is similar to (11. ip . but now m is the mass of a typical star (!), and / is the density 
of stars in phase space; moreover the first equation of (11.21) and the relation (II. 3p 
should be replaced by 



where Q is the gravitational constant, E the gravitational field, and p the density 
of mass. The books by Binney and Tremaine [HI HO] constitute excellent references 
about the use of the Vlasov-Poisson equation in stellar dynamics — where it is 
often called the "coUisionless Boltzmann equation", see footnote on [TUl p. 276]. On 
"intermediate" time scales, the Vlasov-Poisson equation is thought to be an accurate 
description of very large star systems [27], which are now accessible to numerical 
simulations. 

Since the work of Lynden-Bell [53] it has been recognized that Landau damping, 
and wilder coUisionless relaxation processes generically dubbed "violent relaxation" , 
constitute a fundamental stabilizing ingredient of galactic dynamics. Without these 
still poorly understood mechanisms, the surprisingly short time scales for relaxation 
of the galaxies would remain unexplained. 

One main difference between the electrostatic and the gravitational interactions 
is that in the latter case Landau damping should occur only at wavelengths smaller 
than the Jeans length [TOl Section 5.2]; beyond this scale, even for Maxwellian 
velocity profiles, the Jeans instability takes over and governs planet and galaxy 
aggregation^ 

On the contrary, in (classical) plasma physics. Landau damping should hold at 
all scales under suitable assumptions on the velocity profile; and in fact one is in 
general not interested in scales smaller than the Debye length, which is roughly 
defined in the same way as the Jeans length. 

""^or at least would do, if galactic matter was smoothly distributed; in presence of "microscopic" 
heterogeneities, a phase transition for aggregation can occur far below this scale [53]. In the 
language of statistical mechanics, the Jeans length corresponds to a "spinodal point" rather than 
a phase transition [7^ . 
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Nowadays, not only has Landau damping become a cornerstone of plasma physic^, 
but it has also made its way in other areas of physics (astrophysics, but also wind 
waves, fluids, superfiuids,. . . ) and even biophysics. One may consult the concise 
survey papers [71], [751 [HS] for a discussion of its influence and some applications. 

1.2. Interpretation. True to his legend. Landau deduced the damping effect from 
a mathematical-style studjH, without bothering to give a physical explanation of the 
underlying mechanism. His arguments anyway yield exact formulas, which in prin- 
ciple can be checked experimentally, and indeed provide good qualitative agreement 
with observations |54j . 

A first set of problems in the interpretation is related to the arrow of time. In 
the thermodynamic argument, the exterior field is awkwardly imposed from time 
— oo on; moreover, reconciling a positive energy dissipation with the reversibility 
of the equation is not obvious. In the dispersion argument, one has to arbitrarily 
impose the location of the singularities taking into account the arrow of time; at 
mathematical level this is equivalent to a choice of principal value: 

— ^ = P-V- ( - ) + ^7r(5o. 
z — zO \z J 

This is not so serious, but then the spectral study requires some thinking. All in 
all, the most convincing argument remains Landau's original one, since it is based 
only on the study of the Cauchy problem, which makes more physical sense than 
the study of the dispersion relation (see the remark in [9], p. 682]). 

A more fundamental issue resides in the use of analytic function theory, with 
contour integration, singularities and residue computation, which has played a major 
role in the theory of the Vlasov-Poisson equation ever since Landau [l9l Chapter 
32] [TOl Subsection 5.2.4] and helps little, if at all, to understand the underlying 
physical mechanism^ 

The most popular interpretation of Landau damping considers the phenomenon 
from an energetic point of view, as the result of the interaction of a plasma wave 
with particles of nearby velocity [HH p. 18] [IHl p. 412] [H Section 4.2.3] ^ p. 127]. 

^Ryutov |75] estimated in 1998 that "approximately every third paper on plasma physics and 
its applications contains a direct reference to Landau damping" . 

^not completely rigorous from the mathematical point of view, but formally correct, in contrast 
to the previous studies by Landau's fellow physicists — as Landau himself pointed out without 
mercy [15] . 

^Van Kampen [84j summarizes the conceptual problems posed to his contemporaries by Landau's 
treatment, and comments on more or less clumsy attempts to resolve the apparent paradox caused 
by the singularity in the complex plane. 
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In a nutshell, the argument says that dominant exchanges occur with those particles 
which are "trapped" by the wave because their velocity is close to the wave velocity. 
If the distribution function is a decreasing function of \v\, among trapped particles 
more are accelerated than are decelerated, so the wave loses energy to the plasma 
— or the plasma surfs on the wave — and the wave is damped by the interaction. 

Appealing as this image may seem, to a mathematically-oriented mind it will 
probably make little sense at first hearing^ A more down-to-Earth interpretation 
emerged in the fifties from the "wave packet" analysis of Van Kampen [51] and Case 
[H]: Landau damping would result from phase mixing. This phenomenon, well- 
known in galactic dynamics, describes the damping of oscillations occurring when 
a continuum is transported in phase space along an anharmonic Hamiltonian flow 
[ini pp. 379-380]. The mixing results from the simple fact that particles following 
different orbits travel at different ang ulaiS speeds, so perturbations start "spiralling" 
(see Figure 4.27 on [TOl p. 379]) and homogenize by fast spatial oscillation. From 
the mathematical point of view, phase mixing results in weak convergence; from the 
physical point of view, this is just the convergence of observables, deflned as averages 
over the velocity space (this is sometimes called "convergence in the mean"). 

At flrst sight, both points of view seem hardly compatible: Landau's scenario 
suggests a very smooth process, while phase mixing involves tremendous oscillations. 
The coexistence of these two interpretations did generate some speculation on the 
nature of the damping, and on its relation to phase mixing, see e.g. |12] or [lOl 
p. 413]. There is actually no contradiction between the two points of view: many 
physicists have rightly pointed out that that Landau damping should come with 
filamentation and oscillations of the distribution function [Ml p. 962] [IHl p. 141] [H 
Vol. 1, pp. 223-224] [521 pp. 294-295]. Nowadays these oscillations can be visualized 
spectacularly thanks to deterministic numerical schemes, see e.g. [HH] (23 Fig. 3] 
[25] . We reproduce below some examples provided by Filbet. 

In any case, there is still no deflnite interpretation of Landau damping: as noted by 
Ryutov [75l Section 9] , papers devoted to the interpretation and teaching of Landau 
damping were still appearing regularly flfty years after its discovery; to quote just 
a couple of more recent examples let us mention works by Elskens and Escande 
[22| [23l The present paper will also contribute a new point of view. 



Escande Chapter 4, Footnote 6] points out some misconceptions associated with the surfer 
image. 

^"Angular" here refers to action-angle variables, and applies even for straight trajectories in a 
torus. 
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t = 0.16 




-2 2 4 




Figure 1. A slice of the distribution function (relative to a homoge- 
neous equilibrium) for gravitational Landau damping, at two different 
times. 




Figure 2. Time-evolution of the norm of the field, for electrostatic 
(on the left) and gravitational (on the right) interactions. Notice the 
fast Langmuir oscillations in the electrostatic case. 



1.3. Range of validity. The following issues are addressed in the literature [SH 
1121 ESI EH] and slightly controversial: 

• Does Landau damping really hold for gravitational interaction? The case seems 
thinner in this situation than for plasma interaction, all the more that there are 
many instability results in the gravitational context; up to now there has been no 
consensus among mathematical physicists [73] . (Numerical evidence is not conclusive 
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because of the difficulty of accurate simulations in very large time — even in one 
dimension of space.) 

• Does the damping hold for unbounded systems? Counterexamples from [221 ED] 
show that some kind of confinement is necessary, even in the electrostatic case. 
More precisely, Glassey and Schaeffer show that a solution of the linearized Vlasov- 
Poisson equation in the whole space (linearized around a homogeneous equilibrium 
/° of infinite mass) decays at best like 0(t~^), modulo logarithmic corrections, for 
/°(t;) = c/(l + |fp); and like 0((logt)~°) if /° is a Gaussian. In fact. Landau's 
original calculations already indicated that the damping is extremely weak at large 
wavenumbers; see the discussion in 091 Section 32]. Of course, in the gravitational 
case, this is even more dramatic because of the Jeans instability. 

• Does convergence hold in infinite time for the solution of the "full" nonlinear 
equation? This is not clear at all since there is no mechanism that would keep the 
distribution close to the original equilibrium for all times. Some authors do not be- 
lieve that there is convergence as t — > oo; others believe that there is convergence but 
argue that it should be very slow [31], say 0{l/t). In the ffist mathematically rig- 
orous study of the subject. Backus [1] notes that in general the linear and nonlinear 
evolution break apart after some (not very large) time, and questions the validity 
of the linearizationj^ O'Neil [70] argues that relaxation holds in the "quasilinear 
regime" on larger time scales, when the "trapping time" (roughly proportional the 
inverse square root of the size of the perturbation) is much smaller than the damping 
time. Other speculations and arguments related to trapping appear in many sources, 
e.g. [SHI EH]- Kaganovich [H] argues that nonlinear effects may quantitatively affect 
Landau damping related phenomena by several orders of magnitude. 

The so-called "quasilinear relaxation theory" [121 Section 49] [H Section 9.1.2] 
[ISl Chapter 10] uses second-order approximation of the Vlasov equation to predict 
the convergence of the spatial average of the distribution function. The procedure is 
most esoteric, involving averaging over statistical ensembles, and diffusion equations 
with discontinuous coefficients, acting only near the resonance velocity for particle- 
wave exchanges. Because of these discontinuities, the predicted asymptotic state 
is discontinuous, and collisions are invoked to restore smoothness. Linear Fokker- 
Planck equation^ in velocity space have also been used in astrophysics [S31 p. Ill], 

'''From the abstract: "The hnear theory predicts that in stable plasmas the neglected term will 
grow linearly with time at a rate proportional to the initial disturbance amplitude, destroying the 
validity of the linear theory, and vitiating positive conclusions about stability based on it." 

^These equations act on some ensemble average of the distribution; they are different from the 
Vlasov-Landau equation. 



10 



C. MOUHOT AND C. VILLANI 



but only on phenomenological grounds (the ad hoc addition of a friction term leading 
to a Gaussian stationary state); and this procedure has been exported to the study 
of two-dimensional incompressible fluids fi5[ [T6] . 

Even if it were more rigorous, quasilinear theory only aims at second-order cor- 
rections. But the effect of higher order perturbations might be even worse: think of 
something like e~* ^^(£:"^")/^A^T: truncation at any order in e converges exponen- 
tially fast as t ^ oo, but the whole sum diverges to infinity. 

Careful numerical simulation [Sn] seems to show that the solution of the nonlinear 
Vlasov-Poisson equation does converge to a spatially homogeneous distribution, 
but only as long as the size of the perturbation is small enough. We shall call this 
phenomenon nonlinear Landau damping. This terminology summarizes well the 
problem, still it is subject to criticism since (a) Landau himself sticked to the linear 
case and did not discuss the large-time convergence of the data; (b) damping is 
expected to hold when the regime is close to linear, but not necessarily when the 
nonlinear term dominate^; and (c) this expression is also used to designate a related 
but different phenomenon [H Section 10. L3]. 

• Is Landau damping related to the more classical notion of stability in orbital 
sense? Orbital stability means that the system, slightly perturbed at initial time 
from an equilibrium distribution, will always remain close to this equilibrium. Even 
in the favorable electrostatic case, stability is not granted; the most prominent 
phenomenon being the Penrose instability [72] according to which a distribution with 
two deep bumps may be unstable. In the more subtle gravitational case, various 
stability and instability criteria are associated with the names of Chandrasekhar, 
Antonov, Goodman, Doremus, Feix, Baumann, . . . [TOl Section 7.4]. There is a 
widespread agreement (see e.g. the comments in [SH]) that Landau damping and 
stability are related, and that Landau damping cannot be hoped for if there is no 
orbital stability. 

lA. Conceptual problems. Summarizing, we can identify three main conceptual 
obstacles which make Landau damping mysterious, even sixty years after its discov- 
ery: 

(i) The equation is time-reversible, yet we are looking for an irreversible behavior 
as t — >■ -|-oo (or t — ^ — oo). The value of the entropy does not change in time, which 
physically speaking means that there is no loss of information in the distribution 



although phase mixing might stiU play a crucial role in violent relaxation or other unclassified 
nonlinear phenomena. 
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function. The spectacular experiment of the "plasma echo" illustrates this conser- 
vation of microscopic information [HH [SJ]: a plasma which is apparently back to 
equilibrium after an initial disturbance, will react to a second disturbance in a way 
that shows that it has not forgotten the first one0 And at the linear level, if there 
are decaying modes, there also has to be growing modes! 

(ii) When one perturbs an equilibrium, there is no mechanism forcing the system 
to go back to this equilibrium in large time; so there is no justification in the use of 
linearization to predict the large-time behavior. 

(iii) At the technical level. Landau damping (in Landau's own treatment) rests 
on analyticity, and its most attractive interpretation is in terms of phase mixing. 
But both phenomena are incompatible in the large-time limit: phase mixing implies 
an irreversible deterioration of analyticity. For instance, it is easily checked that 
free transport induces an exponential growth of analytic norms as t — oo — except 
if the initial datum is spatially homogeneous. In particular, the Vlasov-Poisson 
equation is unstable (in large time) in any norm incorporating velocity regularity. 
(Space- averaging is one of the ingredients used in the quasilinear theory to formally 
get rid of this instability.) 

How can we respond to these issues? 

One way to solve the first problem (time-reversibility) is to appeal to Van Kam- 
pen modes as in [TOl p. 415]; however these are not so physical, as noticed in [9l 
p. 682]. A simpler conceptual solution is to invoke the notion of weak convergence: 
reversibility manifests itself in the conservation of the information contained in the 
density function; but information may be lost irreversibly in the limit when we con- 
sider weak convergence. Weak convergence only describes the long-time behavior 
of arbitrary observables, each of which does not contain as much information as 
the density function^ As a very simple illustration, consider the time-reversible 
evolution defined by u(t,x) = e^^^Ui^x), and notice that it does converge weakly to 
as t — > ±oo; this convergence is even exponentially fast if the initial datum Ui is 
analytic. (Our example is not chosen at random: although it is extremely simple. 



Interestingly enough, this experiment was suggested as a way to evaluate the strength of 
irreversible phenomena going on inside a plasma, e.g. the collision frequency, by measuring at- 
tenuations with respect to the predicted echo. See [78| for an interesting application and striking 
pictures. 

-"^^In Lynden-Bell's appealing words [551 p. 295], "a system whose density has achieved a steady 
state will have information about its birth still stored in the peculiar velocities of its stars." 
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it may be a good illustration of what happens in phase mixing.) In a way, microso- 
copic reversibility is compatible with macroscopic irreversibility, provided that the 
"microscopic regularity" is destroyed asymptotically. 

Still in respect to this reversibility, it should be noted that the "dual" mechanism 
of radiation, according to which an infinite-dimensional system may lose energy 
towards very large scales, is relatively well understood and recognized as a crucial 
stability mechanism [3], [77] . 

The second problem (lack of justification of the linearization) only indicates that 
there is a wide gap between the understanding of linear Landau damping, and that 
of the nonlinear phenomenon. Even if unbounded corrections appear in the lin- 
earization procedure, the effect of the large terms might be averaged over time or 
other variables. 

The third problem, maybe the most troubling from an analyst's perspective, does 
not dismiss the phase mixing explanation, but suggests that we shall have to keep 
track of the initial time, in the sense that a rigorous proof cannot be based on 
the propagation of some phenomenon. This situation is of course in sharp contrast 
with the study of dissipative systems possessing a Lyapunov functional, as do many 
collisional kinetic equations [5^157] : it will require completely different mathematical 
techniques. 

1.5. Previous mathematical results. At the linear level, the first rigorous treat- 
ments of Landau damping were performed in the sixties; see Saenz [7S] for rather 
complete results and a review of earlier works. The theory was rediscovered and 
renewed at the beginning of the eighties by Degond [20], and Maslov and Fedoryuk 
[58] . In all these works, analytic arguments play a crucial role (for instance for the 
analytic extension of resolvent operators), and asymptotic expansions for the electric 
field associated to the linearized Vlasov-Poisson equation are obtained. 

Also at the linearized level, there are counterexamples by Glassey and Schaeffer 
[29| [30] showing that there is in general no exponential decay for the linearized 
Vlasov-Poisson equation without analyticity, or without confining. 

In a nonlinear setting, the only rigorous treatments so far are those by Caglioti- 
Maffei [13], and later Hwang- Velazquez [3B]. Both sets of authors work in the one- 
dimensional torus and use fixed-point theorems and perturbative arguments to prove 
the existence of a class of analytic solutions behaving, asymptotically as t ^ +oo, 
and in a strong sense, like perturbed solutions of free transport. Since solutions of 
free transport weakly converge to spatially homogeneous distributions, the solutions 
constructed by this "scattering" approach are indeed damped. The weakness of 
these results is that they say nothing about the initial perturbations leading to such 
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solutions, which could be very special. In other words: damped solutions do exist, 
but do we ever reach them? 

Sparse as it may seem, this list is kind of exhaustive. On the other hand, there 
is a rather large mathematical literature on the orbital stability problem, due to 
Guo, Rein, Strauss, Wolansky and Lemou-Mehats-Raphael. In this respect see for 
instance [31] for the plasma case, and [33] for the gravitational case; both sources 
contain many references on the subject. This body of works has confirmed the 
intuition of physicists, although with quite different methods. The gap between 
a formal, linear treatment and a rigorous, nonlinear one is striking: Compare the 
Appendix of [52] to the rest of the paper. In the gravitational case, these works do 
not consider homogeneous equilibria, but only localized solutions. 

Our treatment of Landau damping will be performed from scratch, and will not 
rely on any of these results. 

2. Main result 

2.1. Modelling. We shall work in adimensional units throughout the paper, in d 
dimensions of space and d dimensions of velocity (c/ G N). 

As should be clear from our presentation in Section [T], to observe Landau damping, 
we need to put a restriction on the length scale (anyway plasmas in experiments are 
usually confined). To achieve this we shall take the position space to be the d- 
dimensional torus of sidelength L, namely = / [LZy. This is admittedly a bit 
unrealistic, but it is commonly done in plasma physics (see e.g. [5]). 

In a periodic setting the Poisson equation has to be reinterpreted, since A"~^p is not 
well-defined unless J^a p = 0. The natural solution consists in removing the mean 

value of p, independently of any "neutrality" assumption; in galactic dynamics this 
is known as the Jeans swindle, a trick considered as efficient but logically absurd. 
However, in 2003 Kiessling [H] re-opened the case and acquitted Jeans, on the basis 
that his "swindle" can be justified by a simple limit procedure. In the present case, 
one may adapt Kiessling's argument and approximate the Coulomb potential V 
by some potential exhibiting a "cutoff" at large distances, e.g. of Debye type 
(invoking screening for a plasma, or a cosmological constant for stellar systems; 
anyway the particular choice of approximation has no influence on the result). If 
Wk £ L^(M'^), then VV^ * p makes sense for a periodic p, and moreover 



(VK*p)(x)=/ VV,{x-y)p{y)dy= [ VV^'^^x - y) p{y) dy, 
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where Vk^\z) = ^^^^d Vk{z + £L). Passing to the hmit as k — > yields 

f VV^^)(x-y)p{y)dy= [ VV^^){x-y) {p- {p)){y) dy = -VAl'{p- (p)), 
J [0,1]'' AoM"^ 

where A^^ is the inverse Laplace operator on T^. We refer to [H] for a discussion 
of the physics underlying this limit k —>■ 0. 

More generally, we may consider any interaction potential W on T^, satisfying 
certain regularity assumptions. Then the self-consistent field will be given by 



-VW * p, p{x) = j f{x,v)dv, 



where now * denotes the convolution on T^. 

In accordance with our conventions from Appendix IA.3t we shall write W^^\k) = 
Jjd e~^^'"^'L W{x) dx. In particular, if W is the periodization of a potential M'^ ^ R 

(still denoted W by abuse of notation), i.e., 

W{x) = W^^\x) = J2 + ^L), 

then 

(2.1) w^i^){k) = w(^, 

where W{^^ = J^^ e~'^^'^^''^ W (x) dx is the original Fourier transform in the whole 
space. 

2.2. Linear damping. It is well-known that Landau damping requires some stabil- 
ity assumptions on the unperturbed homogeneous distribution function, say /°(f). 
In this paper we shall use a very general assumption, expressed in terms of the 
Fourier transform 

(2.2) f{r^)= I e-'-^-^-^ f{v)dv, 



the length L, and the interaction potential W . To state it, we define, for t > and 
k e 

(2.3) K\t,k) = -4n'W^^\k)f ^t; 
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and, for any ^ E C, we define a function C via the following Fourier-Laplace trans- 
form of K'^ in the time variable: 

(2.4) C{i,k)= e^^^-t' K'^{t,k)dt, 

Jo 

where ^* is the complex conjugate to ^. Our linear damping condition is expressed 
as follows: 



(L) There are constants Co, A, k > such that for any rj G M"^, |/°('7)| < 
Co e-2'^^l''l; and for any ^ e C with 0<^^<X, 

inf \C{^,k) -1\>K. 

We shall prove in Section [3] that (L) implies Landau damping. For the moment, 
let us give a few sufficient conditions for (L) to be satisfied. The first one can 
be thought of as a smallness assumption on either the length, or the potential, or 
the velocity distribution. The other conditions involve the marginals of f^ along 
arbitrary wave vectors k: 

(2.5) ^k{v)= [ f{w)dw, veR. 

All studies known to us are based on one of these assumptions, so (L) appears as a 
unifying condition for linear Landau damping around a homogeneous equilibrium. 

Proposition 2.1. Let /° = f^{v) be a velocity distribution such that f^ decays 
exponentially fast at infinity, let L > and let W be an interaction potential on T^, 
W E {T'^) . If any one of the following conditions is satisfied: 

(a) smallness: 

(2.6) 47rM max|iy(-^)(A;)| sup / |/°( ra) \ rdr] < 1; 

(b) repulsive interaction and decreasing marginals: for all k eU^ and v G M, 



(2.7) W^^\k)>0; 




V'kiv) > 
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(c) generalized Penrose condition on marginals: for all k G Z*^, 

(2.8) VwgM, v[.iw) = =^ W^^\k) (p.v. f ^^dv] < 1; 

V JrV-w 

then (L) holds true for some Cq, A, k > 0. 

Remark 2.2. [l9l Problem, Section 30] If f^ is radially symmetric and positive, 
and d > 3, then all marginals of f^ are decreasing functions of \v\. Indeed, if 
(p{v) = /jgd-i /(a/w^ + |wP) dw, then after differentiation and integration by parts 
we find 

'if'iv) = -{d-3)v f f{Vv' + \w\^) ^ {d>^) 
y{v) = -27rvf{\v\) {d = 3). 
Example 2.3. Take a gravitational interaction and Mawellian background: 



2T 



Recalling (12. ip . we see that (12. 6p becomes 



(2.9) ^<y0-^^(^'A 

The length Lj is the celebrated Jeans length [IDlIll], so criterion (a) can be applied, 
all the way up to the onset of the Jeans instability. 

Example 2.4. If we replace the gravitational interaction by the electrostatic interac- 
tion, the same computation yields 



(2.10) L<,/^=:L^(T,p°), 

and now is essentially the Debye length. Then criterion (a) becomes quite 
restrictive, but because the interaction is repulsive we can use criterion (b) as soon 
as /° is a strictly monotone function of \v\\ this covers in particular Maxwellian 
distributions, independently of the size of the box. Criterion (b) also applies if d > 3 
and has radial symmetry. For given L > 0, the condition (L) being open, it will 
also be satisfied if f^ is a small (analytic) perturbation of a profile satisfying (b); 
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this includes the so-called "small bump on tail" stability. Then if the distribution 
presents two large bumps, the Penrose instability will take over. 

Example 2.5. For the electrostatic interaction in dimension 1, (12.81) becomes 

(2.11) {fy(w) = o^ fil!!)Mdv< 



This is a variant of the Penrose stability condition [72]. This criterion is in general 
sharp for linear stability (up to the replacement of the strict inequality by the non- 
strict one, and assuming that the critical points of are nondegenerate) ; see [50, 
Appendix] for precise statements. 

We shall show in Section [3] that (L) implies linear Landau damping (Theorem 
13. ip : then we shall prove Proposition 12.11 at the end of that section. The general 
ideas are close to those appearing in previous works, including Landau himself; the 
only novelties lie in the more general assumptions, the elementary nature of the 
arguments, and the slightly more precise quantitative results. 

2.3. Nonlinear damping. As others have done before in the study of Vlasov- 
Poisson [13], we shall quantify the analyticity by means of natural norms involving 
Fourier transform in both variables (also denoted with a tilde in the sequel). So we 
define 

(2.12) ll/IU,, = sup (e^-^l"! e2-^^|7(^)(A;,r/) 

where k varies in /i are positive parameters, and we recall the depen- 

dence of the Fourier transform on L (see Appendix I A. 31 for conventions). Now we 
can state our main result as follows: 

Theorem 2.6 (Nonlinear Landau damping). Let /° : M*^ ^ R+ be an analytic 
velocity profile. Let L > and W : T'l ^ W be an interaction potential satisfying 

(2.13) ykeZ', \W^''\k)\<j^ 

for some constants Cw > 0, 7 > 1. Assume that /° and W satisfy the stability 
condition (L) from Subsection \2.^ with some constants A, k > 0; further assume 
that, for the same parameter \, 

/ ~ \ A" 

(2.14) sup(|r(r/)|e2'^^l''l) <Co, V - ||V:;r|Ui(M.) < Co < +00. 
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Then for any < A' < A, /3 > 0, < /i' < fi, there is e = e{d, L, Cw, Co, k, A, A', /i, fi', 
with the following property: if fi = fi{x,v) is an initial datum satisfying 

(2.15) 5:= ||/^_/0||,^+ ff \f,-f'\e^\''\dvdx<e, 

then 



the unique classical solution f to the nonlinear Vlasov equation 
df 



(2.16) 



dt 



+ V ■ VJ - {VW * p) ■ VJ = 0, 



P 



fdv, 



with initial datum /(O, ■) = fi, converges in the weak topology as t ±oo, with 
rate 0(e~^'^'^''*'), to a spatially homogeneous equilibrium /±oo/ 

• the density p{t,x) = J f{t,x,v)dv converges in the strong topology as t ^ ±oo, 
with rate 0(e~^'''^'l*l), to the constant density 

1 



fi{x,v)dx dv; 



in particular the force F = —VW * p converges exponentially fast to 0. 

• the space average {f)(t,v) = J f{t,x,v)dx converges in the strong topology as 

t ±oo, with rate ©(e^^'^^'l*!), to f±^. 

More precisely, there are C > 0, and spatially homogeneous distributions /+oo(f) 
and f-oo{v), depending continuously on fi and W , such that 



(2.17) 

and 



sup 



f{t,x + vt,v)-f\v) <C6- 

X',fj,' 



V(fc,r/) e X R'^, 

(2.18) Vr G N, 

(2.19) Vr e N, 



L-'^/(^)(t,A;,r/)-/+oo(r/)U=o| =0(e-2-^*) as t ^ +00; 
L-'^7(^)(t,fc,r/)-7_oo(r/)l,=o| =0(e-2-^l*l) as t ^ -00; 

\\p{t, ■) - PooWcr^j,) = 0(e-2"Tl*l) as \t\ 00; 

\\Fit, ■ )||c^(T^) = 0(e~2"i^l*l) as \t\ ^ 00; 
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(2.20) Vr e N, Va > 0, (/(t, -, t;))-/±oo ^^^^^^ = 0(e-2"^l*l) as t ±oo 



In this statement C" stands for the usual norm on r times continuously dif- 
ferentiable functions, and involves in addition moments of order a, namely 
ll/llcj = sup,^/<^,„gigd \f^^'\v) (1 + |f l"^)!- These results could be reformulated in 
a number of alternative norms, both for the strong and for the weak topology. 

2.4. Comments. Let us start with a list of remarks about Theorem 12.61 

• The decay of the force field, statement fl2.19p . is the experimentally measurable 
phenomenon which may be called Landau damping. 

• Since the energy 

E = ^JJ p{x) p{y)W{x-y)dxdy + J f{x,v)^-^dvdx 

(= potential + kinetic energy) is conserved by the nonlinear Vlasov evolution, there 
is a conversion of potential energy into kinetic energy as t ^ oo (kinetic energy 
goes up for Coulomb interaction, goes down for Newton interaction). Similarly, the 
entropy 

S = - [f f\ogf = -([ p\ogp+ f /log^ 



P, 

(= spatial + kinetic entropy) is preserved, and there is a transfer of information 
from spatial to kinetic variables in large time. 

• Our result covers both attractive and repulsive interactions, as long as the linear 
damping condition is satisfied; it covers Newton/Coulomb potential as a limit case 
(7 = 1 in fl2.13p ). The proof breaks down for 7 < 1; this is a nonlinear effect, as any 
7 > would work for the linearized equation. The singularity of the interaction at 
short scales will be the source of important technical problemsEl 

• Condition fl2.14p could be replaced by 

(2.21) \f%r))\ < [ f{v)e^\''\dv < Co- 



But condition fl2.14p is more general, in view of Theorem 14.201 below. For instance, 
f^{v) = 1/(1 + ^;^) in dimension d = 1 satisfies fl2.14p but not fl2.2ip : this distribution 
is commonly used in theoretical and numerical studies, see e.g. [37j. We shall also 



-'^^In a related subject, this singularity is also the reason why the Vlasov-Poisson equation is still 
far from being established as a mean-field limit of particle dynamics (see [35] for partial results 
covering much less singular interactions). 
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establish slightly more precise estimates under slightly more stringent conditions on 
f, see dmi). 

• Our conditions are expressed in terms of the initial datum, which is a consid- 
erable improvement over [T3l [38] . Still it is of interest to pursue the "scattering" 
program started in [13], e.g. in a hope of better understanding of the nonperturba- 
tive regime. 

• The smallness assumption on /j — is expected, for instance in view of the 
work of O'Neil [7D], or the numerical results of [53] • We also make the standard 
assumption that fi — is well localized. 

• No convergence can be hoped for if the initial datum is only close to /° in 
the weak topology: indeed there is instability in the weak topology, even around a 
Maxwellian [T3] . 

• Strictly speaking, known existence and uniqueness results for solutions of the 
nonlinear Vlasov-Poisson equation P, [5l] do not apply to the present setting of 
close-to-homogeneous analytic solutions. (The problem with [B] is that velocities are 
assumed to be uniformly bounded, and the problem with [SI] is that the position 
space is the whole of Mf^; in both papers these assumptions are not superficial.) 
However, this really is not a big deal: our proof will provide an existence theorem, 
together with regularity estimates which are considerably stronger than what is 
needed to prove the uniqueness. We shall not come back to these issues which 
are rather irrelevant for our study: uniqueness only needs local in time regularity 
estimates, while all the difficulty in the study of Landau damping consists in handling 
(very) large time. 

• f{t, ■) is not close to /° in analytic norm as t oo, and does not converge to 
anything in the strong topology, so the conclusion cannot be improved much. Still 
we shall establish more precise quantitative results, and the limit profiles /±oo are 
obtained by a constructive argument. 

• Estimate fl2.17p expresses the orbital "travelling stability" around /°; it is much 
stronger than the usual orbital stability in Lebesgue norms [33] • An equivalent 
formulation is that if (Tt)t^^ stands for the nonlinear Vlasov evolution operator, and 
{T^)t(z^ for the free transport operator, then in a neighborhood of a homogeneous 
equilibrium satisfying the stability criterion (L), oTj remains uniformly close to 
Id for all t. Note the important difference: unlike in the usual orbital stability theory, 
our conclusions are expressed in functional spaces involving smoothness, which are 
not invariant under the free transport semigroup. This a source of difficulty (our 
functional spaces are sensitive to the filamentation phenomenon), but it is also the 
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reason for which this "analytic" orbital stability contains much more information, 
and in particular the damping of the density. 

• Compared with known nonlinear stability results, and even forgetting about 
the smoothness, estimate fl2.17p is new in several respects. In the context of plasma 
physics, it is the first one to prove stability for a distribution which is not necessarily a 
decreasing function of ( "small bump on tail" ); while in the context of astrophysics, 
it is the first one to establish stability of a homogeneous equilibrium against periodic 
perturbations with wavelength smaller than the Jeans length. 

• While analyticity is the usual setting for Landau damping, both in mathematical 
and physical studies, it is natural to ask whether this restriction can be dispended 
with. (This can be done only at the price of losing the exponential decay.) In the 
linear case, this is easy, as we shall recall later in Remark 13. 5t but in the nonlinear 
setting, leaving the analytic world is much more tricky. In Section [131 we shall 
present the first results in this direction. 

With respect to the questions raised above, our analysis brings the following 
answers: 

(a) Convergence of the distribution / does hold for t +oo; it is indeed based on 
phase mixing, and therefore involves very fast oscillations. In this sense it is right to 
consider Landau damping as a "wild" process. But on the other hand, the spatial 
density (and therefore the force field) converges strongly and smoothly. 

(b) The space average (/) does converge in large time. However the conclusions are 
quite different from those of quasilinear relaxation theory, since there is no need for 
extra randomness, and the limiting distribution is smooth, even without collisions. 

(c) Landau damping is a linear phenomenon, which survives nonlinear pertur- 
bation thanks to the structure of the Vlasov-Poisson equation. The nonlinearity 
manifests itself by the presence of echoes. Echoes were well-known to specialists of 
plasma physics [HI Section 35] [1], Section 12.7], but were not identified as a possible 
source of unstability. Controlling the echoes will be a main technical difficulty; but 
the fact that the response appears in this form, with an associated time-delay and 
localized in time, will in the end explain the stability of Landau damping. These 
features can be expected in other equations exhibiting oscillatory behavior. 

(d) The large-time limit is in general different from the limit predicted by the 
linearized equation, and depends on the interaction and initial datum (more precise 
statements will be given in Section [T^ : still the linearized equation, or higher-order 
expansions, do provide a good approximation. We shall also set up a systematic 
recipe for approximating the large-time limit with arbitrarily high precision as the 
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strength of the perturbation becomes smalL This justifies a posteriori many known 
computations. 

(e) From the point of view of dynamical systems, the nonhnear Vlasov equation 
exhibits a truly remarkable behavior. It is not uncommon for a Hamiltonian system 
to have many, or even countably many heteroclinic orbits (there are various theories 
for this, a popular one being the Melnikov method); but in the present case we see 
that heteroclinic/homoclinic orbit^ are so numerous as to fill up a whole neigh- 
borhood of the equilibrium. This is possible only because of the infinite-dimensional 
nature of the system, and the possibility to work with nonequivalent norms; such 
a behavior has already been reported for other systems ^61 W7\ . in relation with 
infinite-dimensional KAM theory. 

(f) As a matter of fact, the nonlinear Landau damping has strong similarities 
with the KAM theory. It has been known since the early days of the theory that 
the linearized Vlasov equation can be reduced to an infinite system of uncoupled 
Volterra equations, which makes this equation completely integrable in some sense. 
(Morrison [HI] gave a more precise meaning to this property.) To see a parallel 
with classical KAM, one step of our result is to prove the preservation of the phase- 
mixing property under nonlinear perturbation of the interaction. (Although there 
is no ergodicity in phase space, the mixing will imply an ergodic behavior for the 
spatial density.) The analogy goes further since the proof of Theorem 12.61 shares 
many features with the proof of the KAM theorem (closest to Kolmogorov's original 
version, see [18] for a complete exposition). 

Thus we see that three of the most famous paradoxical phenomena from twen- 
tieth century classical physics: Landau damping, echoes, and KAM theorem, are 
intimately related (only in the nonlinear variant of Landau's linear argument!). This 
relation, which we did not expect, is one of the main discoveries of the present paper. 

2.5. Interpretation. A successful point of view adopted in this paper is that Lan- 
dau damping is a relaxation by smoothness and by mixing. In a way, phase 
mixing converts the smoothness into decay. Thus Landau damping emerges as a rare 
example of a physical phenomenon in which regularity is not only crucial from the 
mathematical point of view, but also can be "measured" by a physical experiment. 

2.6. Main ingredients. Some of our ingredients are similar to those in [13]: in par- 
ticular, the use of Fourier transform to quantify analytic regularity and to implement 
phase mixing. New ingredients used in our work include 



Here we use these words just to designate solutions connecting two distinct/equal equilibria, 
without any mention of stable or unstable manifolds. 
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• the introduction of a time-shift parameter to keep memory of the initial time 
(Sections m and Ej), thus getting uniform estimates in spite of the loss of regularity 
in large time. We call this the gliding regularity: it shifts in phase space from low 
to high modes. Gliding regularity automatically comes with an improvement of the 
regularity in x, and a deterioration of the regularity in v, as time passes by. 

• the use of carefully designed flexible analytic norms behaving well with respect 
to composition (Section H]). This requires care, because analytic norms are very 
sensitive to composition, contrary to, say, Sobolev norms. 

• "finite-time scattering" at the level of trajectories to reduce the problem to 
homogenization of free flow (Section [5]) via composition. The physical meaning is 
the following: when a background with gliding regularity acts on (say) a plasma, 
the trajectories of plasma particles are asymptotic to free transport trajectories. 

• new functional inequalities of bilinear type, involving analytic functional spaces, 
integration in time and velocity variables, and evolution by free transport (Section 
E]). These inequalities morally mean the following: when a plasma acts (by forcing) 
on a smooth background of particles, the background reacts by lending a bit of its 
(gliding) regularity to the plasma, uniformly in time. Eventually the plasma will 
exhaust itself (the force will decay). This most subtle effect, which is at the heart 
of Landau's damping, will be mathematically expressed in the formalism of analytic 
norms with gliding regularity. 

• a new analysis of the time response associated to the Vlasov-Poisson equation 
(Section [7]), aimed ultimately at controlling the self-induced echoes of the plasma. 
For any interaction less singular than Coulomb/Newton, this will be done by an- 
alyzing time-integral equations involving a norm of the spatial density. To treat 
Coulomb/Newton potential we shall refine the analysis, considering individual modes 
of the spatial density. 

• a Newton iteration scheme, solving the nonlinear evolution problem as a succes- 
sion of linear ones (Section [TOl) . Picard iteration schemes still play a role, since they 
are run at each step of the iteration process, to estimate the scattering operators. 

It is only in the linear study of Section [3] that the length scale L will play a crucial 
role, via the stability condition (L). In all the rest of the paper we shall normalize 
L to 1 for simplicity. 

2.7. About phase mixing. A physical mechanism transferring energy from large 
scales to very fine scales, asymptotically in time, is sometimes called weak turbu- 
lence. Phase mixing provides such a mechanism, and in a way our study shows that 
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the Vlasov-Poisson equation is subject to weak turbulence. But the phase mixing 
interpretation provides a more precise picture. While one often sees weak turbulence 
as a "cascade" from low to high Fourier modes, the relevant picture would rather 
be a two-dimensional figure with an interplay between spatial Fourier modes and 
velocity Fourier modes. More precisely, phase mixing transfers the energy from each 
nonzero spatial frequency k, to large velocity frequences 77, and this transfer occurs 
at a speed proportional to k. This picture is clear from the solution of free transport 
in Fourier space, and is illustrated in Fig. [31 (Note the resemblance with a shear 
flow.) So there is transfer of energy from one variable (here x) to another (here v); 
homogenization in the first variable going together with filamentation in the second 
one. The same mechanism may also underlie other cases of weak turbulence. 

k 

(spatial ihodes) t = h t = h t = h 



»-€^S>' 



(kinetic modes) 

initial configuration 
(t = 0) 

Figure 3. Schematic picture of the evolution of energy by free trans- 
port, or perturbation thereof; marks indicate localization of energy in 
phase space. 

Whether ultimately the high modes are damped by some "random" microscopic 
process (collisions, diffusion, . . . ) not described by the Vlasov-Poisson equation 
is certainly undisputed in plasma physics ^91 Section 410, and is the object of 



""^^See [ini Problem 41]: thanks to Landau damping, collisions are expected to smooth the 
distribution quite efficiently; this is a hypoelliptic problematic. 
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Figure 4. The distribution function in phase space (position, veloc- 
ity) at a given time; notice how the fast oscillations in v contrast with 
the slower variations in x. 

debate in galactic dynamics; anyway this is a different story. Some mathematical 
statistical theories of Euler and Vlasov-Poisson equations do postulate the existence 
of some small-scale coarse graining mechanism, but resulting in mixing rather than 
dissipation [TH [83] . 

3. Linear damping 

In this section we establish Landau damping for the linearized Vlasov equation. 
Beforehand, let us recall that the free transport equation 



has a strong mixing property: any solution of (13. ip converges weakly in large time 
to a spatially homogeneous distribution equal to the space-averaging of the initial 
datum. Let us sketch the proof. 

If / solves (El]) in T'^ x M'^, with initial datum /, = /(O, then f{t,x,v) = 
fi{x — vt, v), so the space- velocity Fourier transform of / is given by the formula 




^ + w ■ V./ = 




f{t,k,r]) = fi{k,r] + kt). 



On the other hand, if foo is defined by 
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then foo{k,ri) = fi{0,ri) lk=o- So, by the Riemann-Lebesgue lemma, for any fixed 
{k, rf) we have 



which shows that / converges weakly to foo- The convergence holds as soon as fi is 
merely integrable; and by (13.21) . the rate of convergence is determined by the decay 
of fiik^T]) as |?7| — > oo, or equivalently the smoothness in the velocity variable. In 
particular, the convergence is exponentially fast if (and only if) fi{x,v) is analytic 
in V. 

This argument obviously works independently of the size of the box. But when 
we turn to the Vlasov equation, length scales will matter, so we shall introduce a 
length L > 0, and work in = W^/{Ll/). Then the length scale will appear in 
the Fourier transform: see Appendix IA.3I (This is the only section in this paper 
where the scale will play a nontrivial role, so in all the rest of the paper we shall 
take L = 1.) 

Any velocity distribution /° = f^{v) defines a stationary state for the nonlin- 
ear Vlasov equation with interaction potential W. Then the linearization of that 
equation around /° yields 



Note that there is no force term in (13.31) . due to the fact that does not depend 
on X. This equation describes what happens to a plasma density / which tries to 
force a stationary homogeneous background equivalently, it describes the reaction 
exerted by the background which is acted upon. 

Theorem 3.1 (Linear Landau damping). Let /° = f^{v), L > 0, W : T'l ^ M. such 
that ||VW^||li < Cw < +oo, and fi{x,v) such that 

(i) Condition (L) from Subsection \2.2\ holds for some constants X, k, > 0; 

(ii) \/7] e R"^, \f{v)\ < Co e-^'^^l^'l for some constant Cq > 0; 

(ill) VA; G Z"', Vr/ G M"', \f}^\k,r])\ < e-^™!*)! for some constants a > 0, 



f{t,k,v) - foo{k,r]) 



\t\~too 



-0, 



(3.3) 




a>o. 
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Then as t ^ +00 the solution f{t, ■) to the linearized Vlasov equation (13.31) with 
initial datum fi converges weakly to foo = {fi) defined by 

/00(f) = Td I fi{x,v) dx; 
and p{x) = J f{x,v)dv converges strongly to the constant 

P°° = Td 11 fii^^v) dxdv. 



More precisely, for any X' < min{A ; a}, 



Vr e N, 



''00 I I Qr 



2tv\' I 

e L I 



') 



f^'^\t,k,r^)-f^\k,r^) =0{e 



Iktl 



Remark 3.2. Even if the initial datum is more regular than analytic, the conver- 
gence will in general not be better than exponential (except in some exceptional 
cases [36]). See fiUi pp. 414-416] for an illustration. Conversely, if the analyticity 
width a for the initial datum is smaller than the "Landau rate" A, then the rate of 
decay will not be better than 0(e~"*). See [TJ [19] for a discussion of this fact, often 
overlooked in the physical literature. 

Remark 3.3. The fact that the convergence is to the average of the initial da- 
tum will not survive nonlinear perturbation, as shown by the counterexamples in 
Subsection [HI 

Remark 3.4. Dimension does not play any role in the linear analysis. This can be 
attributed to the fact that only longitudinal waves occur, so everything happens "in 
the direction of the wave vector" . Transversal waves arise in plasma physics only 
when magnetic effects are taken into account [H Chapter 5]. 

Remark 3.5. The proof can be adapted to the case when and fi are only C°°; 
then the convergence is not exponential, but still 0{t~°°). The regularity can also be 
further decreased, down to 1^'^, at least for any s > 2; more precisely, if /° G W^°'^ 
and fi G W^"^ there will be damping with a rate 0{t~'^) for any k, < max{so — 2 ; Sj}. 
(Compare with [H Vol. 1, p. 189].) This is independent of the regularity of the 
interaction. 



28 



C. MOUHOT AND C. VILLANI 



The proof of Theorem 13. II relies on the following elementary estimate for Volterra 
equations. We use the notation of Subsection 12. 2[ 

Lemma 3.6. Assume that (L) holds true for some constants Cq, k, X > 0; let Cw = 
\\W\\ii(jd^ and let he defined by (12.31) . Then any solution (p{t,k) of 

(3.4) ^{t,k) =a{t,k)+ [ K^{t - r.k) ^{r.k) dr 

Jo 

satisfies, for any k E'L'^ and any \' < \, 
sup(|^(t,A;)|e2-^'^*) < 

Here C(A, A', k) = C {1 + hi~^{l + (A — A')"-*^/^)) for some universal constant C. 

Remark 3.7. It is standard to solve these Volterra equations by Laplace trans- 
form; but, with a view to the nonlinear setting, we shall prefer a more flexible and 
quantitative approach. 

Proof of Lemma \3. 6[ li k = this is obvious since K^{t, 0) = 0; so we assume k ^ 0. 
Consider A' < A, multiply (13.41) by e^'^^ "t*, and write 

$(t, k) = ip{t, k) e2-^'^*, A{t, k) = a{t, k) e^-^'^*; 
then (13. 4p becomes 

(3.5) ^t,k) = A{t,k)+ K\t-T,k)e^^^'^-^^'"^H{T,k)dT. 

Jo 

A particular case: The proof is extremely simple if we make the stronger assump- 
tion 

/ \K%T,k)\e^''^ -^^dT<l-K, «G(0,1). 
Jo 

Then from fl331) . 

sup |$(t. A;) I < sup \A{t, k)\ 

0<t<T 0<t<T 

+ sup ( f \K\t - r, /t)| e^-^'^^*-") dA sup |$(r, k)\, 

0<t<T \Jo J 0<T<T 



1 + CoCwC{X,X',k) 



sup 



a(t, k)\ 
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whence 



sup |y4(r, A;)| sup |74(r, /c)| 

sup |$(r,A;)| < -—^ ■ < - - 



f+OO 



and therefore 



sup fe2-^'^*|<^(t,A;)|) < (-\ supf |a(t, e^-^'^*). 



The general case: To treat the general case we take the Fourier transform in the 
time variable, after extending A and $ by at negative times. (This presentation 
was suggested to us by Sigal, and appears to be technically simpler than the use of 
the Laplace transform.) Denoting the Fourier transform with a hat and recalling 
<^^, we have, for i = \' + iujL/\k\, 

$(cu, k) = A{uj, k) + C{i, k) $(cj, k). 
By assumption C{^, k) ^ 1, so 

A{uj, k) 



^{uj,k) 



c{^,ky 



From there, it is traditional to apply the Fourier (or Laplace) inversion transform. 
Instead, we apply Plancherel's identity to find (for each k) 

ll^ll ^ \\A\\LHdt) _ 

and then we plug this in the equation (13.51) to get 

(3.6) IML^idt) < \\A\\L^idt) + i|A^°e2'^^'^*|U2(rfi) IMi^dt) 
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It remains to bound the second term. On the one hand, 



(3.7) 



\A\ 



|a(t, k)\' e^""^ dt 



1/2 



< 



-47r(A-A')^t 



L 



.1/2 




I supl a(t, 


/ t>0 ^ 





On the other hand, 



(3. 



„2TTX'^t\ 



\L'2{dt) 



47r|A:|(A-A') 



sup( \a{t, k) \ e 



\kl 

L2 



g4^A'Mt 



1/2 



AttX'u I 



\f{au)\'u'du 



1/2 



where a = k/\k\. The estimate follows immediately. (Note that the factor |A;| in 
fl371) cancels with \k\^/' in ([33]).) 

It seems that we only used properties of the function £ in a strip 3?^ ~ A; but this 
is an illusion. Indeed, we have taken the Fourier transform of $ without checking 
that it belongs to (L^ + L'){dt), so what we have established is only an a priori 
estimate. To convert it into a rigorous result, one can use a continuity argument after 
replacing A' by a parameter a which varies from — e to A'. (By the integrability of 
and Gronwall's lemma, (p is obviously bounded as a function of t; so (p{k,t) e"'^''^'*/^ 
is integrable for any e > 0, and continuous as e — > 0.) Then assumption (L) 
guarantees that our bounds are uniform in the strip < 3?.^ < A', and the proof goes 
through. □ 

Proof of Theorem \3.1\ Without loss of generality we assume t > 0. Considering 
(13. 3p as a perturbation of free transport, we apply Duhamel's formula to get 



(3.9) f{t,x,v) = f\{x~vt,v) + 
Integration in v yields 



/ [iyW * p) ■ V./°] (r, X - v(t -t),v) dr. 
Jo 



(3.10) pit,x) 



f,{x-vt,v)dv+ [ [ [{VW*p)-V^f°]{T,x-v{t-T),v) dvdr. 

Jo JR'i 



Of course, / p(t,x) dx = JJ fi{x,v) dx dv. 
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For k ^ 0, taking the Fourier transform of (13.101) . we obtain 

p(L)(t,k)= [ I fi{x-vt,v)e-^'''T^-^dvdx 

Jo Jt^ Jr'^ 
= [ [ fi{x,v)e-^'''i-''e-^'''i-''Uvdx 

+ /* / / [i^W*p)-Vj'']{T,x,v)e-^'^T-e-^'^i-''^'-^UvdxdT 
Jo Jjf Jr'' 

= (k, I) + I\VW « pp)(r, k) ■ Vjo (^^) dr 

In conclusion, we have estabhshed the closed equation on p^^^: 
(3.11) p<-)(U) = F)(fc,^) 

- 4vr^ W^(^) (A:) p^^) (r, fc) f (^^) ^ " ^) 
Recalhng (12.31) . this is the same as 

f^'\t, k) = f^'^ (^k, + 1^ K%t - r, k) p<^)(r, k) dr. 
Without loss of generality, \ < a. By Assumption (L) and Lemma [3. 6 [ 

In particular, for k ^ we have 

Vt>l, |p^^)(t,fc)|=0(e-^*e-^^l^l); 

2-!t\" 

so any Sobolev norm of p — poo converges to zero like 0{e ^*), where A" is arbi- 
trarily close to A' and therefore also to A. By Sobolev embedding, the same is true 
for any C" norm. 
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Next, we go back to fl3.9p and take the Fourier transform in both variables x and 
f , to find 




fi{x - vt, v) e-^'^T^-'' e-^'^"^-" dxdv 




m{VW * p) (r, x-v{t- r)) • Vj\v) e-^*"*'" e-^^""''' dx dv dr 
-d 

v) e-^'^i-"" e-^'^i-""' e-'^^'^''' dx dv 

miyW * p){t, x) ■ VJ%v) g-^-r-C*-) g-^^-'?-' dxdv dr 

jd 



kt 



So 

(3.12) 



kt 



fiL)[t,k,rj-^)=fi''\k,rj)+ I VW''''{k)p('^\T,k)-VJ^iv-^] dr. 



In particular, for any rj G Mf^, 

(3.13) p\t,o,v) = F\o,v); 

in other words, (/) = J f dx remains equal to (/j) for all times. 
On the other hand, if /c 7^ 0, 

(3.14) 

{t,k,v-j)\<\n'\k,v)\ 



kr 



L 



+ I \VW^'^\k)\ \p^''\T,k) 
Jo 

< C i e~2™l''l - 



kr 



dr 



kr 

T 



-^'^^l""^! dr 



27rA'Mr -2tt 



(a'+A) I _feT 
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where we have used A' < (A' + A)/2 < A, and C only depends on Cw, Ci, A, A', k. 
In the end, 







< 



7r(A - A') 

Plugging this back in flXTIj) . we obtain, with A" = A' - (A - A')/2, 



< Ce 



~2ttX" 



(3.15) 

In particular, for any fixed r] and k ^ 0, 

\P\t,k,v)\ < C e~''''"\'^+'^\ = 0{e-'<\'\). 

We conclude that /^^^ converges pointwise, exponentially fast, to the Fourier trans- 
form of (fi). Since A' and then A" can be taken as close to A as wanted, this ends 
the proof. □ 

We close this section by proving Proposition 12.11 

Proof of Proposition \2.1\ First assume (a). Since /° decreases exponentially fast, 
we can find A, k > such that 

47r2 max|iy(^)(/t)| sup / |/°(r(T)| r e^^^^ c/r < 1 - k. 

|o-|=l Jo 

Performing the change of variables kt/L = ra inside the integral, we deduce 



A7r^\W^^\k)\ 



\k\H 
L2 



and this obviously implies (L). 

The choice u; = in (12.81) shows that Condition (b) is a particular case of (c), so 
we only treat the latter assumption. The reasoning is more subtle than for case (a). 
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First we note that 
k) 



Att^ W{k) 



-Air 



2lTT- 



\k\^W{k)t 
k\W{k) 





\k\' 




L2 










/2i7r \k\t 


r 


V L 





t dv 
t dv 



L 



^'f,{v) e-^'^-t*^ dv. 



Then, for ,^ = 7 + ztu, using the formula 



e-'' e'"^' dt 



we get from (12.41) 



C{i,k) = W{k) / ^'^{v) 



s + iuj 



[v + uj) 



27 



(if. 



[v + ci;)2 + 72 

(To be rigorous, one may first estabhsh this formula for 7 < 0, and then use analyt- 
icity to derive it for 7 e [0, A).) 

As 7 — >■ 0, this expression approaches, uniformly in k and u;. 



(3.16) 



C{iu,k) = W{k) 



W{k) p.v. 



V + u + iO 



■ dv 



V + LO 



dv -i'KW{k)^'k{-uj) 



(Plemelj formula for the Cauchy transform). The problem is to show that C{iuj, k) 
stays away from 1 as varies in M; then the same will be true for ^ = 7 + ic^; with 
7 small enough. Equation (I3.16P shows that the imaginary part of C{iuj, k) vanishes 
only in the limit W{k) — (but then also the real part approaches 0), or in the 
limit \uj\ 00 (but then also the real part approaches 0), or if (p'f,{—uj) = 0; but 
then by (EBD 



C{iio,k) = W{k) 



V + u 



dv < 1, 



so even in this case C cannot approach 1. Case (c) of Proposition 12.11 follows by a 
compactness argument. □ 
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4. Analytic norms 

In this section we introduce some functional spaces of analytic functions on M*^, 
-jd _ jjd.^ most importantly T*^ x W^. (Changing the sidelength of the torus 
only results in some changes in the constants.) Then we establish a number of 
functional inequalities which will be crucial in the subsequent analysis. At the end 
of this section we shall reformulate the linear study in this new setting. 

Throughout the whole section d is a. positive integer. Working with analytic 
functions will force us to be careful with combinatorial issues, and proofs will at 
times involve summation over many indices. 

4.1. Single- variable analytic norms. Here "single- variable" means that the vari- 
able lives in either or T*^, but d may be greater than 1. 

Among many possible families of norms for analytic functions, two will be of 
particular interest for us; they will be denoted by C^'^ and J^^'^. The C^'^ norms 
are defined for functions on or T'^, while the J^^'^ norms are defined only for 
T'^ (although we could easily cook up a variant in W^). We shall write Nq for the 
set of (i-tuples of integers (the subscript being here to insist that is allowed). If 
n e Nq and A > we shall write A*^ = A'*^'. Conventions about Fourier transform 
and multidimensional differential calculus are gathered in the Appendix. 

Definition 4.1 (One-variable analytic norms). For any p G [1,cxd] and A > 0, we 
define 

A- / . 

(4.1) ii/iic.. := E 7^ wf^'^^w^^-^ yy^- ■= E 1/(^)1' 

the latter expression standing for sup^.(e^'^''*''^' |/(A;)|) if p = oo. We further write 

(4.2) C^'°° = C^, J^^'^ = J^^. 

Remark 4.2. The parameter A can be interpreted as a radius of convergence. 

Remark 4.3. The norms C'^ and are of particular interest because they are 
algebra norms. 

We shall sometimes abbreviate || ■ Wc^w or || ■ \\jr\;p into || • ||a;p when no confusion 
is possible, or when the statement works for either. 

The norms in (14. ip extend to vector-valued functions in a natural way: if / is 
valued in M'^ or T'^ or define = (/;"\ . . . , fP), f{k) = {f,{k),..., Uk))- 
then the formulas in (14. ID make sense provided that we choose a norm on or T'^ 
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or . Which norm we choose will depend on the context; the choice will always be 
done in such a way to get the duality right in the inequality |a ■ 6| < ||a|| For 
instance if / is valued in Z"^ and g in T*^, and we have to estimate / ■ (?, we may norm 
if" by \k\ = T'^ by \x\ = sup This will not pose any problem, and 

the reader can forget about this issue; we shall just make remarks about it whenever 
needed. For the rest of this section, we shall focus on scalar-valued functions for 
simplicity of exposition. 

Next, we define "homogeneous" analytic seminorms by removing the zero-order 
term. We write = M;^ \ {0}, = Z'^ \ {0}. 

Definition 4.4 (One- variable homogeneous analytic seminorms). For p G [1, C)o] 
and A > we write 

\n I ^ 

It is interesting to note that affine functions x ^ a ■ x + b can be included in 
_ (jA;oo^ even though they are unbounded; in particular ||a-a; + 6||(jA = A|a|. On 
the other hand, linear forms x i — > a ■ x do not naturally belong to JF'*', because their 
Fourier expansion is not even summable (it decays like 1/k). 

The spaces C^'^ and JF'*''?' enjoy remarkable properties, summarized in Propositions 
14.51 14.81 and 14.101 below. Some of these properties are well-known, other not so. 

Proposition 4.5 (Algebra property), (i) For any A > 0, andp,q,r G [1, +cxd] such 
that 1/p + 1/q = 1/r, we have 

\\f 9\\c^'- < WfWc^'P hWc^-^i- 
(ii) For any A > 0, and p,q,r E [1, +oo] such that 1/p + 1/q = 1/r + 1, we have 

||/^||^A;r < ||/||^A;p ||^||^A;,. 

(Hi) As a consequence, for any A > 0, C'^ = C'^'°° and = J^^'^ are normed 
algebras: for either space, 

\\f9\\x<\\f\\x\\9\\x- 
In particular, \\f^\\\ < ||/||a fo'r any n G Nq, and \\e^\\\ < e"-'^"^. 




Of course all norms arc equivalent, still the choice is not innocent when the estimates are 
iterated infinitely many times; an advantage of the supremum norm on is that it has the 
algebra property. 
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Remark 4.6. Ultimately, property (iii) relies on the fact that and are normed 
algebras for the multiplication and convolution, respectively. 

Remark 4.7. It follows from the Fourier inversion formula and Proposition 14.51 that 
ll/llc^ ^ II/IIj^^ (and \\f\\cx < II/IIjfa); this is a special case of Proposition HTSl (iv) 



below. The reverse inequality does not hold, because ||/||oo does not control 

Analytic norms are very sensitive to composition; think that if a > then 
11/ o (ald)llcAip = a~'^^P ll/llc^ip; so we typically lose on the functional space. This is 
a major difference with more traditional norms used in partial differential equations 
theory, such as Holder or Sobolev norms, for which composition may affect constants 
but not regularity indices. The next proposition controls the loss of regularity im- 
plied by composition. 

Proposition 4.8 (Composition inequality), (i) For any A > and any p G [1, +oo] , 

-i||i/p 

1 1 oo 



< lUdetWHY 



\\foH\\cX:j, 

where H is possibly unbounded. 

(a) For any A > 0, any p G [1, oo] and any a > 0, 



\H\ 



/o(aId +G) 
(iii) For any A > 0, 

/o(Id +G) 



< a 



'd/p 



V 



\\G\V 



< 



(iv) For any A > and any a > 0, 



/ o (a Id + G) 



< 



IGI 



JC-A- 



\G\ 



Remark 4.9. Inequality (iv), with C on the left and T on the right, will be most 
useful. The reverse inequality is not likely to hold, in view of Remark 14.71 

The last property of interest for us is the control of the loss of regularity involved 
by differentiation. 

Proposition 4.10 (Control of gradients). For any X > X, any p G [1, +oo], we have 

1 
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The proofs of Propositions 14.51 to 14.101 will be preparations for the more compli- 
cated situations considered in the sequel. 

Proof of Proposition \4 ■ 5\ (i) Denoting by || ■ ||a;p the norm of C^'^, using the multi- 
dimensional Leibniz formula from Appendix IA.2t we have 



feNj; rn<e 

£EE(™)ii/""'iu.iif""-iu4 

ml U-m)\ 



(ii) Denoting now by || ■ ||a;p the norm of JF^'*', and applying Young's convolution 
inequality, we get 



\ k £ 

J2 imi^e^^^^^'-A ' (j2 i^wre'"^"'" 



l/r 



X\£\ 



< 



□ 

Proof of Proposition 4-8. Case (i). We use the (multi-dimensional) Faa di Bruno 
formula: 



(/oi7)(«)= 



mi! . . . m„! 



^j(mi+...+m„) o if) JJ 



so 



n! 



mi! . . . m^ 



II " H^i'' 

y(mi+...+m„) Q ^ TT 
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thus 

n>l \ k=l 



y ^ n 

mi! . . . mJ. 

j:]^Jm,=n, j:]^,7n,=k " j=l 



/ 

ll(detVif)-|i::^ 

V 



fc>i ■ \\e\>i ' j 



where the last step follows from the multidimensional binomial formula. 
Case (ii). We decompose h{x) := f{ax + G{x)) as 

d 



and we apply V*^: 

= 

ki\k2 \ n\ 



fci+fc2=A:,eNg nGNg 

Then we take the norm, multiply by X^/k\ and sum over k: 

\ki+k2 |„|fci 

ki,k2,n>0 



fci,n>0 



Klin! 



k-i ! n! 

fci,n>0 



ml 

where Proposition 14.51 (iii) was used in the but-to-last step. 
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SO 



Case (iii). In this case we write, with Go = G{0), 

h{x) = fix + Gix)) = J2 e''"'"" e'*"'^-^" ^2inkiG{x)~Go). 

k 



Then (using again Proposition 14.51) 



k i 



^2iwk-{G-Go) 



< 



< 



k 

E |7(A;)| e^'^^l^l gll27rfc.(G-Go)IU 

k 

|J(^^^|g27rA|fc|g27r|fc|||G-Go|U 

k 

A+l|G-Go|U = II/IIa+!^, ^ = \\G\\jrx. 



Case (iv). We actually have the more precise result 

(4.5) ||/oif||c. <$^|/(A;)|e2-l'=lll^^llc.. 

Writing f o H = ^ f{k) e'^'^'"^'^ ^ we see that 04.51) follows from 

(4.6) lle^'^IIcA < eli^llcA. 

To prove (14. 6p . let P„ be the polynomial in the variables Xm {fn < n) defined by the 
identity (c'^)^"-' = -Pn((/*'™'^)m<n) e-^; this polynomial (which can be made more ex- 
plicit from the Faa di Bruno formula) has nonnegative coefficients, so || (e*-^)*^"'-' ||oo < 
-Pn((||/^™^||)m<n)- The couclusiou will follow from the identity (between formal se- 
ries!) 

(4.7) ^^^^ P„((X^)„,<„) = exp 5^ — X J . 

To prove (14.71) . it is sufficient to note that the left-hand side is the expansion of 
in powers of A at 0, where g{\) = J2keNi fcT -^fc- ^ 
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Proof of Proposition \4-10[ (a) Writing 



|A;p 

A" 



we have 



n 

where di = d/dxi. If Ij is the d-uple of integers with 1 in position i, then (n + Ij)! < 
(|n| + l)n!, so 



and the proof of (14.31) follows easily, 
(b) Writing || ■ ||a;p = || ■ \\j^^;p, we have 



|m|>l 



1/p 



p g27rAp |fc| 



< 



sup {\k\e 



27r(A-A) \k\ 



1/p 



and (144|) follows. 



□ 



4.2. Analytic norms in two variables. To estimate solutions and trajectories of 
kinetic equations we will work on the phase space x M^, and use three parameters: 
A {gliding analytic regularity); fi (analytic regularity in x); and r (time-shift along 
the free transport semigroup). The regularity quantified by A is said to be gliding 
because for r = this is an analytic regularity in v, but as r grows the regularity is 
progressively transferred from velocity to spatial modes, according to the evolution 
by free transport. This catch is crucial to our analysis: indeed, the solution of a 
transport equation like free transport or Vlasov cannot be uniformly analytic^ in v 
as time goes by — except of course if it is spatially homogeneous. Instead, the best 
we can do is compare the solution at time r to the solution of free transport at the 
same time — a kind of scattering point of view. 

The parameters A, fi will be nonnegative; r will vary in M, but often be restricted 
to M+, just because we shall work in positive time. When r is not specified, this 
means r = 0. Sometimes we shall abuse notation by writing ||/(x, t")]! instead of 
to stress the dependence of / on the two variables. 



""^^By this we mean of course that some norm or seminorm quantifying the degree of analytic 
smoothness in v wiU remain uniformly bounded. 
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Putting aside the time-shift for a moment, we may generahze the norms and 
J-"^ in an obvious way: 

Definition 4.11 (Two-variables analytic norms). For any A,/i > 0, we define 
(4.8) A /i 



^ ^ n! m! ^ 



(4.9) 



Of course one might also introduce variants based on or norms (with two 
additional parameters p, q, since one can make different choices for the space and 
velocity variables). 

The norm (14. 9 p is better adapted to the periodic nature of the problem, and is very 
well suited to estimate solutions of kinetic equations (with fast decay as |f | oo); 
but in the sequel we shall also have to estimate characteristics (trajectories) which 
are unbounded functions of v. We could hope to play with two different families 
of norms, but this would entail considerable technical difficulties. Instead, we shall 
mix the two recipes to get the following hybrid norms: 

Definition 4.12 (Hybrid analytic norms). For any A,yU > 0, let 



(4.10) 



A" 



EE 

n! 



V^/(£,tO 



More generally, for any p G [1, oo] we define 

\ 

(4.11) WfWzK.. ^ ^ 



V«/(£,t;) 



LP 



Now let us introduce the time-shift r. We denote by {S^)r>o the geodesic semi- 
group: (S'°)(a:, v) = {x + VT, v). Recall that the backward free transport semigroup 
is defined by (/ o S^)r>o, and the forward semigroup by (/ o 



!t>Q- 



Definition 4.13 (Time-shift pure and hybrid analytic norms) 



(4.12) 



EE 



A" fi' 



V^(V. + rV.)"/ 
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(4.13) 
(4.14) 



i. = \\foSl\ 



\f%ri)\e 



27rA|fcT+r?| 27r^|fc| 



A" 



Zr" 



^ — ^ ^ — ' n! 



(4.15) 



EE 



LP 



This choice of norms is one of the cornerstones of our analysis: first, because of 
their hybrid nature, they will connect well to both periodic (in x) estimates on the 
force field, and uniform (in v) estimates on the "scattering transforms" studied in 
Section [51 Secondly, they are well-behaved with respect to the properties of free 
transport, allowing to keep track of the initial time without needing ridiculous (and 
inaccessible) amounts of regularity in x as time goes by. Thirdly, they will satisfy the 
algebra property (for p = oo), the composition inequality and the gradient inequality 
(for any p G [1, oo]). Before going on with the proof of these properties, we note the 
following alternative representations. 

Proposition 4.14. The norm Z^'^^'^ admits the alternative representations : 



(4.16) 

(4.17) 

where 
(4.18) 



A" 



LP 



2^'t^'P 



E 



A^ 



(V, + rV,)"/ 



I LP 



4.3. Relations between functional spaces. The next propositions are easily 
checked. 



Proposition 4.15. With the notation from Subsection 4-3. for any r G M, 
(i) if f is a function only of x then 



ll/l 



CA|t|+m, 
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(ii) if f is a function only of v then 



(Hi) for any function f = f{x, v), if {■) stands for spatial average then 

||(/)||c^^p < ||/||^A,m;p; 

(iv) for any function f = /(x, v), 



fdv 



■fX\T\+IJ. 



< II / I|^a,m;1. 



Remark 4.16. Note, in Proposition 14.151 (i) and (iv), how the regularity in x is 
improved by the time-shift. 



Proof of Proposition \4-l3[ Only (iv) requires some explanations. Let p(a;) = J f{x,v)dv. 
Then for any /c G Z'^, 



so for any n G Nq, 



Pik) = / f{k,v)dv; 



{2i-KtkY p{k) = {2i-KtkYf{k,v)dv 



(V„ + 2nxtk)'' f{k, v) dv. 



Recalling the conventions from Appendix lA.ll we deduce 

.1 |27rAtA;r 



k,n 



,27r/x|fc 



\p{k)\ < Ve^'^^l'^l ^ [ (V, + 2i'KtkYf{k 
^-^ n\ J 

k,n 



dv 



?A,fi;l . 



□ 



Proposition 4.17. With the notation from Subsection \4-2 

A < A', < /i' =^ 



Moreover, for r, r G M, and any p G [1, oo], 



(4.19) 



^A,m;p < 11/ II yA,M+A|--T|;p. 

T T 
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Remark 4.18. Note carefully that the spaces Z^'^ are not ordered with respect 
to the parameter r, which cannot be thought of as a regularity index. We could 
dispend with this parameter if we were working in time 0(1); but (14.191) is of course 
of absolutely no use. This means that errors on the exponent r should remain 
somehow small, in order to be controllable by small losses on the exponent /i. 

Finally we state an easy proposition which follows from the time- invar iance of the 
free transport equation: 

Proposition 4.19. For any X G {C, JF, Z}, and any t, r G M, 

ll/°5?llx^' = ll/llx,Y/ 

Now we shall see that the hybrid norms, and certain variants thereof, enjoy prop- 
erties rather similar to those of the single-variable analytic norms studied before. 
This will be sometimes technical, and the reader who would like to reconnect to 
physical problems is advised to go directly to Subsection I4.11[ 

4.4. Injections. In this section we relate norms to more standard norms 

entirely based on Fourier space. In the next theorem we write 

(4.20) = = sup sup e2-H^^le2-A|.+fcr| \J{^k,7^)\_ 

Theorem 4.20 (Injections between analytic spaces), (i) If X,fi > and r G M then 

(4.21) WfWyy < 



(u) If < X < X, < fi< Jl < M , T e M, then 

(in) If < X < X < A, < j2 < jL < M, b < (3 < B, then there is C 
C(A, M, b, B, d) such that 



^.,,.1 <C-MX-A;7.-M} \\f\\^j-,+ 



max 



|/(x,w)|e^l^ldt;dx) ; ( / / \f{x,v)\e'^^''^dvdx^^ 
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Remark 4.21. The combination of (ii) and (iii), plus elementary Lebesgue inter- 
polation, enables to control all norms ^^'^-p, 1 < p < oo. 

Proof of Theorem \4.20[ By the invariance under the action of free transport, it is 
sufficient to do the proof for r = 0. 

By integration by parts in the Fourier transform formula, we have 

f{k, v) 6-2-"- dv = j V^fik, V) ^—^ dv. 

So 

and therefore 



e^-^H \f{k,v)\<e''''^'\Y.-^\vnf{k,v)\ 

n 

^e'"''"' E^y \^if{k,v)\dv. 



This establishes (i). 

Next, by differentiating the identity 



f[Kv)= I f{k,r])e'^-^-'^dv, 

we get 

(4.23) V™/(A;, v) = / /{k, 7]) (2i7rr/)™ e^^'^'''^ dr]. 

Then we deduce (ii) by writing 

J2e'-^\'\ — \\V-f{k,v)h^^,^) 
^-^ ml 

k,m 



< 

k 



k 

- fj^^e"'"^''"''^!''' j (y e-'"(^-^)l''l rfr/^ (^sup e^-^l*?! e^^'^l'^l |7(A;, r/)| 
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The proof of (iii) is the most tricky. We start again from (14.231) . but now we 
integrate by parts in the rj variable: 



V^[/(A:,r/) (2z7rr/)™] -— — rf^, 

where q = q{v) is a multi-index to be chosen. 

We split into 2^ disjoint regions A(zi, . . . , where the ij are distinct indices 

in 

A(/) = jt; e M^; \vi\ > 1 Vi G /, \vi\ < 1 Vi ^ /}. 

If V G A(ii,...,i„) we apply (14.241) with the multi-index q defined by qj = 2 if 
j G {ii, . . . , in}, Q'j = otherwise. This gives 

/ \V-Kk,v)\ dv < f-i^ / sup V^[7(fc,r/) (2^vrr/)-] 

Summing up all pieces and using the Leibniz formula, we get 

/ \V':;:f{k,v)\dv<C{d){l + m^'') sup sup \S/lf{k,7])\\2Ti7]r-'i. 

J k,ri \q\<2d 

At this point we apply Lemma 14.221 below with 

1 . ( X — \ u — u 
e = - mm < — = — : 

4 \ \ -p 

and we get, for q < 2d, 
I V^7(/t, v)\ < C(rf)"''"{^ • ^} K{b, B) 



(sup e \^\\f{k, r^) I j max | |^sup j ; |^sup 

Of course, 

— 7^ — <(w) y^j/(^,.)i^c^. 

</ |/(a;,^)|(27r/5)^^rft;< / \fix,v)\ e'^^^^^ dv. 

. had V. ./md 



2e 

oo 
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So, all in all, 
^-^ ml J 



< Y C(ci,A,M,6,S)™{^-^^T^-"} 

\q\<2d 



' m\ 

l—e { / r- \ £ 



Since 

V^'"'^ + '">"l^"'l""<Cto.A)e^-^l-l 
^-^ ml 

m 

and 

fc k 

we easily end up with the desired result. □ 

Lemma 4.22. Let f -.W^ ^ C, and let a > Q, A > I, q e Let (3 such that 
< b < P < B . If \f{x)\ < y4e~"l^l for all x, then for any e G (0, 1/4) one has 

(4.25) I VV(a;) I < df^ K{b, B) A^-' e-^^"'")"!"! 



I ( or 11^ / 1 1 oo A / 1 1 V / 1 1 oo 

sup max < p — ; P 



7*! / V T! 



Remark 4.23. One may conjecture that the optimal constant in the right-hand side 
of fl4.25p is in fact polynomial in if this conjecture holds true, then the constants 
in Theorem 14.201 (iii) can be improved accordingly. Mironescu communicated to us 
a derivation of polynomial bounds for the optimal constant in the related inequality 

II fWii , ^ < r(iA II f i|i/(fc+2)|| .(fc+i)||('=+i)/(fc+2) 

based on a real interpolation method. 
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Proof of Lemma \4 ■ 22[ Let us first see / as a function of xi, and treat x' = (x2, . . . , Xd) 
as a parameter. Thus tlie assumption is < (Ae""'^ ') e~"I^^L By a more 

or less standard interpolation inequality [21] Lemma A.l], 



(4.26) < 2VAe-'^\-Ue-^\'-^\\\dlf{xi,x')\\i = 2V Ae-^\'-\^J\\^ff\\oo. 

Let Cq-^^^ri be the optimal constant (not smaller than 1) such that 



(4.27) \drfix,,x')\ < C,,,r, (Ae""l^l)^" ||9;7(a:i, x')||5 • 

By iterating (I4.26p . we find C^^^j-i < '^^/Cq^^i^n Cgi+i,ri- It follows by induction that 

Next, using (14.271) and interpolating according to the second variable X2 as in 
f|T26D . we get 

We repeat this until we get 

(4.28) |VV(x)| < {C,,,r. . . . C,„,J (Ae""l^l)(^"^)-(^-^) 

Choose (1 < i < (i) in such a way that 

e Qi 2e 
-< — <—; 

d ri d 

this is always possible for e < d/A. Then Cg^^r, < (2'^''?)^/^, and f l4.28p imphes 

|VV(x)|<(2'^l«l^)^/M^e-"l^'l)i- max{||V7roc; l|V7ll^} • 

Then, since 2(r + g)£: < 3dq we have, by a crude application of Stirling's formula (in 
quantitative form), for s < r + q, 

'/3N|VVI|oo\' fs\ 



l|VVIIL< 



e 



< sup j C{p, q, d) e 



■3dq 
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and the result follows easily. 



□ 



4.5. Algebra property in two variables. In this section we only consider the 
norms Z^''^'^; but similar results would hold true for the two-variables C and JF 
spaces, and could be proven with the same method as those used for the one- variable 
spaces J^^ and respectively (note that the Leibniz formula still applies because 
Vx and (V^ -|- rV^) commute). 

Proposition 4.24. (i) For any A,/i > 0, r G M and p,q,r G [1, +oo] such that 
l/p+l/q = l/r,we have 



\\f g\\^x,t,:r < ll/II^A.Mip lls'll^^'" 



(a) As a consequence, Z^'^ = Z^''^'°° is a normed algeb 



ra: 



WfgW 



In particular, \\f"\\^x,^, < 



for any n G No, and We-^W^^x.t^ < e ^^'^ . 



Proof of Proposition 4-24' First we note that (with the notation (14.181) ) ||| ■ \\\n-r sat- 
isfies the ''{p,q,r) property": whenever p, g, r G [1,+cxd] satisfy 1/p + 1/g = 1/r, we 
have 



£e2 



Y^f(k,-)g{i-k, 



iim;p I US' III a*;?- 
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Next, we write 



A" 



A" 



n m<n 



^ E ^ E ( ^ ) + III(V. + Tv.r~'-g\ 



n m<n 



\ m 



A^ 



I ^A,/j,;p |||5'||| ^■^'f^'il ' 

(We could also reduce to r = by means of Proposition I4.19[ ) 



□ 



4.6. Composition inequality. 

Proposition 4.25 (Composition inequality in two variables). For any X,n> and 
any p e[l,oo], r eR, a eR, a eR \ {0}, beR, 



(4.29) 

where 



f{x + bv + X{x, v), av + V{x, v)) 



(4.30) 



a 



A a + \\V\\ j.A,^, 



P = ^ + \ \b + T — aa\ + \\X — aV\\^\,p 



Remark 4.26. The norms in (14.301) for X and V have to be based on L°°, not just 
any L^. Also note: the fact that the second argument of / has the form av + V (and 
not av + cx + V) is related to Remark 14.71 



Proof of Proposition \4 ■ 25\ The proof is a combination of the arguments in Proposi- 
tion In a first step, we do it for the case r = o" = 0, and we write || ■ ||a,/x;p = 
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From the expansion /(x, v) = Y^ f{k, v) e^*'^^'^ we deduce 



h{x, v) := f{x + hv + f ), av + V{x, f ) j 
= f{k, av + V) e2^-fc-(-+''-+^) 

k 

ym 



. ^ 2i7rk-x 2i7rk-bv 2i7rk-X 

LL U } ' O O O • 



k m 



Taking the Fourier transform in x, we see that for any i G Z*^, 



h{e, v) = Y,Y. ^™/(^' ^'"'■'^ E ^^^^^ (e^^-'^-^Kf -k-j). 

km j ' 



Differentiating n times via the Leibniz formula (here apphed to a product of four 
functions), we get 



rii! rio! ?T-3' ?T-4' 

k,m.,j ni+n2+ns+n4=n 
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Multiplying by A" e^'^'^'^' / n\ and summing over n and i, taking norms and using 
< ||/||lp||5'IU-, we finally obtain 



< \a\ 



E 



fcj/£Zg; m,n, ni+n2+n3+n4=n>0 



ni\ 77-2! n^}. 



jm+n-i 



f{k, 



I LP 



-d/p 



V^He'*"'^-^)(£- A;-j) (27r|6| |A;|)" 

^ni+n2+n3+n4 g27r/i|fc| g27r/i|j| g27r/i|£-fc-j| IqI"! 



fcjjfGZJj, m,ni,n2,n3,n4>0 



rii! 77-2! 77-3! 71.4! 



iivr"v(^,-)iiLP 



E\ni I „|ni 



k,ni ,m 



rill 



ml 



\LP 



\n3 



^27r/i|fc| 



00 

1 



"2 J 



E^^"'""'iiK'(e"'"-"r(ft)iL) (E 



(2jrA|()||A|)»' 



V rig, ft, 



n4 



774! 



-<i/p ^ (A|a|)"i g27rAt|fc| llym+mj^^ )||lp ^^^""^^^'^ [[g2»7rfc-X||^ ^ g27rA|fe||fc| 



k,p,m 



Till 



<\a\-'^/P (^MZe27r(M+A|6|)|fc| ||yni+mj('^^ All H H A,M g27r|fc| ||X|U. 



ml 



k,ni ,m 



ni 



«I"'^"E^(^i«i + II^IU'/^)"II^"/(^'-)IIl.^^ 

k,n 



ml 



27r|fc|(At+A|fe| + ||X||A,^) 



-d/p 



A|a| + ||V||A,M.M+Alfe| + ||X|U,^ 
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Now we generalize this to arbitrary values of a and r: by Proposition 14.191 

/( X + vi}) + r) + X{x + VT, v), av + V(x + vt, v) ) 
f o S'^o S°^^(x + v{b + T) + X{x + vT,v), av + V{x + VT,v)) 
{f o S^)(^x + v{b + T - aa) + {X - aV){x + vt,v), av + V{x + vt,v) 
(/ o 5°) (^x + v{b + T-aa) + Y{x, v), av + W{x, v)^ 



where 



W{x, v) = Vo S^{x, v), Y{x, v) = {X - aV) o S^{x, v). 

Applying the result for r = 0, we deduce that the norm of h{x,v) = f{x + bv + 
X{x,v), av + V{x,v)) in Z^'f^ is bounded by 



where 



and 



a = X\a\ + \\VoS^ 



^||^A,M = |a|A + ||\^||^A,M, 



P = fi + \\b + T - aa\ + \\{X - aV) o S°\\z>~,^^ = fj. + \\b + r - aa\ + \\X - aV\ 
This establishes the desired bound. 

4.7. Gradient inequality. In the next proposition we shall write 



□ 



(4.31) 



^ ^ n! 



,2n^^\e\ 



{W, + 2i7rTirfii,v] 



This is again a homogeneous (in the x variable) seminorm. 

Proposition 4.27. For X > X > 0, JI > fi > 0, we have the functional inequalities 



C{d) 



(/i-/i) 

<^^llfll - 

- Alog(A/A) ll^^-- 
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In particular, for r > we have 



Alog(A/A)y \{fi-fi)J 

The proof is similar to the proof of Proposition I4.10( the constant Cid) arises in 
the choice of norm on M*^. As a consequence, if 1 < A/ A < 2, we have e.g. the bound 



l|v/iu,„<cw(^ + i±l) 



I 2^,m;p ■ 

4.8. Inversion. From the composition inequahty follows an inversion estimate. 

Proposition 4.28 (Inversion inequality), (i) Let X, jj, > 0, r G M, and F : T^xM.^ 
j^d ^ Then there is e = e{d) such that if F satisfies 

||V(F-Id)||^v,,/ <e{d), 



where 

A' = A + 2\\F - Id II^A,^, /x' = /X + 2(1 + |r|) ||F - Id 



then F is invertible and 



(4.32) - Id II^A.M < 2 ||F - Id || ^a,^. 



(a) More generally, if F and G are functions T'' x M'^ ^ T'' x M.'^ such that 

(4.33) ||V(F-Id)||^AV' <£(ci), 
where 

A' = A + 2\\F - G\\^x,^, fi' = 11 + 2(1 + \r\) \\F - G\\^x,^, 
then F is invertible and 

(4.34) o G - Id II^A,^ < 2 ||F - GH^a,,.. 

Remark 4.29. The conditions become very stringent as r becomes large: basically, 
F — Id (or F — G in case (ii)) should be of order o{1/t) for Proposition 14.281 to be 
applicable. 

Remark 4.30. By Proposition 14.271 a sufficient condition for fl4.33p to hold is that 
there be X",fi" such that A < A" < 2A, < /x", and 

IIF — Id II „a",m" < , ^ ^ min{A" — A' ; fi" — fi'}. 

1 + T 

However, this condition is in practice hard to fulfill. 
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Proof of Proposition \4 ■ 2^ We prove only (ii), of which (i) is a particular case. Let 
/ = F - Id , /i = F-i o G - Id , ^ = G - Id , so that Id + ^ = (Id + /) o (Id + /i), or 
equivalent ly 

h = g-fo{\d+h). 

So is a fixed point of 

$ : Z I — > g - f o{\d + Z). 

Note that $(0) = g- f. If $ is (l/2)-Lipschitz on the ball S(0,2||/ - ^||) in Z^^i", 
then (14.341) will follow by fixed point iteration as in Theorem IA.2I 
So let Z, Z be given with 



^11^^"' ll^ll^^''' < 2||/ - fi-II^A.M. 



We have 



$(Z) - ^{Z) = /(Id + Z) - /(Id + Z) 



v/(id + {i-e)z + ez) de] -{z-z). 



By Proposition 14.241 

i|$(Z)-$(Z)||^.,. < (^'||v/(ld+(l- 
For any 9 G [0, 1], by Proposition 14.251 

vf(id + {i-e)z + ez 

where 



)Z + ^Z) 



C?^ I \\Z — 2'||^A,n. 



<I|V/|L 



A = A + max{||Z||; ||Z||} < A + 2||/ - (7] 
and (writing Z = [Z^, Z^), Z = [Z^, Z^)) 

Ji = fi + max{||Zx. - tZi,\\ ; - tZ^\\} < + 2(1 + |r|) ||/ - g\\z^,^^. 
If F and G satisfy the assumptions of Proposition I4.28[ we deduce that 

||*||Lip(B(o,2)) < C{d)e{d), 
and this is bounded above by 1/2 if e{d) is small enough. □ 
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4.9. Sobolev corrections. We shall need to quantify Sobolev regularity corrections 
to the analytic regularity, in the x variable. 

Definition 4.31 (Hybrid analytic norms with Sobolev corrections). For A,/i, 7 > 0, 
r G M, p G [1, 00], we define 

\n ^ 

^.,. = Ee"^""(i + i^ir 1/(^)1- 

Proposition 4.32. Lei A,/i, 7 > 0, r G M anc? j» G [1,+cxd]. We have the following 
functional inequalities: 

('^) |l/ll;7>'.(f.7);P = ||/o5'°|| a,(m,7);p; 

("m^ 1/g = 1/r =^ \\fg\\-z^,(^^,i).r < ||/||^A,(M,^);p llfi'll^A.c^.^);, and therefore m 

particular Zr'^^''^^ = Zr'^^''^^'°° is a normed algebra; 

(Hi) If f depends only on x then ||/|La,(m,7) = ||/||jrA|r|+M,7; 

('i'V) ||/L^(M.7);P < ||/|L^,(M + A|r-T|,7);p; 

(v) for any a eR, a eR \ {0}, beR, pe[l, 00], 



f (^x + bv + X(x,v), av + V{x,v)j 



where a = X\a\ + \\V\\ „a,(m,7) and /5 = /i + A|6 + r — acr| + ||X — frV || _a,(h,7) . 
(^mj Gradient inequality: 

II ,A,(M,7);p < _ ^ ^ 11/11 yA,(7r,7);p, 

^ — ^ ^ 



||V/|La,(m,7);p < C(ci) 



1 1 + r 



^A,c.,7,P ^ -vu; ^^-^ ^ ^— ^ y ||^A,(jr,7) 

(^mj Inversion: If F and G are functions T'^ x M*^ — T'^ x Mf^ such that 

||V(F-Id)||^A',(.',7) <e(c^), 

where 

X' = X + 2\\F - G||^A,(.,7), /i' = /i + 2(1 + |r|) ||F - G||2a,(..7), 
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then 
(4.35) 



\F o G* — Id II _a,(m,7) ^ 2 ||i^ — Gil „a,{m,7)- 



Proof of Proposition 4-32. The proofs are the same as for the "plain" hybrid norms; 
the only notable point is that for the proof of (ii) we use, in addition to e^'^'^''^' < 
g27rA|fc-f| g27rA|f|^ ^Yie inequality 

(l + |A;|r< (1 + |A;-£|)^(1 + |£|)^. 

□ 



Remark 4.33. Of course, some of the estimates in Proposition 14.321 can be "im- 
proved" by taking advantage of 7; e.g. for 7 > 1 we have 

\\'^J\\zy:.<C{d) ||/||^.,(...);.. 

4.10. Individual mode estimates. To handle very singular cases, we shall at times 
need to estimate Fourier modes individually, rather than full norms. If / = f{x,v), 
we write 

(4.36) {P,f){x,v) = f{k,v)e''-'-'. 

In particular the following estimates will be useful. 

Proposition 4.34. For any A,/i > 0, r G M, Lebesgue exponents 1/r = l/p+ 1/q 
and k G Z*^, we have the estimate 



k~i9\\ z^''^''' 



Proposition 4.35. For any A > 0, /i > /i > 0, r G M, p G [I, 00] and k G Z'^, we 
have the estimate 



f{x + X{x,v),v)] ^^^^^ < ^ e-2-('^-'')l^-^l ||P,/||^., 



u = n+ \\x\ 



These estimates also have variants with Sobolev corrections. Note that when 
^ = Ji, Proposition 14.351 is a direct consequence of Proposition 14.251 with = 0, 
6 = and a = 1: 



Pk 



f{x+X{x,v),v) 



Z^'t^''P 



< 



f{x+X{x,v),v) 



< 



Z/ = 11^ II ^-^./J 



Proof of Propositions \4 ■ 34\ and \4-3^ The proof of Proposition 14.341 is quite similar 
to the proof of Proposition 14.241 (It is no restriction to choose r = because Pk 
commutes with the free transport semigroup.) Proposition 14.351 needs a few words 
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of explanation. As in the proof of Proposition 14. 251 we let h{x,v) = f{x + X{x,v),v), 
and readily obtain 



I Pkh\\ 



\n 2nfj,\k\ 



SEE 

At this stage we write 



^2nfj,\k\ 



A 



nl 



\m>0 



L°°(dv) 



and use the crude bound 



^27r7l|fc-^| 



< 



L°°{dv) 



The rest of the proof is as in Proposition 14.251 



□ 



4.11. Measuring solutions of kinetic equations in large time. As we already 
discussed, even for the simplest kinetic equation, namely free transport, we cannot 
hope to have uniform in time regularity estimates in the velocity variable: rather, 
because of filamentation, we may have ||Vt,/(t, ■)|| = ||V^/(t, ■)|| = O(t^), etc. 

For analytic norms we may at best hope for an exponential growth. 

But the invariance of the "gliding" norms Z^'^ under free transport (Proposition 
14.191) makes it possible to look for uniform estimates such as 



(4.37) 



ll/(r,-)|| 



X.fi 



0(1) as r ^ +00. 



Of course, by Proposition 14.271 fl4.37p implies 
(4.38) ||V,/(r, Oll^v,.' = 0(r), A' < A, /i' < /i, 

and nothing better as far as the asymptotic behavior of Vt,/ is concerned; but fl4.37p 
is much more precise than (14.381) . For instance it implies ||(V^ + rVa;)/(T, ■)||^a',m' = 
0(1) for A' < A, /i' < /i. 

Another way to get rid of filamentation is to average over the spatial variable x, 
a common sense procedure which has already been used in physics |^ Section 49]. 
Think that, if / evolves according to free transport, or even according to the lin- 
earized Vlasov equation (13.31) . then its space-average 



(4.39) 



{f)ir,v) 



f{T,x,v) dx 
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is time-invariant. (We used this infinite number of conservation laws to determine 
tlie long-time behavior in Theorem 13.11 ) 

The bound fl4.37p easily implies a bound on the space average: indeed, 

(4.40) ||(/)(r,-)||c^ = ||(/)(r,-)b.,M < ||/(r, O^.,. = 0(1) as r ^ oo; 
and in particular, for A' < A, 

(4.41) ||(V,/)(r,-)||cA' =0(1) asr^oo. 
Again, (14.371) contains a lot more information than (14.411) . 

Remark 4.36. The idea to estimate solutions of a nonlinear equation by comparison 
to some unperturbed (reversible) linear dynamics is already present in the definition 
of Bourgain spaces X"^'^ [12]. The analogy stops here, since time is a dummy variable 
in X^'^ spaces, while in spaces it is frozen and appears as a parameter, on which 
we shall play later. 

4.12. Linear damping revisited. As a simple illustration of the functional analy- 
sis introduced in this section, let us recast the linear damping (Theorem 13. ip in this 
language. This will be the first step for the study of the nonlinear damping. For 
simplicity we set L = 1. 

Theorem 4.37 (Linear Landau damping again). Let /° = f^{v), : T'^ — s> M such 
that llVl^ll^i < Cw, o-nd fi{x,v) such that 

(i) Condition (L) from Subsection \2.S\ holds for some constants Co, A, k > 0; 

(^^) lirilc^a <Oo; 

(Hi) ||/i||^A,M;i < S for some fi > 0, 6 > 0; 

Then for any A' < A and p' < p, the solution of the linearized Vlasov equation (13. 3p 
satisfies 

(4.42) sup \\f{t,-)\\ ,>y.,,<CS, 

for some constant C = C{d, Cw, Co, A, A', ) . In particular, p = J f dv satisfies 

(4.43) sup ||p(t, ■)|Uvi,+.' <0(5. 

As a consequence, as \t\ — > oo, p converges strongly to poo = // fi{x,v) dx dv , and 
f converges weakly to (/j) = J fidx, at rate 0(e~'*'"'*') for any \" < A'. 
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If moreover \\f^\\c>^-p < C*o (in-d ||/i||2^,M;p < S for all p in some interval [l,p], then 
(14.421) can be reinforced into 



(4.44) 



sup ||/(t, 



< C6, 



I < p < p. 



Remark 4.38. The notions of weak and strong convergence are the same as those 
in Theorem 13. 1[ With respect to that statement, we have added an extra analyticity 
assumption in the x variable; in this hnear context this is an overkill (as the proof 
will show), but later in the nonlinear context this will be important. 



Proof of Theorem \4.31\ Without loss of generality we restrict our attention to t > 0. 
Although fj03|) follows from fj02l) by Proposition KT5[ we shall establish 
first, and deduce (14.421) thanks to the equation. We shall write C for various con- 
stants depending only on the parameters in the statement of the theorem. 
As in the proof of Theorem 13.11 we have 

p{t,k) = fi{k,kt) + [ K'^{t-T,k)p{T,k)dT 
Jo 

for any t > 0, A; G Z'^. By Lemma [3.61 ^ot any A' < X, p' < fi, 



-27^{t^-^,')\k\ 



sup(5^|p(t,fc)|e^-(^'*+-')l'=l 
<CiX,X',K) (j^e 

V k 

7^ sup [y2\fi{k,kt) 



C(A, A', k) 



sup sup \fi{k, kt) I e 

^27r(At+A»)|fc| 



27r(A't+^)lfc| 



t>0 



(/i-/i 

Equivalently, 

(4.45) sup ||p(t, OLaw < C sup 

By Propositions 14.151 and 14.191 

j fioS'Ltdv 
This and flOSD imply flCTj) 



fi o S\ dv 



< \\fi O S'°i||^A,M;l = WftW^X.^;! < 6. 
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To deduce (14.421) . we first write 



fit, 



where pr = p{t, ■ ). Then for any A" < A' we have, by Propositions 14.24 1 and 14.151 
for all t > 0, 

(4.46) 

" II V.r II w,,,cir 









Jo 



(/ ll'^^*P-^ll.^^"-+'^' l|V^/°llc^'M- 



Since VH^(O) = 0, we have, for any r > 0, 

II VW^ * P.||^vv+. < e~^-^^"-^> \\VW * p.||^v.+.' 

<||VW^|U.e-2-(^"-^>||p.||^v.+.' 

in particular 

(4.47) ri|VW^*PrLv'+.' < 







Also, by Proposition I4.10[ for 1 < A'/A" < 2 we have 

C ,, „n,, . C Cq 



(4.48) llV.rilc^". < ^^-^ lir lie.. - A - A"- 

Plugging (14.471) and (14.481) in (14.461) . we deduce (14.421) . The end of the proof is an 
easy exercise if one recalls that {fit, • )) = {fi) for all t. □ 

5. Scattering estimates 

Let be given a small time-dependent force field, denoted by e F(t, x), on T*^ x Mf^, 
whose analytic regularity improves linearly in time. (Think of eF as the force created 
by a damped density.) This force field perturbs the trajectories S^^ of the free 
transport (r the initial time, t the current time) into trajectories Sr,t- The goal of this 
section is to get an estimate on the maps flt,T = St^r ° (so that St^r = ^t,T ° S^^). 
These bounds should be in an analytic class about as good as F, with a loss of 
analyticity depending on e; they should also be (for < r < t) 

• uniform in t > r; 
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• small as r ^ oo; 

• small as r — t. 

We shall informally say that Qt.r is a scattering transform, even though this 
terminology is usually reserved for the asymptotic regime t — > ±00. 

Remark 5.1. The order of composition of the free semigroup and perturbed semi- 
group is dictated by the need to get uniformity as t — > 00. If we had defined, say, 
At^r = S^t ° St^T, SO that St^T = S^^ o At,,-, and if the force was, say, supported in 
< t < 1, we would get (denoting St^r = {^t,T, ^i,r)) 



At,o{x, v) = (^Xi,o(x - v{t - l),v) + tVi,o{x - v{t -l),v), Vi,o{x - v{t - 1), v) 
which does not converge to anything as t — > 00. 

5.1. Formal expansion. Before stating the main result, we sketch a heuristic per- 
turbation study. Let us write a formal expansion of Vo,j(x, as a perturbation 
series: 

Vo,t{x, v) =v + e v^^\t, X, v) + v^'^\t, x,v) + ... 

Then we deduce 

Xo^t{x,v) = X + vt + e / v^^\s,x,v) ds + e'^ / v^'^\s,x,v) ds + . . . , 



with = 0) = 0. 
So 

On the other hand, 



eF{t,Xp^t) = e Y^Fit, k) e^i-k:^2i.k-vt^2i.k\e p.vW ds+e^ J^^.i^) ds+..] 
k 

= e V Fit, k) gSi-fc-^e^^"^-''* [1 + 2inek ■ ! v^^^ ds + 2ine'^k ■ ! v^^^ ds 

'- Jo Jo 



(27r)V j\^^Us^\... 



By successive identification. 



Fit, k) e2^-fc-^e2*"^-'^* 
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I. Jo 



dt 



J2 F{t, k) e'^'^^-^'e^'' 



2iiTk-vt 



t \ 2 

ds 



etc. 



In particular notice that 



dt 



<Ek\nt,k)\,so 



/•oo f)„Al) foc ^ poo ^ 

/ ^ dt< ^\F{t,k)\dt< / ^\F{t,k)\e^^^'e-^^^'dt 



27r/i 



So, under our uniform analyticity assumptions we expect Vo,j(x, to be a uniformly 
hounded analytic perturbation of v. 

5.2. Main result. On T^. we consider the dynamical system 

eF{t,xy, 



dt^ 



its phase space is T*^ x M*^. Although this system is reversible, we shall only consider 
t > 0. The parameter e is here only to recall the perturbative nature of the estimate. 

For any {x,v) G T"^ x M"^ and any two times r, t G M+, let Sr^t be the transform 
mapping the state of the system at time r, to the state of the system at time t. In 
more precise terms, Sr,t is described by the equations 

Sr,t{x,v) = {Xr,t{x,v),Vr,t{x,v)); 
Xr^rix, V) = X, Vr^rix, v) = V] 



(5.1) 



^;Xr,t{x,v) = Vr,tix,v), -^K-,t(a;, f) = E F {t, Xr,t{x, v)) . 



dt ' ' dt 

From the definition we have the composition identity 

in particular St^r is the inverse of Sr^t- 

We also write 5"°^ for the same transform in the case of the free dynamics [e = 0); 
in this case there is an explicit expression: 

(5.3) S'{x,v) = {x + v{t-T),v), 
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where x + v{t — t) is evaluated modulo Z*^. Finally, we define the "scattering trans- 
forms associated with eF" : 

(5.4) nt,r = St,r O 5°,. 

(There is no simple semigroup property for the transforms Qt,T-) 
In this section we establish the following estimates: 

Theorem 5.2 (Analytic estimates on scattering transforms in hybrid norms). Let 

e>0 and let F = F{t, x) on M+ x T'^ satisfy 

(5.5) F(t,0)=0, snpi\\F{t, ■)\\jr\t+^ + WVxFlt, ■)\\jr\t+^ < Cf 
for some parameters X, fi > and Cp > 0. Let t > t > 0, and let 

be the scattering transforms associated with eF. Let < X' < X, < fi' < fi and 
t' > be such that 

(5.6) X'{r'-r)<^J^^. 
Let 

/2i(r,t) = CFe-2-(^-^')- min{(t-r); (27r(A - A'))"^} ; 

i?2(r, t) = Cf e"2-(^-^')r _ ^)2/2 ■ (27r(A - A'))"'} . 

Assume that 

(5.7) VO<r<t, eR2{T,t) < ~ ^ \ 
and 

Then 

(5.9) VO<r<t, \\nXtr-ld\\^x'.,' <2eR2iT,t) 
and 

(5.10) VO<r<t, \\nVtr-ld\\^x'.^.' <eRi{T,t). 

Remark 5.3. The proof of Theorem 15. 21 is easily adapted to include Sobolev correc- 
tions. It is important to note that the scattering transforms are smooth, uniformly 
in time, not just in gliding regularity (r' = is admissible in (15.61) ). 
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Proof of Theorem \5.^ For a start, let us make the ansatz 

St,r{x,v) = (^x - v{t - r) + e Zt,T{x,v), v + e drZt^r{x,v)j , 

with 

Ztt{x,v) = 0, drZt^r {X,V) = 0. 

T=t 

Then it is easily checked that 

- Id = £ (Z, drZ) O St,- 

in particular 

||nt,^-id|| v,m' = £ ||(z, 9^z)|La',m' • 

t' t + r'-T 

To estimate this we shall use a fixed point argument based on the equation for 
St,T, namely 

or equivalently 



d'z 



f(j,x - v{t - t) +eZt,-^. 



dr2 

So let us fix t and define 



such that (Z^ is the solution of 

(5.11) <! dr^ 

Zt,t = 0, {drZt,r 



F(T,x-v{t-r)+eWt,r 



= 0. 

T=t 



What we are after is an estimate of the fixed point of \E'. We do this in two steps. 
Step 1. Estimate of ^(0). Let Z° = $(0). By integration of (l5lTil (for W = 0) 
we have 



Z^, = J^{s-t)F{s,x- v{t - s)) ds. 



Let a such that \'a < (// — lJi'')/'2.- We apply the Z^^l^ norm and use Proposi- 
tion SH 

ll^t^rll^^''"' - / (-^ ~ ^) 11-^ (-5' ■ ) IU^'.m' ds= I (s - r) ||F (s, ■ ) II^a'.+av+m' ds. 



-'t + CT 
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Of course AV + /x' < /i, so in particular 

X' s + X'a + 1^' < -(A - A') s + A s + /i. 
Combining this with the assumption F{s, 0) = yields 



So 



(t - r)2 



mm 



(2vr(A-A0)^ 



< R2{r,t). 



With t still fixed, we define the norm 
(5.12) 



1 (^t,^)o<r<t 


:= sup 1 







< r < t; + t > 0; AV < ^- — — 



The above estimates show that ||||^E'(0) |||| < 1. (We can assume t + a > since 
t + (r' — r) > t — r > 0, and we aim at finally choosing a = t' — r.) 

Step 2. Lipschitz constant of We shall prove that under our assumptions, 
^ is 1/2-Lipschitz on^the ball 5(0,2) in the norm |||| ■ ||||. Let W,W e 5(0,2), and 
Z = \i/(Vr), Z = "^iW). By solving the differential inequality ioi Z — Z we get 



1 rt 



Jt 



[s-t) V,F(^s,x-t;(t-s) 

+ e(eWt,s + {i-e)WtA]dsde 



We divide by R2{r,t), take the Z norm, and note that i?2(s,t) < i?2(T, t); we get 



Zt,T — Ztr 



0<T<t 



< e 



0<s<t 



A{t) 
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with 



A{t) = sup 

a,T JQ Jt 



1 rt 



V^F[s,x-v{t-s)+e[eWt,s+{'^-e)Wt,. 
By Proposition 14.251 (composition inequality), 



ds d9. 



A{t) < / (s-r) ||V.F(s, 



^A',/j' + e(t,s,(T) 



ds 



[s-t) \\S/n,F{s, ■)\\jr\u+x'a+t.'+,(t,s,<r) ds, 



with 



e{t,s,a):=e 9Wt,s + {1 - 9)Wt,s ^, < 2 e i?2(s, t) < 2 e i?2(r, t) 



Using (15.71) . we get 

X's + X'a + /i' + e(s, t, a) < X's + X'a + fi' + 2e i?2(r, t) 

< X's + fi={Xs + n)-{X- X')s. 

Using again the bound on V^-f^ and the assumption -F(s, 0) = 0, we deduce 

A{t) < sup l\s - r) Cf e-^'^^^-^') ^ ds < R^iO, t) < yy - 
Using (15. Sp . we conclude that 



0<T<t 



1 

< - 

- 2 



0<s<t 



So is 1/2-Lipschitz on -8(0, 2), and we can conclude the proof of (15. 9p by applying 
Theorem IA.2I and choosing a = t' — t. 

It remains to control the velocity component of Q, i.e., establish (I5.10p : this 
will follow from the control of the position component. Indeed, if we write Qt^r = 
e~^{QVt^r — Id)(a;,f), we have 

Qt,T = j^F[s,x-v{t- s)+e Wt,)j ds 
so we can estimate as before 

llQt.r |La',m' < / \\F{S,-)\\-p\l!, + \l(Tl-T) + ^l + ^(t,8,Tl-T) ds 

t + {r'-r) Jt 
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to get 

Thus the proof is complete. □ 

Remark 5.4. Loss and Bernard independently suggested to compare the estimates 
in the present section with the Nekhoroshev theorem in dynamical systems theory 
[nni EZl- The latter theorem roughly states that for a perturbation of a completely 
integrable system, trajectories remain close to those of the unperturbed system for 
a time growing exponentially in the inverse of the size of the perturbation (unlike 
KAM theory, this result is not global in time; but it is more general in the sense that 
it also applies outside invariant tori). In the present setting the situation is better 
since the perturbation decays. 

6. Bilinear regularity and degay estimates 

To introduce this crucial section, let us reproduce and improve a key computation 
from Section [3l Let G be a function of v, and R a time-dependent function of x 
with R{0) = 0; both G and R may be vector-valued. (Think of G{v) as Vyf{v) and 
of R{t,x) as VW * p{t,x).) Let further 

a{t,x)= [ [ G{v) ■ R{t,x -v{t-T))dvdT. 

Jo JR'i 

Then 

a{t,k)=[ [ [ G{v)-R{T,x-v{t-T))e-^'^''-^dvdxdT 

Jo Jjd. JlRd 

= [ [ [ G[v)-R{T,x)e-^'''''-''e-^'^^-''^'-^UvdxdT 
Jo Jt'^ Jr'' 

= [ G{k{t-T)) ■R{T,k)dT. 

Jo 

Let us assume that G has a "high" gliding analytic regularity A, and estimate a 
in regularity Xt, with A < A. Let a = a{t, r) satisfy 



0<a(t,r) < (A-A) (t-r); 
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then 



< fsUpe2-^l''l|G(r/)A f sup ||i?(r, OI^a.-.) f\-2nl(X-X)(t-r)-a] 

\ V J \0<r<t J Jo 

where we have used 

k^O^ 27r(A(t-r) + < 27rA|A:|(t - r) - ((A - A)(t - r) - a) . 



Let us choose 



/ , (A -A) 
"(^' V = — n mmjl ; t - r j; 



then 



7r(A - A) 
So in the end 

ll^ll A 

7r(A — Aj o<T<t 

where \\G\\^j = sup^(e2-^l''l|G(r/)|). 

In the preceding computation there are three important things to notice, which 
he at the heart of Landau damping: 

• The natural index of analytic regularity of a in x increases linearly in time: 
this is an automatic consequence of the gliding regularity, already observed 
in Section m 

• A bit a{t, r) of analytic regularity of G was transferred from G to R, however 
not more than a fraction of (A — A)(t — r). We call this the regularity 
extortion: if / forces /, it satisfies an equation of the form dtf + v ■ Vxf + 
^[f] ' ^vf = S, then / will give away some (gliding) smoothness to p = 
Ifdv. 

• The combination of higher regularity of G and the assumption -R(O) = 
has been converted into a time decay, so that the time-integral is bounded, 
uniformly as t ^ oo. Thus there is decay by regularity. 
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The main goal of this section is to estabhsh quantitative variants of these effects 
in some general situations when G is not only a function of v and R not only a 
function of t, x. Note that we shall have to work with regularity indices depending 
on t and r! 

Regularity extortion is related to velocity averaging regularity, well-known in ki- 
netic theory [IQ] ; what is unusual though is that we are working in analytic regularity, 
and in large time, while velocity averaging regularity is mainly a short-time effect. 
In fact we shall study two distinct mechanisms for the extortion: the first one will 
be well suited for short times {t — r small), and will be crucial later to get rid of 
small deteriorations in the functional spaces due to composition; the second one will 
be well adapted to large times (t — r ^ oo) and will ensure convergence of the time 
integrals. 

The estimates in this section lead to a serious twist on the popular view on Landau 
damping, according to which the waves gives energy to the particles that it forces; 
instead, the picture here is that the wave gains regularity from the background, and 
regularity is converted into decay. 

For the sake of pedagogy, we shall first establish the basic, simple bilinear estimate, 
and then discuss the two mechanisms once at a time. 

6.1. Basic bilinear estimate. 

Proposition 6.1 (Basic bilinear estimate in gliding regularity). Let G = G{t, x, v), 
R = R{t, X, v), 





Then 





and 



(6.2) 
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Proof of Proposition \6.1[ Obviously (16. 2p follows from (16.11) . To prove (16. ip we apply 
successively Propositions 14. 15^ 14.191 and 14.241 



< 



/ {G-R)oSUdv 
(G ■ i?) o Si 



\G-R\\ < IIGII l|i?II.A,M. 



□ 



6.2. Short-term regularity extortion by time cheating. 



Proposition 6.2 (Short-term regularity extortion). Let G = G{x,v), R = R{x,v), 
and 



(G-R) {x-v{t-T),v) dv. 



Then for any X, fi,t > and any b > —1, we have 
(6.3) \\pyxt+, < \\G\ ^ 

Moreover, if Pk stands for the projection on the kth Fourier mode as in (14.361) . one 
has 



I „A(l+6)./j;l LR _A(1+6),M. 

^ bt bt 



(6.4) 



^27r {\t+^J,)\k 



\m\<J2\\p^G\\ 



_A{1+6),m;1 Pfc-f-R ^A{l+6),M. 

^ bt bt 



' 1+6 



"1+6 



Remark 6.3. If R only depends on t, x, then the norm of R in the right-hand side 
of (16. 3p is ll-Rll^F" with 

bt 



v = X{\^b) 



1 



+ /i = (Ar + /i) - 6(t - r), 



as soon as r > 6t/(l + 6). Thus some regularity has been gained with respect to 
Proposition 16. 11 Even if R is not a function of t, x alone, but rather a function of t, x 
composed with a function depending on all the variables, this gain will be preserved 
through the composition inequality. 

Proof of Proposition \6.2[ The proof presented here relies on commutators involving 
Vt,, Vx and the transport semigroup, all of them classically related to hypoelliptic 
regularity and velocity averaging. Separating the different components of R and G, 
we may assume that both are scalar-valued. 
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Let S = S^_-i., so that R o S{x, v) = R{x — v{t — t),v). By direct computation, 
(6.5) 

tV^{RoS) = (tV^i?)o^ = r((r-6(t-r))V^. + (l + 6)V^)i?j oS-{l + b)V^{RoS). 
Let 

D = Dr,t,b ■■= [r - hit - r)) V,. + (1 + 6)V,. 

Then (16.51) becomes 

(6.6) tV^(i? oS) = (DR) oS-{l+ b)V.,{R o S). 

Since V^; commutes with V„ and D, and with the composition by S as well, we 
deduce from (16. 6p that 

td^^ [(1 + b)''\/tiiD'R) o S)] = (1 + b)''\/';{td^SD'R) o S) 

= (1 + 6)^V^ ((D'+^'i?) o 5) - (1 + 6)'=V^ ((1 + b)d,^ [D'R o S 
= [{l + b)W.]\{D'+^^R)oS) - [{l + b)W,f^^\{D'R)oS). 

So by induction, 

{tV^TiR °S)=Y,r )[-{! + h)V,Y\{D"R) o S). 



(6.7) 



Applying this formula with R replaced by G • i? and integrating in v yields 



(tV^)" / {G-R)o dv = / D"(G ■ R) o 



/ D''{G-R)dv. 



It follows by taking Fourier transform that 



{2nitkYP{k) = / [Z)"(G ■ R)] dv 

= [ ({l + b)\/y + 2m{r-b{t-r))ky {G- R){k,v)dv, 
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whence 

k,n 



k,n 



2iTfi\k 



(A(l + 6))' 



V,, + 2z7r r 



bt 



1 + b 



(G ■ R){k, v) 



L^(dv) 



(1+6), 



6f 
' 1+6 



and the conclusion follows by Proposition 14.241 

Inequality (16 ■4p is obtained in a similar way with the help of Proposition 14.341 □ 



Let us conclude this subsection with some comments on Proposition 16.21 When 
we wish to apply it, what constraints on r) (assumed to be nonnegative to fix 
the ideas) does this presuppose? First, h should be small, so that A(l + 6) < A given. 
But most importantly, we have estimated Gr in a norm Z^i instead of (this is 
the time cheating), where |r' — r| = ht/{l + h). To compensate for this discrepancy, 
we may apply (I4.19p . but for this to work ht/{l + h) should be small, otherwise we 
would lose a large index of analyticity in x, or at best we would inherit an undesirable 
exponentially growing constant. So all we are allowed is h{t, r) = 0(1/(1 + t)). This 
is not enough to get the time-decay which would lead to Landau damping. Indeed, 
iiR = R{x) with ^(0) = 0, then 



\R\\-f\T+iJ.-\b{t-T) < e ^^^^ '^^ \\R\\jr\T+ti] 
so we gain a coefficient e~'^^(*~'^\ but then 



\R\\ _A(l+6),fi 

^T~bt/(l + b) 



Jo 



e-^<^) dr 



1 



which of course diverges in large time. 

To summarize: Proposition 16. 2l is helpful when (t — r) = 0(1), or when some extra 
time-decay is available. This will already be very useful; but for long-time estimates 
we need another, complementary mechanism. 

6.3. Long-term regularity extortion. To search for the extra decay, let us refine 
the computation of the beginning of this section. Assume that Or = V^Qt, where 
(fi'r)r>o solves a transport-like equation, so G(r, k,r]) = 2inrig{T, k,ri), and 



|G(r,A;,r/)|<27r|r/| 



-27r/i|fe| g— 27rA|r;+feT 
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Up to slightly increasing A and fi, we may assume 

(6.8) \G{t, k,ri)\<{l + r) e"'"^!^'! e-'^^l^+'^^l . 

Let then p(r, x) = f /(r, x, v) dv, where also / solves a transport equation, but has 

a lower analytic regularity; and R = VW * p. Assuming |VVr(/i;)| = 0{\k\~^) for 
some 7 > 0, we have 

(6-9) \Rir,k)\< ' 
Let again 

a{t,x)= / / G{t,x — v{t — t),v) ■ R{t,x — v{t — t)) dv dr. 
Jo Jr-^ 

As t ^ +00, G in the integrand of a oscillates wildly in phase space, so it is not 
clear that it will help at all. But let us compute: 

a{t,k)= [ [ [ G{T,x-v{t-T),v) ■ R{T,x-v{t-T))e~'^'^^-''dxdvdT 

Jo Ju'i Jj'i 
t 



G(r, X, v) ■ R{t, x) e"2*"'= " g-'^'^-''^*-") dx dv dr 

Jw^ Jf^ 

= [ [ Cfllir,k,v)e-^'^''''^'-^UvdT 
Jo Jr'^ 

= 11 y2G{T,i,v)-R{T,k-i)e-^'^''''^'~^UvdT 

Jo Jr'' g 

= f VG(r,£,A;(t-r)) ■R{T,k-i)dT. 
Jo ^ 

At this level, the difference with respect to the beginning of this section lies in the 
fact that there is a summation over i e Z*^, instead of just choosing i = 0. Note 
that a{t, 0) = /q // G{t, x, v) ■ R{t, x) dx dv dr = 0, because G is a w-gradient. 
From fl6.8p and (16.91) we deduce 

^e2-(^*+'')l^^l|a(t,A:)| 

k 

/t _ -2TT\T\k-l\ 
[l+r) ^ e^'"^!''! ^2-K\t\k\ ^-27T-p\e\ ^-27TX\k{t-T)+eT\ ^-27T,i\k-e\ J: 



\k-i\^ 
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Using the inequalities 

^-2■K^l\k-l\ g27rAt|fe| g-27rp|£| ^ ^-2■K[Jl-^l)\i\ 

and 



E -2tt{JI~^)\1\ rt 



g-27rAr|fc-^| g27rAi|fe| g-27rA|A:(t-T)+^T| ^ g-27r(A-A)|fc(t-T)+^r| 

we end up with 

-2T7{ji~n)\i\ ft 

|(T(t)||^At+M < ^ ' " 

If it were not for the negative exponential, the time-integral would be O(t^) as 
t oo. The exponential helps only a bit: its argument vanishes e.g. for d = 1, 
k > 0, i < and r = {k/ {k + Thus we have the essentially optimal bounds 

(6.10) f e-2-(A-A)|A:(t-r)+£r| < ^ L 

^ ' Jo - 7r(A - A) 



and 



7r(A-A)|A;-£| 

1 / 1 \ \k\t 



\k-i\ 



(6.11) /" e-2-(A-A)|fc(t-r)+£r| ^ < ^ ^ _1 ^ 

Jo - 27r2(A- A)2|A;-£|2 Vvr(A - A)/ 

From this computation we conclude that: 

• The higher regularity of G has allowed to reduce the time-integral thanks to a 
factor e-"l'=(*-^)+^^l j b^t this factor is not small when r/t is equal to k/{k — i). As 
discussed in the next section, this reflects an important physical phenomenon called 
(plasma) echo, which can be assimilated to a resonance. 

• If we had (in "gliding" norm) ||Gt|| = 0(1) this would ensure a uniform bound 
on the integral, as soon as 7 > 0, thanks to fl6.10p and 

1^ (l + |fc- £1)1+7 <+^- 

• But Gr is a velocity-gradient, so — unless of course G depends only on v — 

diverges like 0(r) as r ^ 00, which implies a divergence of our bounds in large 
time, as can be seen from (16.111) . If 7 < 1 this comes with a divergence in the k 

variable, since in this case Xlfc^ (i+|fc_f|)i+7 = +00. (The Coulomb case corresponds 
to 7 = 1, so in this respect it has a borderline divergence.) 

The following estimate adapts this computation to the formalism of hybrid norms, 
and at the same time allows a time-cheating similar to the one in Proposition 
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Fortunately, we shall only need to treat the case when R = R{t, x); the more general 
case with R = R{t, x, v) would be much more tricky. 

Theorem 6.4 (Long-term regularity extortion). Let G = G{t,x,v), R = R{t,x), 
and 

a{t,x)= / / G{t,x — v{t — t),v) ■ R{t,x — v{t — t)) dv dr. 

Jo jRd 

Let A, A,/i,7I, /i' = n'it.T), M >l such that (1 + M)A >A>A>0, 71>/i'>/i>0, 
7 > an(i 6 = h{t, r) > 0. Then 



(6.12) 
where 
(6.13) 



\(y{t,-)\\pxt+^ < / K^{t,r)\\Rr\\T-dT + / K^{t,T)\\Rr\\ dr, 
Jo Jo 



p = max {\t + jj! — ^h{t — t) ] 



(6.14) 



K^{t,T) = e 



-2nl 



G{t, X, ■) dx 



CA(1 + 6);1 



(6.15) 
(6.16) 

Klit.T) 



\G. 



Kf{t, r) = sup 

0<T<t 



T II ^\(l + b),-p 
^T-bt/(l + b) 



l + T 



Kl{t,T) 



[l+r) sup 

k=iO, li^O 



2'^(^)l^l e-2'^(w)lMt-r)+^r| ^_2^ [(M'-/.) + ^(i-r)] |fc-^| ' 

1 + \k-iY< 



Remark 6.5. It is essential in (16.121) to separate the contribution of Gir^ 0, v) from 
the rest. Indeed, if we removed the restriction £ 7^ in (16.161) the kernel K\ would 
be too large to be correctly controlled in large time. What makes this separation 
reasonable is that, although in cases of application Gir^ x, v) is expected to grow like 
0(r) in large time, the spatial average / G{r,x,v) dx is expected to be bounded. 
Also, we will not need to take advantage of the parameter 7 to handle this term. 



Proof of Theorem \6.4 Without loss of generality we may assume that G and R are 
scalar- valued. (E.g. choose HGH = sup;^<j<^ llCil, \\R\\ = Ei II ^il, where G = 
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{G\ . . . , R = {R\ ...,R'^); then it suffices to bound G' R' for all i.) First 
assume G{t, 0, v) = 0, and we write as before 



(6.17) \a{t,k)\< I I G{Tj,v)e-^'^^-''^'-^Uv 

Jo ^_^Q J 



\R{T,k-i)\dT. 



Next we let t' = t — b{t — r) and write 



(6.18) e^''^^*^''^!''! < e"^''^^"^^l^l g-27rA(r-r')|fc-£| g-27r(M'-A.)|fe-^| g-27r(A-A) |fc(t-r')+^T'| 

g27r77|£| g27r(Ar+At')|fc-^| g27rA|fe(t-T')+^T'| 



Since < A — A < MA, we have 



-27r{X-\)\k{t-r')+er'\ < g-2^(A^) |fc(t-r)+^r| 2^A(^_^')|fe_/| . 



SO (16.181) implies 



(6.19) e2^(^*+^)l*^l < e-2'r(P-^)'^l e"^^(^)l''(*"^)+^^l g-2^[(M'-M) + t(r-r')] 



g2.pKI ^2niXr+^')\k-e\ ^ | (2mA) (A:(t-rO+fr^)|" 
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For each n G Nq, 
(2i7rA) {k{t - r') + ir') 



G(r, i, v) e-^^-'^-C*-) dv 



r 



G(r, £, v) [2t7r{k{t - r') + £t')] " e-^^^'^'^^*"^) dv 

t - T 



< 



t - 


t' 


t- 




t - 


t' 


t - 




t - 


t' 


t - 




t - 


t' 



G(r,£,w) {2m)[k{t - t) + £r' 

t - T 



t-T' 



-2iTTk-v{t-T) 



n\ \t — T 



Vv + 2i7rir 
V„ + 2z7r£ ( r 



V„ + 2i7iiT' 

,[t-T 



-2i-Kk-v{t-T) 



t-T' 

, f t-T 
t-T' 

bt 



t - t' 

G{t, i, v) e-2^-'=-(*--) dv 



1 + b 



G{T,i,v] 



G{T,i,v) 



L^{dv) 



Combining this with (16.171) and fl6.19p . summing over /c, we deduce 



< 



kin 



A"(1 + fe)" 



\R{T,k 



i + \k-e\^ 

V„ + 2i7r£ ( r 



,2n-p\e\ 2n[[XT+f,')-^b[t-T)] \k-l\ 



bt 



1 + b 



GiT,e,v] 



dT, 



L^{dv) 



and the desired estimate follows readily. 

Finally we consider the contribution of G{t,0,v) = J G{t, x,v) dx. This is done 
in the same way, noting that 

g-27r(A-A)|fc|{t-T) 

sup 



k^O 



l + \k\' 



^ -27r(A-A)(t-r)_ 



□ 
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To conclude this section we provide a "mode by mode" variant of Theorem I6.4t 
this will be useful for very singular interactions (7 = 1 in Theorem 12. 6p . 

Theorem 6.6. Under the same assumptions as Theorem 6^, for all k ^l/" we have 
the estimate 

(6.20) e2'^(^*+'^)l*^' A;) | < f K^{t, r) (e^^^l^'l \R{t, k)\) dr 

Jo 

where Kq is defined by fl6.14p . z/ by (16.131) . and 



\\G\\ A(i+b),7r;i 



Kk/t^ r) = sup 

0<T<t 



-6t/(l+6) 



1 + r 



^fc/(i, t) 



l + \k-l\^ 



Proof of Theorem \6.(A The proof is similar to the proof of Theorem 16.41 except that 
k is fixed and we use, for each the crude bound 



,27rjl|£| 



V„ + 2iTli. T 



bt 



1 + b 

2n-p\j 



G{T,i,v] 



L^{dv) 



+ 2i7rj T 



bt 



1 + b 



GiT,j,v] 



□ 



7. Control of the time-response 



To motivate this section, let us start from the linearized equation (13.31) . but now 
assume that f^ depends on t, x, f and that there is an extra source term S, decaying 
in time. Thus the equation is 



dt 



V ■ VJ - {VW * p) ■ V,J° = S, 
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and the equation for the density p, as in the proof of Theorem 13.11 is 
(7.1) 

p{t,x) = j fi{x—vt,v)dv+ 



[ [ VJ^{T,x-v{t-T),v)-{VW*p){T,x-v{t-T))dvdT 
Jo Jr'' 

S{t, X — v(t — t),v) dvdr. 



Hopefully we may apply Theorem 16.41 to deduce from (17.11) an integral inequality 
on ip{t) := ||p(t)||jpAt+M, which will look like 

(7.2) if{t)<A + c [ K{t,T)(p{T)dT, 

Jo 

where A is the contribution of the initial datum and the source term, and K(t, r) a 
kernel looking like, say, (16.161) . 

From (17. 2p how do we proceed? Assume for a start that a smallness condition of 
the form (a) in Proposition 12.11 is satisfied. Then the simple and natural way, as in 
Section [HI would be to write 

ip{t) <A + c([ K{t,r) dr] ( sup ^pir] 



0<T<t 



and deduce 

A 



(7.3) ifit) < 



c Jq K{t, t) dr 



(assuming of course the denominator to be positive). However, if K is given by 
(16.161) . it is easily seen that K(t,T) dr > nt as t oo, where k > 0; then 
(17. 3p is useless. In fact (17. 2p does not prevent ip from going to +oo as t — > cxd. 
Nevertheless, its growth may be controlled under certain assumptions, as we shall 
see in this section. Before embarking on cumbersome calculations, we shall start 
with a qualitative discussion. 



7.1. Qualitative discussion. The kernel K in (16.160 depends on the choice of 
fi' = fi(t, t). How large /i' — /i can be depends in turn on the amount of regularization 
offered by the convolution with the interaction VW . We shall distinguish several 
cases according to the regularity of the interaction. 
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7.1.1. Analytic interaction. If 'VW is analytic, there is cr > such that 

Vi^ > 0, Hp * vw^||^.-+<T < c iipIIjti.; 

then in ( 16. 16^ we can afford to choose, say, fi' — fi = a, and 7 = 0. Thus, assuming 
b = B/{l + t) with B small enough so that (/i' — jj) — \b(t — r) > a/2, K is bounded 
by 

(7.4) Z^"^ (t, r) = (1 + r) sup e""!'! e""!^-'! ^-a\kit-r)+ir\ ^ 

where a = ^minjA — A ; /I — /i ; a}. To fix ideas, let us work in dimension d = 1. 
The goal is to estimate solutions of 

(7.5) ip{t) <a + c K^'"^ (t, r) (^(r) dr. 

Jo 

Whenever r/t is a rational number distinct from or 1, there are k,i & Z such 

that \k{t — t) + irl = 0, and the size of ° (t, r) mainly depends on the minimum 
admissible values of k and k — i. Looking at values of r/t of the form l/(?7, + 1) or 
n/{n + 1) suggests the approximation 

(7.6) Z^"^ < (1 + r) min {e-"(^) e'^" ; e'^^C^) e""}. 

But this estimate is terrible: the time-integral of the right-hand side is much larger 

than the integral of i^' .In fact, the fast variation and "wiggling" behavior of K 
are essential to get decent estimates. 




Figure 5. the kernel K (t, r), together with the approximate upper 
bound in ([71]), for a = 0.5 and t = 10, t = 100, t = 1000. 
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To get a better feeling for K , let us only retain the term in/c = l,£=— 1; this 
seems reasonable since we have an exponential decay as /c or £ go to infinity (anyway, 
throwing away all other terms can only improve the estimates). So we look at 
^H(t,r) = (1 + r) e-3°e-°l*-2r|. Let us time-rescale by setting kt{e) =tK^''\t,te) 
for 9 G [0, 1] (the t factor because dr = td6); then it is not hard to see that 

t ' 2a ^ 

This suggests the following baby model for 07.51] : 

(7.7) ip{t)<a + ct^(^^ 

The important point in (17.71) is that, although the kernel has total mass 0{t), this 
mass is located far from the endpoint t = t; this is what prevents the fast growth 
of ip. Compare with the inequality (p(t) < a + ctip{t), which implies no restriction 
at all on ip. 

To be slightly more quantitative, let us look for a power series $(t) = "^i^dkt^ 
achieving equality in (17.71) . This yields = a, a^+i = ca^ 2~^, so 



(7.8) *W = aE 



2fc(fc-i)/2 ■ 

fe=0 

The function $ exhibits a truly remarkable behavior: it grows faster than any poly- 
nomial, but slower than any fractional exponential exp(ct''), v G (0,1); essentially 
it behaves like (as can also be seen directly from (17. 7p ). One may conjecture 

that solutions of (17.51) exhibit a similar kind of growth. 

Let us interpret these calculations. Typically, the kernel K controls the time 
variation of (say) the spatial density p which is due to binary interaction of waves. 
When two waves of distinct frequencies interact, the effect over a long time period 
is most of the time very small; this is a consequence of the oscillatory nature of 
the evolution, and the resulting time-averaging. But at certain particular times, the 
interaction becomes strong: this is known in plasma physics as the plasma echo, 
and can be thought of as a kind of resonance. Spectacular experiments by Malmberg 
and collaborators are based on this effect [311 [55] . Namely, if one starts a wave at 
frequency I at time 0, and forces it at time r by a wave of frequency k — i, a strong 
response is obtained at time t and frequency k such that 

(7.9) k{t -T)+iT = 
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(which of course is possible only if k and i are parallel to each other, with opposite 
directions). 

In the present nonlinear setting, whatever variation the density function is subject 
to, will result in echoes at later times. Even if each echo in itself will eventually 
decay, the problem is whether the accumulation of echoes will trigger an uncontrolled 
growth (unstability). As long as the expected growth is eaten by the time-decay 
coming from the linear theory, nonlinear Landau damping is expected. In the present 
case, the growth of 07.81] is very slow in regard of the exponential time-decay due to 
the analytic regularity. 



7.1.2. Sobolev interaction. If VW only has Sobolev regularity, we cannot afford in 
(16.161) to take fi'(t, r) larger than fi + rjit — T)/t (because the amount of regularity 
transferred in the bilinear estimates is only 0{{t — T)/t), recall the discussion at the 
end of Subsection 16. 2p . On the other hand, we have 7 > such that 

Vz/ > 0, \\VW * pW^..-, <C\\p\\:f., 

and then we can choose this 7 in (16.161) . So, assuming b = B/{1 + 1) with B small 
enough so that (/i' — jj) — \b(t — t)> rjit — r)/ (2t), K in (I6.16P will be controlled by 



(7.10) i^(")''^(t,r) = (1 + r) sup 



where a = ^ min{A — A ; fi — fi; rj}. The equation we are considering now is 

(7.11) (p{t)<a+ [ K^''^'"'{t,T)y^{T)dT. 

Jo 

For, say, r < t/2, we have K^"'^ < K^"^'^\ and the discussion is similar to that in 
17.1.11 But when r aproaches t, the term exp(— a(^) — £|) hardly helps. Keeping 
only A; > and i = —1 (because of the exponential decay in £) leads to consider the 
kernel 

„-a\kt-{k+l)T\ 

i^(")(t,r) = (l + r) sup ■ 

Once again we perform a time-rescaling, setting kt{0) = t K^°'\t,t9), and let t 
00. In this limit each exponential exp(— a|A;t — [k + 1)t|) becomes localized in a 
neighborhood of size 0{l/kt) around 9 = k/{k + 1), and contributes a Dirac mass 
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aXe = k/{k + 1), with amplitude 2/{a{k + 1)); 

t t-.oc' a 1 + (A; + 1)7 (k + 1)2 l"feTT' 

This leads us to the following baby model for ( 17. 11^ : 
(7.12) ^W<« + ^^E^^((l-l)^)- 

If we search for Yl (^nt"' achieving equality, this yields 

ao = a, a„+i = cl (l-^ 

\k>l ^ 

To estimate the behavior of the Ylk above, we compare it with 

^ ^ ~ lo " " ^ ^ function) 

_r(7)r(n + l) 



r(n + 7 + l) \n'r ^ 

All in all, we may expect ip in (17.111) to behave qualitatively like 



$(t) = a 5^ 



Notice that $ is subexponential for 7 > 1 (it grows essentially like the fractional 
exponential exp(t^/''')) and exponential for 7=1. In particular, as soon as 7 > 1 we 
expect nonlinear Landau damping again. 

7.1.3. Coulomb/Newton interaction (7 = 1 When 7 = 1, as is the case for Coulomb 
or Newton interaction, the previous analysis becomes borderline since we expect 
(I7.12p to be compatible with an exponential growth, and the linear decay is also ex- 
ponential. To handle this more singular case, we shall work mode by mode, rather 
than on just one norm. Starting again from (17. ip . we consider, for each k G Z'', 

y,,(t) = e2-(^*+^)l'=l \p{t,k)\, 

and hope to get, via Theorem 16. 6[ an inequality which will roughly take the form 

(7.13) ^k{t)<Ak + c Kk/jt^ 't) ^k-i{r) dr. 

Jo 
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(Note: summing in k would yield an inequality worse than f l7.1ip .) To fix the ideas, 
let us work in dimension d = 1, and set k > 1, i = —1. Reasoning as in subsection 
fl7.1.2p . we obtain the baby model 

ct f kt 

(7.14) ^,(i)<^, + __^,,,(^_ 

The gain with respect to fl7.12p is clear: for different values of /c, the "dominant 
times" are distinct. From the physical point of view, we are discovering that, in some 
sense, echoes occurring at distinct frequencies are asymptotically well separated. 
Let us search again for power series solutions: we set 



{t) — ak,m t^ 



Ofc,0 — ^k- 



m>0 



By identification, a^,™ = a'k+i,m-i c{k + 1) [k/ (k + l))'" ^, and by induction 



A 



k+m C 



k\ 



1+7 



Xk + m)\\ {k + l){k + 2) . . .{k + m) 
We may expect Ak+m ^ then 



A 



k+m 



k\ 



{k + m)\ 



7+2 



^m-1 



a-k,: 



< 



A{ke 



-ak\ im 171 



(m!)T+2' 



and in particular 



Mt)<Ae-^'/'J2 



{ckt)"' 

1^7+2 



a 



ml 



7 



This behaves like a fractional exponential even for 7 = 1, and we can now believe in 
nonlinear Landau damping for such interactions! (The argument above works even 
for more singular interactions; but in the proof later the condition 7 > 1 will be 
required for other reasons, see pp. 11331 and 11441 ) 

7.2. Exponential moments of the kernel. Now we start to estimate the kernel 
}<[{^)n from fl7.10p . without any approximation this time. Eventually, instead of 
proving that the growth is at most fractional exponential, we shall compare it with 
a slow exponential e^*. For this, the first step consists in estimating exponential 
moments of the kernel e~^* j K{t, r) e^"^ dr. (To get more precise estimates, one can 
study e""^*" J K{t, r) e'^'^" dr, but such a refinement is not needed for the proof of 
Theorem 12.61 ) 

The first step consists in estimating exponential moments. 
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Proposition 7.1 (Exponential moments of the kernel). Let 7 G [l,oo) be given. 
For any a G (0, 1), let K^°'^''^ be defined by (17.101) . Then for any 7 < 00 there is 
a = 0(7) > such that if a <a and e G (0, 1), then for any t > 0, 



-et 



[ K^''^'"'{t,T)e'^dT 
Jo 



< C 



Inl 

+ ^ + 



1 



where C = C {'-/). In particular, 



^2^1+7^1+7 



+ 



1 In^ 

-7 + ° 



6 4 + 



a t 

e 2 



a- 



If-f>l ande<a, then e""* [ K^''^'^(t,T) e'^ dr < 

Jo a-^e^ 

7/7 = 1 then e""* f K^"^'"'(t, r) e"" c/r < ( - + ^ ) . 

Jo \e eHj 



-1 ' 



Remark 7.2. Much stronger estimates can be obtained if the interaction is analytic; 
that is, when A'(")''>' is replaced by K^"'' defined in (17.41) . A notable point about 
Proposition 17. II is that for 7 = 1 we do not have any time-decay as t — > 00. 



Proof of Proposition \7.1\ To simplify notation we shall not recall the dependence of 
K on 7. We first assume 7 < 00, and consider r < t/2, which is the favorable case. 
We write 

K^''\t,T) < (1 + r) sup sup e-°l^le-tl'=-^le-°l^(*-")+^"l. 

Since we got rid of the condition i 0, the right-hand side is now a nonincreasing 
function of d. (To see this, pick up a nonzero component of k, and recall our norm 
conventions from Appendix lA.li ) So we assume d = 1. By symmetry we may also 
assume A; > 0. 



Explicit computations yield 



rt/2 

/ e-'I'^^'-'^+^^l (l+r)dr < < 
Jo 



a{i-k) a^£-k)^ 
t 



-akt 



2 



if 



a\k-(.\ 
2 



2kt 



ii(i>k 
l = k 

ii-k<£<k 
1 



a\k — i\ a\k — i]"^ a'^\k 



iil< -k. 
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In all cases, 

ft/2 



-a\k(t-T)+eT\ 



(1 + T)dT < 



3 1 2t 



a\k-e\ a^\k-i\^ a\k 



+ e 



-akt 



t e 



So 



-et 



t/2 



-a\k{t-T)+lT\ 



(1 + r) e"" dr 



_£t 

< e 2 



< e 4 



a\k 



3 1 2t 



-akt 



t e\ 

2 + ~ I l^=fc 



a2|A;-£|2 

where z = sup(xe 

ft/2 



lfc5^£ + e" 



ta Z O Z 



a a" 



-x\ — Q i_ xhen 



Jo 



Using the bounds (for a ~ 

1 



\k~£\ 



3 1 

+ 



a\k-i\ a'^\k~i\^ ae\k-i\ 



ta V \ 



-a\i\ 



Z »Z 
- + — 

a a" 



a 



£2 



0(1), 



we end up, for a < 1/4, with a bound like 



_ , 1 , 1 1 1,1 
Ce"- (— ln- + — + — In- 



^ - , 1 1 



< C 



_ ei 

e -1 



1 1,1 
— + — In- 



e 2 



(Note that the last term is 0{t ^), so it is anyway negligible in front of the other 
terms if 7 < 4; in this case the restriction e < a can be dispended with.) 
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• Next we turn to the more delicate contribution of r > t/2. For this case we 
write 

^-a\k{t~T)+tT\ 



(7.15) /s:(")(t,r) < (1 + r) sup e' 



■a\i\ 



sup 



1 + |A;-£h ' 



and the upper bound is a nonincreasing function of so we assume d = 1. Without 
loss of generality we restrict the supremum to £ > 0. 

The function x i — > (1 + |x — £|''')~-^ ^-a\x{t-T)+tT\ jg (decreasing for x > i, increasing 
for X < —ir/it — r); and on the interval its logarithmic derivative goes 

from 



1. 

tt 



So if t > 7/a there is a unique maximum at z = —ir/ {t — t), and the supremum in 
(17.15P is achieved for k equal to either the lower integer part, or the upper integer 
part of —iT/{t — r). Thus a given integer k occurs in the supremum only for some 
times r satisfying k — 1 < —iT/{t — t) < k + 1. Since only negative values of k 
occur, let us change the sign so that k is nonnegative. The equation 

k-1 < < A; + 1 

t — r 



is equivalent to 



k-1 \ f k+1 , 

t <T < { t. 



k + i-lj \k + i+l 

Moreover, t > t/2 implies k > i. Thus, for t > 7/a we have 
(7.16) 



*/2 e>i k>£ 

For t <'-f/a we have the trivial bound 

-et 



e 



f K(-\t,T)e^^dT<^; 
Jt/2 2a 

so in the sequel we shall just focus on the estimate of fl7.16p . 

To evaluate the integral in the right-hand side of (17.161) . we separate according to 
whether r is smaller or larger than kt/{k + i); we use trivial bounds for e^"^ inside 
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the integral, and in the end we get the exphcit bounds 



-£t 



\ k+l J 
fc-1 ' 

k+e-1 , 



kt 



-et 



' fc+i 



1 kt 



_a{k + i) a{k + iy a^{k + i)\ ' 



All in all, there is a numeric constant C such that (17.161) is bounded above by 



(7.7) cj:e-'j:{-^ 

i>l k>l ^ ^ 



kt 



g 1 1 



together with an additional similar term where e~^^^^^''~^^^ is replaced by e-^^*/(*:+^+i)^ 
and which will satisfy similar estimates. 

We consider separately the three contributions in the right-hand side of (17.171) . 
The first one is 



a2 



sit 

e k+i 



£>1 k>i 

To evaluate the behavior of this sum, we compare it to the two-dimensional integral 



I{t) 



1 



{x + y) 



2+7 



dydx. 



We change variables {x,y) — {x,u), where u{x,y) = ext/{x + y). This has Jacobian 
determinant {dx dy) / {dx du) = {ext)/v?, and we find 



m 



1 



a2 £1+7^1+7 x^+^ 



dx 



e-" u^du = 



a2£l+7tl+7 



The same computation for the second integral in the right-hand side of (I7.17P 
yields 



dx 



et/2 



e-''u^-Uu = 



In- 



ae-y t-y 



(The logarithmic factor arises only for 7 = 1.) 

The third exponential in the right-hand side of (17.171) is the worse. It yields a 
contribution 



(7.18) 



t ^ e k 

^(k + iy+'r' 
e>i k>i ^ ' 
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We compare this with the integral 

t r r e-—yy 



£ Xt 



dx dy, 



and the same change of variables as before equates this with 

dx I e"" n^^"^ du / dx / e~" u' du 



= O 



Ini 



( Again the logarithmic factor arises only for 7 = 1.) 

The proof of Proposition 17.11 follows by collecting all these bounds and keeping 
only the worse one. □ 

Remark 7.3. It is not easy to catch (say numerically) the behavior of fl7.18p . be- 
cause it comes as a superposition of exponentially decaying modes; any truncation 
in k would lead to a radically different time-asymptotics. 




5000 10000 15000 20000 25000 20 40 60 80 



Figure 6. The function (17.181) truncated at £ = 1 and k < K, for 
K = 5, 10, 100, 1000. The decay is slower and slower, but still ex- 
ponential (picture on the left); however, the maximum value occurs 
on a much slower time scale and slowly increases with the truncation 
parameter (picture on the right, which is a zoom on shorter times). 



From Proposition 17. II we deduce exponential bounds: 
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Corollary 7.4 (L^ exponential moments of the kernel). With the same notation as 
in Proposition \7.1\ 



(7.19) r r)' e^'^ dr < 

Jo 



a 



4^1+27 ^2(7-1) 



Proof of Corollary \7.4\ This follows easily from Proposition 17.11 and the obvious 
bound 

r)2 < C (1 + t) J^('")'2^(t, r). 

□ 

7.3. Dual exponential moments. 

Proposition 7.5. VFzi/i the same notation as in Proposition \7.1\ for any 7 > 1 we 

have 

/ 1 ln-\ 

(7.20) .up I e- X<").^«, r) rf* < C(-,) + ^) • 

Remark 7.6. The corresponding computation for the baby model considered in 
Subsection 17.1.21 is 



^ fc>i ^ ^ ^ 



So we expect the dependence upon e in fl7.20l) to be sharp for 7 — 1. 
Proof of Proposition \775\ We first reduce to d = 1, and split the integral as 



poo poo p2T 

/ e-"*K(")'^(t,r)rft = e"W e""* ^("^'^(t, r) rft + e""* ^("^'^(t, r) rft 

Jt J2T Jt 



=:h + h. 

The first term Ii is easy: for 2t < t < +00 we have 

i^'(")-^(t,r)<(l + r) 5^ e-a|^|-f|fc-£|<^^l + ^^ 

fc>l, ^7^0 
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and thus 



POO 



e 



We treat the second term I2 as in the proof of Proposition I7.lt 

A'(")'^(t,r) e""* dt 



p~a\k{t~r)~ir\ 



t>l k>l " \~k+i 



£>1 /c>£ ^ ^ 



We compare this with 

2(1 + t) r r e-'y^ 



a 



POO POO -e-r 



ae^ \ J Ji JO 
^ C \n{l/a) 

where we used the change of variables u = exr/y. The desired conclusion follows. 
Note that as before the term ln(l/a) only occurs when 7 = 1, and that for 7 > 1, 
one could improve the estimate above into a time decay of the form 0(r~^'''~^)). □ 

7.4. Growth control. To state the main result of this section we shall write = 
Z*^ \ {0}; and if a sequence of functions t) {k G Zf, t G M) is given, then 
\\^{t)\\x = EfcC^''^'''' We shall use A:(s)$(t) as a shorthand for (A:(A;, s) t))fc6Zd, 

etc. 

Theorem 7.7 (Growth control via integral inequalities). Let p = P{v) and W = 
W{x) satisfy condition (L) from Subsection \2. 2\ with constants Cq, Aq, in particular 
\fiv)\ < Coe-2'^^«l''l. Let further 

Cw = max \ V \W{k)\, sup \k\ \W{k)\ 
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Let A>0,fi>0,Xe (0, A*] with < A* < Aq. Let t))fcg2d^ t>o be a continuous 
function oft>0, valued in C^*, such that 



(7.21) Vt>0, 



<l>(t)- [ K\t-T)^T)dT 
Jo 



Xt+fi 



< A 



Ko{t,T)+Ki{t,T) 



Co 



1 + r 



|$(r)||A.+;. dr, 



where Cq > 0, m > 1 and KQ{t,T), Ki{t,T) are nonnegative kernels. Let (p(t) = 
\mt)\\xt+^.- Then 

(i) Assume 7 > 1 and Ki = cK^°''>'"' for some c > 0, a G (0,a(7)), where K^"''^'"' 
is defined by ( IT.lOp . and a (7) appears in Proposition Then there are positive 
constants C and x, depending only on 7, A*, Aq, Cq, Cw, uniform as 7 1, such 
that if 



(7.22) 

and 



sup / Ko{t,T) dr <x 
t>o Jo 



(r-t \ 1/2 ^00 

/ Ko(t,r)^rfr ) +sup / Ko{t,T)dt 
Jo J r>0 Jr 

then for any e G (0, a), 



< 1, 



(7.24) Vt > 0, ifit) <CA ii±ill e^"» 



where 

(7.25) T = C max 



1 + — ] e^' e 

a e 



CT ^Cc{l+T^) ^et 



£2+7 



1 

7-1 



C \ / Cq 

e 



(a) Assume Ki = X]i<i<Ar /^'^ some ai G (0,«(1)), where a(l) appears 



in Proposition 7. 1 : then there is a numeric constant F > such that whenever 



N 



a 



3 ' 
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one has, with the same notation as in (i), 



(7.26) Vt>0, (^(t) <CA^^^^:^4^e^'^e^'=(^+^')e^*, 



where 



N 



i=l 




Remark 7.8. Let apart the term cq/ (1 + r)"^ which will appear as a technical cor- 
rection, there are three different kernels appearing in Theorem 17.71 the kernel K^, 
which is associated with the linearized Landau damping; the kernel Ki, describing 
nonlinear echoes (due to interaction between differing Fourier modes); and the ker- 
nel Kq, describing the instantaneous response (due to interaction between identical 
Fourier modes). 

We shall first prove Theorem 17.71 assuming 

(7.27) Co = 
and 

(7.28) / sup|if°(A;,t)|e2"^°l'=l*dt< k G (0, 1), 

Jo k 

which is a reinforcement of condition (L). Under these assumptions the proof of 
Theorem 17.71 is much simpler, and its conclusion can be substantially simplified too: 
X depends only on k; condition fl7.23p on Kq can be dropped; and the factor e^'^{l + 
c/(q;£^/^)) in fl7.24p can be omitted. If < (as for gravitational interaction) 
and /° > (as for Maxwellian background), these additional assumptions do not 
constitute a loss of generality, since 07.281) becomes essentially equivalent to (L), 
and for co small enough the term co(l + t)~™ can be incorporated inside Kq. 

Proof of Theorem \7. 7| under (17.231) and (17.281) . We have 

(7.29) ifit) <A + J^{\K'\{t - r) + Ko{t, r) + K,{t, r)) ifir) dr, 

where |A'''(t)| = sup^ \K^{k,t)\. We shall estimate ^9 by a maximum principle 
argument. Let tplt) = i?e^*, where B will be chosen later. If ip satisfies, for some 



96 



C. MOUHOT AND C. VILLANI 



T > 0, 
(7.30) 



'(p{t)<tp{t) forO<t<r, 



^ (\K'\{t,T) + Ko{t,T) + Ki{t,T)^^iT)dT for t>T, 



then u{t) := tp{t) — ip{t) is positive for t <T, and satisfies u{t) > K{t, r) u{t) dr 
for t > T, with K = \K^\ + Kq + Ki > 0; this prevents u from vanishing at later 
times, so u > and < ip. Thus it is sufficient to estabhsh fl7.30p . 

Case (i): By Proposition 17. H and since J{\K^\ + Kq) dr < 1 — k/2 (for x < '^/2), 



(7.31) A+ I [\K^\{t,T)+Ko{t,T)]ijiT)dT + C I K^''^'''{t,T)ij{T)dT 

cC{-f) 



<A + 



1-1 

2 



a 



3 ^1+^ ^^-1 



Be' 



et 



For t>T := (4 cC (a^ ^1+7^)) 1/(7-1)^ this is bounded above by A + (1 - K/A)Be 
which in turn is bounded by B e^* as soon as 5 > 4 A/k. 
On the other hand, from the inequahty 



et 



sup (/9(r) + c (1 + t) / lp{t) dr 

0<T<t 



we deduce 



m< {l + t)e 



In particular, if 



with c' = c'{c, k) large enough, then for < t < T we have ip{t) < il){t)/2, and fl7.30p 
holds. 

Case (ii): Ki = ^ q /T'-"'-''^. We use the same reasoning, replacing the right-hand 
side in f l7.3ip by 



A 



/ N N 



. a^e ^ af t 

1=1 * 1=1 ' 



Be 



Et 
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To conclude the proof, we may first impose a lower bound on e to ensure 



TV 



(7.32) 



i=l 



^ - 8' 



and then choose t large enough to guarantee 

AT 



(7.33) 



this yields (ii) 



< -: 



i=l 



a: 



feH - 8 



□ 



Proof of Theorem 7.7 in the general case. We only treat (i), since the reasoning for 
(ii) is rather similar; and we only establish the conclusion as an a priori estimate, 
skipping the continuity/approximation argument needed to turn it into a rigorous 
estimate. Then the proof is done in three steps. 

Step 1: Crude pointwise bounds. From fl7.2ip we have 



(7.34) y.{t) = Y,mk,t)\e'- 



< A + V f\K\k,t-T)\ e2-{^*+^)l'=l \<^{t,T)\dT 

I. Jo 



+ 



<A + 



Ko{t,T)+K^{t,T) + 



Co 



1 + r 



(/?(r) dr 



sup \K^{k,t — r)| e 

k 



+ Ki{t,r)+Ko{t,T) 



27rA(f-T)|fc| 



Co 



'1 + r) 



(P{t) dr. 



We note that for any k eTJ^ and t > 0, 

|i^°(A;,t - r)| e'"^!^'^*-") < An^ \W{k)\CQe~^^^^''~^'^\^\' \k\'' t 

CCo f I ; I ifT?/; \A / Cq Ciy 

- — ^ 1^1 1^"^^)! - — r- 

■^0 — J Ao — A 
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where (here as below) C stands for a numeric constant which may change from hne 
to hne. Assuming J K^it, t) dr < 1/2, we deduce from (17.341) 



^{t) <A + l(snp vir)) +C [ ( 
^ \0<T<t J Jo \ 



and by Gronwall's lemma 

(7.35) ^(t)<2Ae^('^*+^(*+*')+'=°^" 
where Cm = Ji^{l + r)"™ dr. 



Step 2: bound. This is the step where the smallness assumption (I7.22p will be 
most important. For all k G Ijf, t > 0, we define 

(7.36) ^fc(t) = e""* $(A;, t) e2-(^*+^)l'=l , 

(7.37) /C°(t) = e-^* K°{k, t) e2-(^*+A')l*^l, 

(7.38) Rk{t) = e-=* (^$(fc, t)- K\k, t - r) $(fc, r) d^ e^-^^^+^^l'^l 

= (vl/,-M/,*/C°)(t), 

and we extend all these functions by for negative values of t. Taking Fourier 
transform in the time variable yields -R^ = (1 — /C°) ^k'-, since condition (L) implies 
|1 — Al^^l > K, we deduce ||^'a:||l2 < ll-RfclUa, ^-c, 

(7 qQ^ Ml II ^ \\Rk\\L^{dt) 

Plugging (I7.39P into (I7.38p . we deduce 



(7.40) VA; G Z^, \\^>k - Rk\\^,^^^) < ''^"'J^'''^"'^ \\Rk\Widt). 
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Then 



(7.41) 



-et\ 



< 



< 



El*' 

k 
k 



+ ^ \\Rk — ^k\\L2{dt) 
L^idt) k 

L^(dt) \ t&i j 



(Note: We bounded \Ri\ by || which seems very crude; but the decay of 

K^y. as a function of k will save us.) Next, we note that 



1^2 



k\\L^{dt) 



Jq 



< 471^ \W{k)\ 



Co 



(Ao - A)^ 



so 



Co 



A)^ 



Plugging this in (17.411) and using (17.211) again, we obtain 
(7.42) 

II /,\ -et II ^ / 1 , C Co Cw 

\\ip{t)e 



< 1 



^ ' k(Ao- A)2 



< 1 



CCoC 



'^w 



«:(Ao-A)2 



Co 



'1 + t) 



2 p 

(/9(r) (ir dt 



We separate this (by Minkowski's inequality) into various contributions which we 
estimate separately. First, of course 



(7.43) 



A 



2e 



100 



C. MOUHOT AND C. VILLANI 



Next, for any T > 1, by Step 1 and /* Ki{t, r) dr < Cc{l + t)/ 



a, 



(7.44) 



„-2et 



Ki{t,T)ip{T)dT ] dt 



2 ^l 



sup ip(t) 

0<t<T 



T / rt 

„-2et 



Ki{t,T)dT^ dt 



Invoking Jensen and Fubini, we also have 



(7.45) 



oo / rt 

„-2et 



2 ^1 



Ki{t,T)(p{T)dT] dt 



oo / rt \ 2 2 

Ki(t,r)e-"(*-")e-"Xr)rfr ) dt 

oo / rt 



ft 



< ( sup / e-"*is:i(t,r)e"^rfr 
.t>r Jo 





oo rt 




Ki{t, r) e-"(*-") e-'"Xr)2 rfr rft 

'T Jo J 



t \ 2 / /"OO r + oo N 

^+ ^^ ^sr^^\ I I / ^^^^^ e""(*"") e-2"Xr)2 rft 

J max{r ; T} / 



sup / e-'' Ki{t,T)e'^ dr 

t>T Jo 



< ( sup / e""* (t, r) e"" dr sup / e"" A'l (t, r) e""* rft 

,t>T io / \r>0 Jr 



oo \ I 

e~^'^ ^{rf dr 
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(7.46) 



~2et 



Ko{t,T)(p{T)dT] dt 







< ( sup / e-'^Ko{t,T)e'^dT 
.t>o Jo 



sup 

T>0 Jt 



e'^Koit,T) e""* dt 



^{rf dr 



t>o Jo 



< I sup / A'o(t, r)rfr I sup / Ko{t,T)dt 



T>0 Jt 



e-""Xr)"rfr 



The last term is also split, this time according to r < T or r > T: 



2 



(1 + r) 

1 

< Co — ^ e LV ^0-^ J 'i C„ 
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and 
(7.48) 



= Co 

< Co 

< Co 
= Co 
= Co 



Co (p{t) dr 
(1 + r)™ 



dt 



oo / pt 



(1 + r)'" 



dr I dt 



t g-2e(t-r) 

r (1 + r)2- 
e-^''ip{tfdt 



dr 



LpirYdT dt 



foo g_2e(t-T) 




JT 

oo 



(l + r)2'^ 
1 



(l + r 



dr 



2m. 



-2^(*-")dt dr 



(l + r 



2m 



-2"ds 



'-^2m Co 



^pitfdt 



Gathering estimates (17.431) to (I7.48p . we deduce from (I7.42p 
(7.49) 



-et\ 



\L'2{dt) 



< 1 



"^H^ 



fi:(Ao - A)^ 



CA 



Vae /J 



lL2(dt)' 



where 



«:(Ao-A)2 



t>r Jo 



sup / e"^* JsTi (t, r) e^^ dr sup / e^^ (t , r) e"^* dt 



T>0 



t>o Jo 



+ ( sup / Ko{tiT)dT\ I sup / Ko{t,T)dt\ + 



T>0 



'^2m Co 



Using Propositions 17.11 (case 7 > 1) and 17.51 as well as assumptions (I7.22p and 
( ]7.23p . we see that a < 1/2 for x small enough and T satisfying ( 17.25^ . Then from 
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(^rm follows 

hit) 



-et\ 



< 1 



C Co Cw \ C A 



«:(Ao-A)2 



a e 



+ Co C„ 



Step 3: Refined pointwise bounds. Let us use (17.211) a third time, now for t > T: 
(7.50) 

e-''^{t) <Ae-'' + (^sup |i^°(A;, t - r)| e^"^^*"")!'^!^ V'(r)e~""c/r 



+ 



i^o(t,r) + 



Co 



L 

t 



(1 + r)™ 
/ (e~"* Ki {t, t) e"") (^(r) e""" cir 



sup |KO(A;,t-r)|e'"^(*-")l^l ) dr 

VfcGZ^ 



+ / Ko{t,TfdT 



/ V^O 

e-2^*A'i(t,r)2e2^"dr 





+ 



2m 



^(r)'e-'""dr 



We note that, for any G Z^, 



\K\k,t)\e 



2nX\k\t 



< 167r^ |1?(A;)|2 \f{kt)\^ | A;!^^ g4.A|fc|i 

< C Co' \W{k)\^ e-MAo-A)|/=|t |^|4 ^2 



< 



(Ac - A)^ 



\Wik)\'e 



2 ^-27T(Xo-X)\k\t |^|2 



ri2 

^ ^ ri2 ^-27r(Ao-A)|fc|t 

^ ~ TTTT i-^W fc- 



< 



(Ac - A) 

cci 

(Ao - A)^ 



2 



(^2 -27r(Ao-A)t. 
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SO 

2 

,2^\(t-r)\k\ \ ^ ^ ^O^W 



[ (sup \K\k,t-T)\e'^ 
Jo ykezi 



(Ao - A)3 

Then the conclusion follows from (I7.50p . Corollary 17.41 conditions (I7.25P and (17.231) . 
and Step 2. □ 

Remark 7.9. Theorem 17.71 leads to enormous constants, and it is legitimate to 
ask about their sharpness, say with respect to the dependence in e. We expect the 
constant to be roughly of the order of 

sup (e'^'^*''^''^ e~^*) ~ exp (^e~^^ . 

Our bound is roughly like exp(£:~'^^"''^'>')/*^'^~^)); this is worse, but displays the expected 
behavior as an exponential of an inverse power of e, with a power that diverges like 
0((l-7)~^) as 7^1. 

Remark 7.10. Even in the case of an analytic interaction, a similar argument 
suggests constants that are at best like (In l/e)''^^/"^, and this grows faster than any 
inverse power of 1/e. 

To obtain sharper results, in Section [TT] we shall later "break the norm" and 
work directly on the Fourier modes of, say, the spatial density. In this subsection we 
establish the estimates which will be used later; the reader who does not particularly 
care about the case 7 = 1 in Theorem 12.61 can skip them. 

For any 7 > 1, a > 0, A;, ^ G Z'^ \ {0} = and < r < t, we define 

(7.51) r) = ""^"'^ \ + \k-i/ • 

We start by exponential moment estimates. 

Proposition 7.11. Let 7 G [1, 00) be given. For any a G (0, 1), k,i G Zf, let K^j''^ 
be defined by (17.511) . Then there isa = 0(7) > such that if a <a and e G (0, a/4) 
then for any t > 



(7.52) sup Ve^^* f Kif"'{t,T)e'^dT < 

k&i , Jo 



a 



(7.53) sup Ve"^*f/ 4t^(t,r)2e2-dr') 

k&i ^g^d \Jo ' J 



1 



Q,d ^7+2 p 2 ' 
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(7.54) sup V sup e"" / ^ e""* dt < 



Proof of Proposition \7.1i\ We first reduce to the case d = 1. Monotonicity cannot 
be used now, but we note that 

Kif^it, r)< J2 e"""'' e-"l'^l . . . e-"l^-il Kj^^^it, r) e-^'^^+^l . . . e""'^^!, 

l<i<'i 

where Ku e- stands for a one- dimensional kerneL Thus 



sup 

k 



t / \ t 

E / "(^' dr < sup 5: e-l™l E E / 

- ^^^1 sup sup E f e-* (t, r) c/r. 



In other words, for fl7.52p we may just consider the one-dimensional case, provided 
we allow an extra multiplicative constant C{d)/a'^~^. A similar reasoning holds for 
(17.531) and (17.541) . From now on we focus on the case d = 1. 

Without loss of generality we assume k > 0, and only treat the worse case £ < 0. 
(The other case k,i > is simpler and yields an exponential decay in time of the 
form e^'^ ™"^^°'''^^*). For simplicity we also write K^^e = K^j'^ . An easy computation 
yields 



/ KkAt.r)e'^dT 
Jo 



Ce-"l^l / 1 \k\t 1 \ 

< \ — \ e TFT? 

- l + \k-i\^ \a\k -i\ a\k - £|2 a^\k - 



Then for any k > 1, we have (crudely writing = 0(a)) 
(7.55) E / e~''KkAt.r)e'^dT 



^ I ^ a2(fc + £)i+7 + 2-. a (A; + £)2+7 ) • 
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For the first sum in the right-hand side of f l7.55p we write 



(7.56) 



E 



£> 



1+7 



< 



< 



i>i 



eit 
k + 1 



1+7 



e k+e 



(sty- 



For the second sum in the right-hand side of fl7.55p we separate according to 
1 <i <k OT i>k + l: 



(7.57) 



E 

l<£<k 



—ML e gf 



(A; + £) 



2+7 



kt<Y^ 



i<e<k 



(fc + 1) 



2+7 



C 
< — 

a 



< 



et 
k + l 



1+7 



k 



C 



k+l J {ety+^ 



(7.58) 



£>/c+l 



(k + i) 



2+7 



kt<C J2 



e 2 



€>fc+l 



^2+7 



kt 



^ C e-^'l^ ^ 



C 



The combination of ( 17.55^ ( ]7.56p . ( 17.57^ and ( 17.58^ completes the proof of (E 
Now we turn to ( 17.53^ . The estimates are rather similar, since 

KkAt^ry <C{l + t)Kk,eit,T) 
with 7 ^ 27 and a —>■ 2a. So fl7.55p should be replaced by 



(7.59) $^e-^M / KkAt^rfe'^^dT 



< 



E 

e>i 



—a£ - 

e e 



^(A;t)i/2(l + t)V2 



a 



1/2 {k + iy+^ 



ON LANDAU DAMPING 



107 



For the first sum we use (17.561) with 7 replaced by 7 — 1/2: for t > 1, 

(7.60) {i+ty/'Yl 



(A; + ^)l+(^-V2) - (£t)7+l/2 - £7+1/2^7- 

For the second sum in the right-hand side of (I7.59P we write 



i<e<k 



e '=+1 



et 



k + 1 



7+1/2 



(1 + A;)V2 (5t)7+l/2 



< 



C 



a 



and 



E 



„-T^ ^1/2^ 



(fc + £)i+^ 



With this (17.531) is readily obtained. 

Finally we consider (I7.54p . As in Proposition 17.51 one easily shows that 

C 



POO ^ 

sup sup e^^ / e"^* Kk,i{t, r) dr < — - > 



-a\l\ 



< 



Then one has 



2t 



a^k aek ak 



So the problem amounts to estimate 



sup 



:i + r) 



e e k 



ak{k + i)'r 



< 



_a ei{k + i)^ 



e k 



< c 



e 



and the proof is complete. 

We conclude this section with a mode-by-mode analogue of Theorem 17. 7[ 



□ 



Theorem 7.12. Let p = f^{v) and W = W{x) satisfy condition (L) from Sub- 
section with constants Cq,\q,k; in particular \f^{r])\ < Cq e"^'^''*"''^'' . Further 
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let 



Cw = max I 1^(^)1' sup \k\ \W{k)\ 

Let {Ak)k&i, > 0, A G (0, A*] with < A* < Aq. Let {^{k,t))i^^id^ (>o he a 
continuous function oft>0, valued in C^*, such that for allt>0 and k G Z^, 



(7.61) 



^{k,t)- [ K\k,t-T)^{k,r)dr < Ak+ [ A'o(t, r) e^'^^'+^^l^l |$(A;, r)| dr 
Jo Jo 



^27r(Ar+Aj)|fe-£| 



mk-i,T)\dT, 



where c > 0, q > G Tj'^), m > 1, 7 > 1, Ko{t, r) is a nonnegative kernel, K^"^'"' 
are defined by fl7.5ip . a < 0(7) defined in Proposition \7.11 Then there are positive 
constants C and x, depending only on 7, A*, Aq, c := max{^^ q, c^)"'^''^}, 
Cw, such that if 



(7.62) 

and 



sup / KQ{t, t) dr < X 
t>o Jo 



(ft \ 1/2 poo 

/ Ko{t,TfdT] +sup / Koit,T)dt 
Jo J T>0 Jt 

then for any e G (0, a/4) and for any t > 0, 



(7.64) sup (e2-(^*+'^)l'=l \^{k,t)\) <CA^^^'''^ 



where A := (sup^Afc) and 
(7.65) T = (7 max 



1 + 



< 1, 



3+2d g-7+2 



0,^57+2 



Proof of Theorem 7.12. The proof is quite similar to the proof of Theorem 17.71 so 
we shall only point out the differences. As in the proof of Theorem 17.71 we start 
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by crude pointwise bounds obtained by Gronwall inequality; but this time on the 
quantity 

(^(t) = sup |$(A;,t)|e'"(^*+^)l*^l. 

k 

Since K^^iit, r) = 0((1 + r)/a), we find 

(7.66) ifit) <2A e^(^flF^*+^(*+*')+""^™) . 



Next we define /C°, Rk as in Step 2 of the proof of Theorem 17. 71 and we deduce 



(17:391) and ([730]). Let 
(7.67) 



(t) = e2-(^*+'^)l^l \^{k,t)\, 



then 



yk(t) e 



-et\ 



LHdt) ^ W^kh^dt) ( 1 + 
< ||-Rfe||L2(dt) ( 1 



C Cw C*o 



whence 
(7.68) 

\Wk{t) e 



-£t\ 



iL'^idt) 



< 1 



C Co Cw 



-2et 



Ak+ Ko{t,T)(pk{T)dr 



iPk~i{r) dr dt 



We separate this into various contributions as in the proof of Theorem 17.71 In 
particular, using (\7Mh and jlY.iKk,idT = 0{{l + t)/a^), we find 

(7.69) 



t) ipk-eir) dr j dt 



< 



sup sup if kit) 
k 0<t<T 



f£^i^fc,,(^,r)rfr'j dt 



< CA- 



2 ^3/2 



e 
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Also, 



-2et 



2 1 

2 



Kk/{t, r) (pk~i{T) dr dt 



and the last term inside parentheses is 



^ / / KkAt. r) e-^(*--) rft ^^-lirf dr 

£ Jo \Jmax{T;T} / 

^ (^S^ '^P/ sup j e-^'^Mrfdr 

The computation for Kq is the same as in the proof of Theorem 17. 71 and the terms 
in (1 + r)~"^ are handled in essentially the same way: simple computations yield 



(7.70) 



-2et 



T 



2 .1 

2 



dr dt 



(1 + r)" 
< ( sup supipeir) ] \ e 

\0<r<T I / [Jo 

1^ 



(1 + r)™; 



dt 



and 
(7.71) 



-2et 



* J2iCe(Pk-e{r) dr 
(1 + r)™ 



< <j ^sup i^j^ e-'^^ Mr? dr) [J^ce 



roo pt ^-2e{t-T) 




JT 



(1 + r) 



2m 



(ir dt 



) 



< c 



2rn 



^ ^2m~l 



sup 

e. Jo 



+ 00 



(fe{r?dr 
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All in all, we end up with 
(7.72) 



sup \\(pkit) e 

k 



-et\ 



< 1 + 



'^W 



C A 



e 



+ cC„ 



a sup ||v5fc(t) e 

k 



-et\ 



where 



C Cq Cw 



k(Ao - A)^ 



sup^ I 

, t>T Jo 



^sup f 



t>0 Jo 



+ ( sup / Ko{t,T)dT\ (sup / Ko{t,T)dt 



T>0 



Applying Proposition I7.1H we see that a < 1/2 as soon as T satisfies fl7.65p . and 
then we deduce from f l7.72p a bound on sup;. ||v5fc(t) e~'^%L'^(dt)- 

Finally, we conclude as in Step 3 of the proof of Theorem 17.71 from (17.611) , 

(7.73) 



sup 

Vfcez^ 



|ir°(A;,t-r)|e2"^(*-")l^^l ) dr 

dr 



Ko{t, ry dr 



2 m 



+ e-'^'K,,,{t,Tre'^^dr 
and the conclusion follows by a new application of Proposition 17.111 



sup / (pk[T) e 

k Jo 



□ 



8. Approximation schemes 



Having defined a functional setting (Section H]) and identified several mathemat- 
ical/physical mechanisms (Sections [5] to [7]), we are prepared to fight the Landau 
damping problem. For that we need an approximation scheme solving the nonlinear 
Vlasov equation. The problem is not to prove the existence of solutions (this is much 
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easier), but to devise the scheme in such a way that it leads to relevant estimates 
for our study. 

The first idea which may come to mind is a classical Picard scheme for quasilinear 
equations: 

(8.1) dtr^' + V ■ v.r+i + Fin ■ v.r+^ = o. 

This has two drawbacks: first, f'^'^'^ evolves by the characteristics created by Flf^], 
and this will deteriorate the estimates in analytic regularity. Secondly, there is no 
hope to get a closed (or approximately closed) equation on the density associated 
with f''+^. More promising, and more in the spirit of the linearized approach, would 
be a scheme like 

(8.2) dtr^' + V ■ v,.r+^ + ■ v,r = o. 

(Physically, /""^^ forces /"", and the question is whether the reaction will exhaust 
jn+i large time.) But when we write (18. 2p we are implicitly treating a higher 
order term (Vyf) of the equation in a perturbative way; so this has no reason to 
converge. 

To circumvent these difficulties, we shall use a Newton iteration: not only will 
this provide more flexibility in the regularity indices, but at the same time it will 
yield an extremely fast rate of convergence (something like 0(£:^")) which will be 
most welcome to absorb the large constants coming from Theorem 1 7. 71 or Theorem 

[713 

8.1. The natural Newton scheme. Let us adapt the abstract Newton scheme to 
an abstract evolution equation in the form 

| = Q(/). 

around a stationary solution (so Q{P) = 0). Write the Cauchy problem with 
initial datum fi ~ in the form 



$(/) := {dJ~Q{f), /(O, ■))-(0,/.). 



Starting from /°, the Newton iteration consists in solving inductively $(/" ^) + 
$'(/"-!) ■ (/" - /"-I) = forn > 1. More explicitly, writing /i" = /" - /"-i, we 
should solve 

'dth^=Q'{n-h^ 

>Ho, ■) = /.- r 
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f dth^+' = Q'in ■ h^+' - [dtr - Qin 

> 1, < 

[/i"+i(0, ■) = 0. 
By induction, for > 1 this is the same as 

/i"+i(0, ■) = 0. 

This is easily apphed to the nonhnear Vlasov equation, for which the nonhnear- 
ity is quadratic. So we define the natural Newton scheme for the nonlinear 
Vlasov equation as follows: 

j-Q _ jO(^^^ ig given (homogeneous stationary state) 
/" = /o + /i^ + . . . + /i", where 



.3) 



'dth^ + V ■ V^h^ + F[h^] ■ V^/° = 

>Ho, ■) = /.- r 



(8.4) 

Vn > 1 



/i"+i(0, = 0. 

Here F[f] is the force field created by the particle distribution /, namely 
(8.5) F[f]{t,x) = - If VW{x-y)f{t,y,w)dydw. 

Note also that all the p" = J dv for n > 1 have zero spatial average. 

8.2. Battle plan. The treatment of ( 18. 3p was performed in Subsection 14.121 Now 
the problem is to handle all equations appearing in fl8.4p . This is much more com- 
plicated, because for n > 1 the background density /" depends on t and x, instead 
of just as a consequence, 

(a) Equation (18. 4p cannot be considered as a perturbation of free transport, be- 
cause of the presence of Vt,/i"^^ in the left-hand side; 

(b) The reaction term F[h'^~^^] ■ Vvf^ no longer has the simple product structure 
(function of x)x (function of v), so it becomes harder to get hands on the homoge- 
nization phenomenon; 
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(c) Because of spatial inhomogeneities, echoes will appear; they are all the more 
dangerous that, Vt,/" is unbounded as t oo, even in gliding regularity. (It grows 
like 0{t), which is reminiscent of the observation made by Backus 

The estimates in Sections [5] to [7] have been designed precisely to overcome these 
problems; however we still have a few conceptual difficulties to solve before applying 
these tools. 

Recall the discussion in Subsection I4.11t the natural strategy is to propagate the 
bound 

(8.6) sup II/^II^a.m;! < +00 

T>0 

along the scheme; this estimate contains in particular two crucial pieces of informa- 
tion: 

• a control of = J fr dv in JF'^^+^ norm; 

• a control of (/t-) = j f^-dx in C'^'^ norm. 

So the plan would be to try to get inductively estimates of each in a norm like 
the one in fl8.6p . in such a way that /i" is extremely small as n ^ oo, and allowing a 
slight deterioration of the indices X, fi as n —>■ oo. Let us try to see how this would 
work: assuming 

VO < k < n, sup Afe,Mfc;i < 6k, 

T>0 " 

we should try to bound /i""*"^. To "solve" (18.41) . we apply the classical method of 
characteristics: as in Section [5] we define (X"^, "K^" ) as the solution of 

jX-,{x, v) = V:,ix, v), ^V:/x, v) = F[n (t, X-/X, v)) 

X"^(X, V) = X, = 

Then (18.41) is equivalent to 

(8.7) j^h-^'{t,X-„V,%x,v)) = J:-^'{T,X-,ix,v),V-^ix,v)), 
where 

(8.8) S"+i(t, X, v) = ■ V„r - Flh""] ■ V,/i". 
Integrating (18. 7p in time and recalling that /i""'"^(0, ■) = 0, we get 
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Composing with {X^q, Vj"q) and using (15.21) yields 



/i"+^(t,x,t;) 



f^Tr^'[r,Xl{x,v),V,l{x,v))dr. 



We rewrite this using the "scattering transforms" 

ni,{x, v) = {XI,, V,l){x + v{t - r), v) = SI, o 5°,; 
then we finally obtain 

(8.9) W'+\t,x,v)= [ oni,){x-v{t-T),v)dT 



t),v) dr 



X — v{t 
X — v{t — t),v) dr. 



Since the unknown /i""*"^ appears on both sides of (18. 9p . we need to get a self- 
consistent estimate. For this we have little choice but to integrate in v and get an 
integral equation on p[h"''^^] = J h"'dv, namely 



^.10) p[/i"+^](t,x 



Git 



o S^_^{x, v) dv dr 
(stuff from stage n) 



where G 



T,t 



By induction hypothesis G^^ is smooth with regularity 
indices roughly equal to A„, /i„; so if we accept to lose just a bit more on the regu- 
larity we may hope to apply the long-term regularity extortion and decay estimates 
from Section [6l and then time-response estimates of Section [TJ and get the desired 
damping. 

However, we are facing a major problem: composition of p[h'^~^^] * VW by Q^, 
implies a loss of regularity in the right-hand side with respect to the left-hand side, 
which is of course unacceptable if one wants a closed estimate. The short-term 
regularity extortion from Section E] remedies this, but the price to pay is that G" 
should now be estimated at time r' = r — bt/ (1 + b) instead of r, and with index of 
gliding analytic regularity roughly equal to A„(l + b) rather than A„. Now the catch 
is that the error induced by composition by fi" depends on the whole distribution 
not just h"", thus, if the parameter b should control this error it should stay of 
order 1 as n — > oo, instead of converging to 0. 
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So it seems we are sentenced to lose a fixed amount of regularity (or rather of radius 
of convergence) in the transition from stage n to stage n + 1; this is reminiscent of 
the "Nash-Moser syndrom" [2j. The strategy introduced by Nash [65] to remedy 
such a problem (in his case arising in the construction of C°° isometric imbeddings) 
consisted in combining a Newton scheme with regularization; his method was later 
developed by Moser [62] for the KAM theorem (see [63| pp. 19-21] for some 
interesting historical comments). The Nash-Moser technique is arguably the most 
powerful perturbation technique known to this day. However, despite significant 
effort, we were unable to set up any relevant regularization procedure (in gliding 
regularity, of course) which could be used in Nash-Moser style, because of three 
serious problems: 

• The convergence of the Nash-Moser scheme is no longer as fast as that of the 
"raw" Newton iteration; instead, it is determined by the regularity of the data, and 
the resulting rates would be unlikely to be fast enough to win over the gigantic 
constants coming from Section [71 

• Analytic regularization in the v variable is extremely costly, especially if we 
wish to keep a good localization in velocity space, as the one appearing in Theo- 
rem U2ni(iii) , that is exponential integrability in v\ then the uncertainty principle 
basically forces us to pay 0{e^^^ ), where e is the strength of the regularization. 

• Regularization comes with an increase of amplitude (there is as usual a trade- 
off between size and regularity); if we regularize before composition by VL"^, this will 
devastate the estimates, because the analytic regularity oi f o g depends not only 
on the regularity of / and but also on the amplitude oi g — Id. 

Fortunately, it turned out that a "raw" Newton scheme could be used; but this 
required to give up the natural estimate (18.61) . and replace it by the |)air of estimates 



(8.11) 



sup ||p^||jrAr+M < +00; 
r>0 



sup 

t>r>0 



fr O ^t,T 



^A{l + 6),7r;l ^ 



Here h = b{t) takes the form const. /(I -|- t), and is kept fixed all along the scheme; 
moreover A, /x will be slightly larger than A, /I, so that none of the two estimates 
in (18. lip implies the other one. Note carefully that there are now two times {t, r) 
explicitly involved, so this is much more complex than (18. 6p . Let us explain why 
this strategy is nonetheless workable. 
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First, the density = j f^dv determines the characteristics at stage n, and 
therefore the associated scattering fl"^. If is bounded in ^-^"'^+/^^ then by Theorem 
15.21 we can estimate Q^^ in Z^r , as soon as (essentially) X'^r' + fi'„ < A„, r + 
A'„ < A„, and these bounds are uniform in t. 

Of course, we cannot apply this theorem in the present context, because A„(l+6) is 
not bounded above by A„. However, for large times t we may afford A„(l+6(t)) < A„, 
while A„(l + fe)(r — bt/{l + b)) < XnT for all times; this will be sufficient to repeat 
the arguments in Section getting uniform estimates in a regularity which depends 
on t. (The constants are uniform in t; but the index of regularity goes down with 
t.) We can also do this while preserving the other good properties from Theorem 
I5.2[ namely exponential decay in r, and vanishing near t = t. 

Figure [7] below summarizes schematically the way we choose and estimate the 
gliding regularity indices. 




t 



Figure 7. Indices of gliding regularity appearing throughout our 
Newton scheme, respectively in the norm of p[hr] and in the norm 
of hj. o Qt^T, plotted as functions of t 

Besides being uniform in t, our bounds need to be uniform in n. For this we shall 
have to stratify all our estimates, that is decompose p[/"] = p[h^] + ■ ■ ■ + p[h"'], and 
consider separately the influence of each term in the equations for characteristics. 
This can work only if the scheme converges very fast. 

Once we have estimates on in a time-varying regularity, we can work with 
the kinetic equation to derive estimates on o and then on all o fi"^, also 
in a norm of time-varying regularity. We can also estimate their spatial average, in 



118 



C. MOUHOT AND C. VILLANI 



a norm C'^^^^''^'^; thanks to the exponential convergence of the scattering transform 
as r ^ oo these estimates will turn out to be uniform in r. 

Next, we can use all this information, in conjunction with Theorem 16.41 to get 
an integral inequality on the norm of p[/i""^^] in JF'^'^+'^, where A and /i are only 
slightly smaller than A„ and fin- Then we can go through the response estimates 
of Section [3, this gives us an arbitrarily small loss in the exponential decay rate, 
at the price of a huge constant which will eventually be wiped out by the fast 
convergence of the scheme. So we have an estimate on p[/i""'"^], and we are in 
business to continue the iteration. (To ensure the propagation of the linear damping 
condition, or equivalently of the smallness of Kq in Theorem 17.71 throughout the 
scheme, we shall have to stratify the estimates once more.) 

9. Local in time iteration 

Before working out the core of the proof of Theorem 12.61 in Section [TUl we shall 
need a short-time estimate, which will act as an "initial regularity layer" for the 
Newton scheme. (This will give us room later to allow the regularity index to 
depend on t.) So we run the whole scheme once in this section, and another time in 
the next section. 

Short-time estimates in the analytic class are not new for the nonlinear Vlasov 
equation: see in particular the work of Benachour [8] on Vlasov-Poisson. His ar- 
guments can probably be adapted for our purpose; also the Cauchy-Kowalevskaya 
method could certainly be applied. We shall provide here an alternative method, 
based on the analytic function spaces from Section HJ but not needing the apparatus 
from Sections E] to [71 Unlike the more sophisticated estimates which will be per- 
formed in Section [T0|, these ones are "almost" Eulerian (the only characteristics are 
those of free transport). The main tool is the 

Lemma 9.1. Let f be an analytic function, X{t) = X — Kt, fi{t) = fi — Kt; let 
T > be so small that A(t),/i(t) > for < t < T . Then for any r G [0,T] and 
any p > 1, 



K 

„A(t),M(t);p < — — II V/|LA(r),M(r);p, 

t=T ^ ^ ' 



where {d^ /dt) stands for the upper right derivative. 

Remark 9.2. Time-differentiating Lebesgue integrability exponents is common prac- 
tice in certain areas of analysis; see e.g. [32]. Time-differentiation with respect to 
regularity exponents is less common; however, as pointed out to us by Strain, Lemma 
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19.11 is strongly reminiscent of a method recently used by Chemin [T7] to derive local 
analytic regularity bounds for the Navier-Stokes equation. We expect that simi- 
lar ideas can be applied to more general situations of Cauchy-Kowalevskaya type, 
especially for first-order equations, and maybe this has already been done. 



Proof of Lemma \9.1\ For notational simplicity, let us assume d = 1. The left-hand 
side of (EU) is 



^g2.Mr)|fc| 2TTfi{T) \k\ ^ (V, + 2inkTYf{k 



n,k 



LP(dv) 



n.k 



(n- 1)! 
A"M 



{V^ + 2mkTyf{k,' 



< -i^y e^-^MI'^l 27r |A;| ^ (V, + 2^7rA;r)" f{k 



n,k 



LP{dv) 



LP(dv) 



n,k 



Kj-^2.,ir)\ki^ (V. + 2^7rA;r)"+V(^ 



LP{dv) 



< 



,27rA»(r)|fc| 



n,k 

Kt 



n,fe 



,27r^»(r)|fc| 



=1 A"(r) 



+ 



n! 



1 + r ^ 

n,fc 



LP{dv) 

{V, + 2i7rkTyVJ{k,v) 
{V^ + 2i'KkTy"v7f{k,v) 



LP{dv) 



LP{dv) 



where in the last step we used ||(V„ + 2mkT)h\\ > (1/(1 + r))(|| V„/i|| - T\\2i7rkh\\) . 
The conclusion follows. □ 



Now let us see how to propagate estimates through the Newton scheme described 
in Section [101 The first stage of the iteration {h^ in the notation of (18. 3p ) was 
considered in Subsection 14. 12^ so we only need to care about higher orders. For any 
k >1 we solve dth^^^ + v ■ Vxh^^^ = where 



S^+i = ■ VJ'' + F[f''] ■ V,/i'=+' + F[h''] ■ V,/i^') 
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(note the difference with (E7!)-(E8D). Recall that f = f + + . . . + . We 
define Xk{t) = \k - 2Kt, fik{t) = fik - Kt, where {\k)k€N, (/^fc)fcGN are decreasing 
sequences of positive numbers. 

We assume inductively that at stage n of the iteration, we have constructed 
(Afc)fc<n, {lJ'k)k<n, {Sk)k<n such that 

\fk<n, sup \\h^{t, ■ )L^fe(t).MfeW;i < 4, 
o<t<r * 

for some fixed T > 0. The issue is to construct A„+i, /i„+i and Sn+i so that the 
induction hypothesis is satisfied at stage n + 1. 

At t = 0, h"'~^^ = 0. Then we estimate the time-derivative of ,A„+i{t),M„+i(t);i. 

Let us first pretend that the regularity indices A^+i and /i„+i do not depend on t; 
then = o dr, so by Proposition Ei 

* Jo * 

t _ 

n+1 1 



and thus 







II or^n + l^/Jji + lil ^ ll^f II o--^n + l.M,i + i;l • 



Finally, according to Lemma WT\ to this estimate we should add a negative multiple 
of the norm of V/i""*"^ to take into account the time-dependence of A^+i, /i„+i. 
All in all, after application of Proposition 14.241 we get 

^ ^ (^5 ■ ) II p-^n + l(*).Mn + l(t);l ^ ||-^[^t^ ] II jr-^n + l* + Mn+l llVu/t 1 1 + 1 .A^n+l ; 1 



dt 



|-^[/r]||j^^n+i'+M„+i II Vt,/?.]; 



n+1 1 



^?i+i.Mn+i;i 



-|- ||i^[/i^ ] ||jpA„_|_it + Mn + l ll^f^J ||2'''n + l'''i+l'-' 



V/,n+l|| 



- K II v./ir+i - K II v./ir+' 



where > 0, t is sufficiently small, and all exponents A„+i and /in+i in the right- 
hand side actually depend on t. 
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From Proposition 14.151 (iv) we easily get ||jrAt+M < C || V/iH^A.^a. Moreover, 

by Proposition 14.101 



II ^Afc+i.Mfe+i; 



f-^min{A,-A„+i; 



k<n k<n 

We end up with the bound 
(1+ 

— m,"' (t, ■ ) A„+i(t),,.„+i{t);l 

at " "^t 



< 



c|y — , ^\-K 

> ^ mmjAfc - \n+l ; - /in+l} ' 



I V/l"^"^ II A„+i{t),M„+i(t);l 



+ - 



mm 



{An — An+1 ; fin ~ l^n+l} 



We conclude that if 
(9-2) E 



, ^ mini 

k<n 



linjAfc — An+i ; /ifc — /^n+i} 

then we may choose 
(9.3) Sn+l 



mm 



ill{A„ — An+i ; /in — /in+l} 



This is our first encounter with the principle of "stratification" of errors, which 
will be crucial in the next section: to control the error at stage n + 1, we use not only 
the smallness of the error from stage n, but also an information about all previous 
errors; namely the fact that the convergence of the size of the error is much faster 
than the convergence of the regularity loss. Let us see how this works. We choose 
Afc — Afe_|_i = yUfc — fik+i = A/fc2, where A > is arbitrarily small. Then for k < n, 
Afc — A,„+i > A/k"^, and therefore Sn+i < S^n^/A. The problem is to check 



(9.4) ^n^^^ 



< +00. 

n=l 



Indeed, then we can choose K large enough for (19.21) to be satisfied, and then T 
small enough that, say A* - 2KT > A", yU* - KT > where A" < A*, /i" < /i* have 
been fixed in advance. 
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If Si = 6, the general term in the series of (19 .4^ is 

^2 ^ (22)2-^ (32)2-^ (42)2-^ ...{{n- 1)2)2 

To prove the convergence for 6 small enough, we assume by induction that 6n < z"-" , 
where a is fixed in the interval (1, 2) (say a = 1.5); and we claim that this condition 
propagates if 2; > is small enough. Indeed, 

^2-" , /z(2-a)a"^2 



J n+1 

n + 1 

and this is bounded above by z"" if z is so small that 



A - V A 

la 
A 



VneN, 2(2-")-^" < 



n2 



This concludes the iteration argument. Note that the convergence is still extremely 
fast — like 0(2;"") for any a < 2. (Of course, when a approaches 2, the constants 
become huge, and the restriction on the size of the perturbation becomes more and 
more stringent.) 

Remark 9.3. The method used in this section can certainly be applied to more 
general situations of Cauchy-Kowalevskaya type. Actually, as pointed out to us 
by Bony and Gerard, the use of a regularity index which decays linearly in time, 
combined with a Newton iteration, was used by Nirenberg [68] to prove an abstract 
Cauchy-Kowalevskaya theorem. Nirenberg uses a time-integral formulation, so there 
is nothing in [HS] comparable to Lemma 19. Ij and the details of the proof of conver- 
gence differ from ours; but the general strategy is similar. Nirenberg's proof was 
later simplified by Nishida [69] with a clever fixed point argument; in the present 
section anyway, our final goal is to provide short-term estimates for the successive 
corrections arising from the Newton scheme. 

10. Global in time iteration 

Now let us implement the scheme described in Section |8l with some technical 
modifications. If / is a given kinetic distribution, we write p[f] = J f dv and F[f] = 
-VW*p[f]. We let 

(10.1) r = f + h^ + ... + h^, 

where the successive corrections /i^ are defined by the natural Newton scheme intro- 
duced in Section [SI As in Section \5\ we define as the scattering from time t to 
time r, generated by the force field = -VW * p[f'']. (Note that = Id.) 
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10.1. The statement of the induction. We shall fix p G [1, oo] and make the 
following assumptions: 

• Regularity of the background: there are A > and Cq > such that 

• Linear damping condition: The stability condition (L) from Subsection l2.2l holds 
with parameters Co, A, (the same as above) and k > 0. 

• Regularity of the interaction: There are 7 > 1 and Cp > such that for any 
u>0, 

(10.2) ||Viy*p||^,,, <Ci.||p||^.. 

• Initial layer of regularity (coming from SectionlH]): Having chosen A" < A, /i" < /i, 
we assume that for all p E [l,p] 

(10.3) VA; > 1, sup (Wh'tWz^^^.^'-. + \\p[ht]\\^.>) < Cfc, 

where T is some positive time, and Cfc converges to zero extremely fast: Cfc = 0{z'^'), 
zj < C 6 < 1, 1 < aj < 2 {fli chosen in advance, arbitrarily close to 2). 

• Smallness of the solution of the linearized equation (coming from Subsection 
14.121) : Given Ai < A^, pi < /i", we assume 



sup ||p[^r]LAi.+Mi < ^1 



(10.4) Vpe[i,p], 



r>0 



sup Ai(i+b),Mi;p < 5i, 



t>T>0 



u 



where 5i <C 5. 

Then we prove the following induction: for any n > 1, 



(10.5) VA;e {l,...,n}, VpG[l,p], 



where 



sup p[/i;] I^A,.+^, < 4 
r>0 



sup 

t>T>0 



T t,T 



< S, 



• {Sk)km is a sequence satisfying < CpCk < 4, and Sk = 0{z°- ), z < zj, 
1 < a < aj {a arbitrarily close to aj), 
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• (Afc,yUfc) are decreasing to (Aoo,/ioo), where (Aoo,/ioo) are arbitrarily close to 
(Ai,/ii); in particular we impose 

(10.6) A« - Aoo < min{l ; /i» - /i^o < min{l ; 

• T is some small positive time in (110.31) : we impose 

(10.7) x*T<t^-^. 

• b = b(t) = , where B G (0, T) is a (small) constant. 

\ ~\~ t 

10.2. Preparatory remarks. As announced in (110. 50 . we shall propagate the fol- 
lowing "primary" controls on the density and distribution: 

(10.8) (E^) VA;e{l,...,n}, sup ||^,,.+,, < 5, 

T>0 

and 



(10.9) (E;^) VA;G{l,...,n}, VpG[l,p], sup 

t>T>0 



h'l o a^;i 



Afc{l+b),Mfc;p — ^f'- 

~- bt 



Estimate (E") obviously implies, via (110. 2p . up to a multiplicative constant, 
(10.10) (E^) yk e {l,...,n}, sup < 4. 



T>0 



Before we can go from there to stage n+1, we need an additional set of estimates 
on the scattering maps {fl^)k=i,...,n, which will be used to 

(1) update the control on — Id ; 

(2) establish the needed control along the characteristics for the background 
(Vvfr) ° ^t'r (same index for the distribution and the scattering); 

(3) update some technical controls allowing to exchange (asymptotically) gradi- 
ent and composition by this will be crucial to handle the contribution 
of the zero mode of the background after composition by characteristics. 
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This set of scattering estimates falls into three categories. The first group ex- 
presses the closeness of to Id : 



Wke{l,...,n}, 



\ sup 




-Id 


1 t>T>0 




1 sup 




-Id 


1 t>r>0 







1+6 



^A*(l+i,),(M*,7) 
1+6 



<7^^(r,^), 



with Afe > a;: > Afc+i, fik> IJ^I> iJ-k+u and 



(10.12) 



,^ 2vr(A, - A*) 



min {(t — r) ; 1} 



mm 



T 



1 . 



The second group of estimates expresses the fact that f2" — is very small when 
k is large: 
(10.13) 



{Ef,) \/ke{0,...,n-l}, { 



sup 

t>T>0 

sup 

t>T>0 

sup 

t>T>0 



_A* (l+6),{(j* ,7) 
1+6 



<2 7^fXr,^), 



.*a+6w . <7^^"(r,t)+7^r(r,t), 

^A* (1 + 6),(m„,7) — 1 \ ' / ^ V ' 



) fc.n / 



6t 
1+6 



J. A* (1+6), n 



with 



(10.14) 



7^^"(r,t) 



27r(A,- - A* 



min {(t — r) ; 1} 



mm 



(Choosing /c = brings us back to the previous estimates (E[^).) 

The last group of estimates expresses the fact that the differential of the scattering 
is uniformly close to the identity (in a way which is more precise than what would 
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follow from the first group of estimates) 
(10.15) 

Wn^'Xt.r - (/, 0) 

I t>T>0 

(E^J VA; = l,...,n 



z 



<2 7^^(r,^), 



sup 

t>T>0 



VQ!'Vt,r - (0,/) < 7^^(r,^) + 7^^(r,^), 



where V = (V^, V„), and / is the identity matrix. 



An important property of the functions 7l\'"'{T,t), T^''^{T^t) is their fast decay as 
r — * oo and as /c — > oo, uniformly in n > /c; this is due to the fast convergence of 
the sequence {5k)k&- Eventually, if r G N is given, we shall have 



(10.16) 



Vr>l, 7^^"(r,t)<c.^;:^(r,t), 7^^'"(r, t) < c.^J(r, t) 



with 



5, 



min {(t — r) ; 1} 



and 



r,2 
k,n 



min{(t-r)V2;l} 



for some absolute constant depending only on r (we also denote cOq]^ = uj^^ and 



r,2 r2\ 
^0,n = J 
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From the estimates on the characteristics and (E[J) will follow the following "sec- 
ondary controls" on the distribution function: 



(10.17) (EJ^) VA;G{l,...,n}, Vpe[l,p], 



sup 

t>T>0 

sup 

t>T>0 



(V. + rV.)(/i(!o(]J;i) 



A(.(l + 6),f.j.;p 



< 4 



Aj.(l+6),Mfc;p 



bt 
1+6 



^Aj.(l+6),Mfe:p 
1+6 



A;,(l + 6),Mfe;p 

r- M 
T+T 



< 5t 



sup 

t>T>0 



1 + r 



(vv/i^ o n 



k-l 



sup (1 + r) 

i>r>0 



r^fe(i+i').Mfc;i 
^__6<_ 

1+6 



<5a 



^Afc(l+6),Mfea 

"r-Jitr 
T+E 



<4. 



The transition from stage n to stage n + 1 can be summarized as follows: 



(An) 



(E^^) + (E^) + (E^„) + (El 



iK) 



(Bn) 



:e 



n+l^ 



fE 



n+l^ 



(E^^) + (E 



The first implication (An) is proven by an amplification of the technique used 
in Section O ultimately, it relies on repeated application of Picard's fixed point 
theorem in analytic norms. The second implication (Bn) is the harder part; it 
uses the machinery from Sections [6] and [3, together with the idea of propagating 
simultaneously a shifted Z norm for the kinetic distribution and an T norm for the 
density. 

In both implications, the stratification of error estimates will prevent the blow up 
of constants. So we shall decompose the force field generated by /" as 

pn ^ p^jn^ = + . . . + E", 

where E'^ = F[h''] = -VW * p[h'']. 
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The plan of the estimates is as follows. We shall construct inductively a sequence 
of constant coefficients 

A« > Ai > A* > A2 > . . . > A„ > a; > A„+i > . . . 

/i" > /il > yU* > /i2 > . . . > yU„ > /i* > /i„+i > . . . 

(where A„, will be fixed in the proof of (An), and \n+i, ^n+i in the proof of (Bn)) 
converging respectively to Aqo and fioo] and a sequence {6k) ken decreasing very fast 
to zero. For simplicity we shall let 

7^"(r,t) =7^^(r,t)+7^^(r,^), 7^'='"(r,^) =7^^"(r,^)+7^2''"(r,^), 

and assume 27r(Aj — A*) < 1; so 
(10.18) 




^'>- "^1 (27r(Aj - A*))2+^ y (1 + r)- 

It will be sufficient to work with some fixed r, large enough (as we shall see, r = 4 
will do). 

To go from stage n to stage n + 1, we shall do as follows: 

- Implication (An) (subsection 110.31) : 

Step 1. estimate fi" — Id (the bound should be uniform in n); 

Step 2. estimate f2" — fi'^ (/c < n — 1; the error should be small when k 00); 

Step 3. estimate Vl^" - /; 

Step 4. estimate (fi'')"^ o f]"; 

- Implication (Bn) (subsection 110. 4p : 

Step 5. estimate h'' and its derivatives along the composition by fi"; 
Step 6. estimate using Sections [H] and [3 

Step 7. estimate from p[/i"+^]; 

Step 8. estimate o fi"; 

Step 9. estimate derivatives of composed with 

Step 10. show that for h"'^^, V and composition by fi" asymptotically commute. 

10.3. Estimates on the characteristics. In this subsection, we assume that esti- 
mate (E") is proven, and we estabhsh (E^) + (EJ^) + (E^j^). Let A* < A„, yU* < /i„ 
to be fixed later on. 
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10.3.1. Step 1: Estimate of Vt^ — Id . This is the first and archetypal estimate. We 
shall bound — x in the hybrid norm 2'^^^^^+'')'^". xhe Sobolev correction 

7 will play no role here in the proofs, and for simplicity we shall forget it in the 
computations, just recall it in the final results. (Use Proposition 14.321 whenever 
needed.) 

Since we expect the characteristics for the force field to be close to the free 
transport characteristics, it is natural to write 

(10.19) Xl^{x,v) =x-v{t-T) + Zl^x.v), 

where Z"^ solves 



92 



■Zl^{x, v) = (r, x-v{t-T) + Zl^{x, v)) 



(10.20) { dr'^' 

Z"(x,t;) = 0, d^Zl^ {x,v)=Q. 



t=T 



(With respect to Section [5] we have dropped the parameter e, to take advantage of 
the "stratified" nature of F^; anyway this parameter was cosmetic.) So if we fix 
t > 0, (Z"t-) is a fixed point of the map 



defined by 
(10.21) 




0<T<t I ^ 0<T<t 

Zt,r = F^(T,X-V{t-T) + Wt,r) 
t,t = 0, QrZt' 



= 0. 

T=t 



,A;(1+6),/^* 
-'t- 



The goal is to estimate Z"^ — a; in the hybrid norm " \t 
We first bound {Z^)t,r = ^(0). Exphcitly, 

{Z^)tA^, v) = j^{s-r) x-v{t- s)) ds. 
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By Propositions 14.151 (i) and 14.191 



(10.22) 



[[(■^O )t,-r|| 7^n(l + '')>A'n 



bt 



< ^ (s-r) F"(s, x-w(t-s)) 



^ bt 



■', bt 
1 + b 



(s-t) \\F^{s,-)\\^.,s,) ds, 



where 



(10.23) 



uis,t) = \:\s-bit~s)\+fj.i 



First case: If s > bt/{l + b), then 



(10.24) 



u{s, t) < A* s + yu* < Afe s + /ifc - (Afc - A* )s (1 < A; < ra) 



Second case: If s < bt/ (1 + 6), then necessarily s < B <T . Taking into account 
(110. 6p . we have 



(10.25) 
(10.26) 



z/(s, t) = \lbt + ^l- A* (1 + b)s 
<\lB + fil-{\k-\l)s. 



(Of course, the assumption A" — Aoo < niin{l, Aoo/2} imphes A^ — A* < A*.) In 
particular, by (110.71) . 



(10.27) 



v{s,t)<^i^-{\k-\l)s {l<k<n). 
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We plug these bounds into (110.221) . then use E'^^s, 0) = and the bounds (110. lOp 
and (110.241) (for large times), and (110. 3p and (110. 27p (for short times). This yields 

(10.28) 

II (^0 )t,-r|| ,A*(l+i,),M^ 
bt 
1+6 

-Y.\ I " ^) ll^''^^' ■ )II^A,.+,,-(A,-A*). ds 

k=i V-^-^m 

I 1 \ J rVv-i-r 



k=l \-^'^'^T+b 

+ j {s-T)e~^-^^'^~^*-)s\\E\s,-)\\^,,ds 

<^^k {s-t) e-^-C^'c-^n)^ ds 
k=i -^^ 

< |:4e---=.' mi„{(^l)!; J^^^J^} < K(r.t). 



Let us define the norm 



0<T<t 



sup 

0<T<t 



^A„(l + b),n„ 
bt 
1 + 6 



7^^5(r,^) 



(Note the difference with SectionO now the regularity exponents depend on time(s).) 
Inequality (110.281) means that ||||^(0)||||„, < 1. We shall check that ^ is (1/2)- 
Lipschitz on the ball i?(0,2) in the norm |||| ■ ||||„. This will be subtle: the uniform 
bounds on the size of the force field, coming from the preceding steps, will allow to 
get good decaying exponentials, which in turn will imply uniform error bounds at 
the present stage. 
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So let W,W e B{0, 2), and let Z = ^{W), Z = ^{W). As in Section [3 we write 



1 




Jr 



(s - r) V.F'^ [s, x-v{t-s)^ [e Wt,s + (1 - ^) H^i,. 



and deduce 



where 



0<T<t 



0<s<t 



A{t) = sup I) X 



/V'^ - r) II V,.F" (s, X - - s) + (0 Wt,, + [l - 0) W*,,) 



yA* (l + b),n* 
bt 



dsde. 



For r < s we have 7^.2 < 7^.2 (T,t). Also, by Propositions 14.251 (applied with 
y = 0, 6 = — (t — s) and a = in that statement) and 14.15] 

A{t)< sup f (s-r) ||V^F"(s, ■)II.F-(»,t)+e(»,t)C?s, 

0<T<t Jr 

where v is defined by (110. 23p and the "error" e(s, t) arising from composition is given 
by 

~ " <27^^(s,^). 



Since 



e(s,t) = sup eWt,s + {i-e)Wt. 

o<e<i 



5, 



■"f bt 

T+b 



mm 



in{(t-.)V2; 1} 



(1 + ^) 



we 



have, for all < s < t, 



(10.29) 

as soon as 
(10.30) 



2 7?.2 (s, t) < ^h{t - s) ls>bt/{i+h) H ^ — ^ ls<w/(i+6) , 



(Ci) Vn>l, 2Cl [Y, 



6, 



- (27r(A. - A*))^ 



< min 



6 ' 2 
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We shall check later in Subsection 110.51 the feasibility of condition (Ci) — as well 
as a number of other forthcoming; ones. 

The extra error term in the exponent is sufficiently small to be absorbed by what 
we throw away in fll0.24p or f 1 1 . 2 5 p - (1 1 . 2 6 1) - (1 1 . 2 71) . So we obtain, as in the estimate 
of Zq, for any k E {1, . . . , n}, 



(z/ + e)(s,t) 



'< XkS + fik- {Xk- Xn)s ioT s>bt/{l + b) 

< - (Afc - a;) s for s < 6t/(l +6), 

and we deduce (using ( 110. lOp and 7 > 1) 



Ait) <snpJ2\f 



[s - t) \\V.j,E^{s, ■ )||^Afe3+^,-(A,-A*). ds 

^ 1+6 



+ j (s-t) \\V^E'{s, ■ )iU.«-(A,-A*> ds 
< sup V^fc / is-T)e-^^^-^"^'ds< sup 7^^(r,t) =7^^(0,^) 

0<r<f Jt 0<T<t 

n j- 

-^(27r(A,-A;))2- 

If the latter quantity is bounded above by 1/2, then \i/ is (l/2)-Lipschitz and we 
may apply the fixed point result from Theorem IA.2I Therefore, under the condition 
(whose feasibility will be checked later) 

we deduce 



bt 
1+6 
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After that, the estimates on the scattering are obtained exactly as in Section [5j 
writing = {fl^Xt^T-, ^^Vt^r), recalhng the dependence on 7 again, we end up with 



(10.32) 




^A*{1+6),(m?;,7) 
1+6 



A* (l+6),(M,*i,7) 
1+6 



<2 7^-(r,^) 

<7^^(r,^). 



10.3.2. Step 2: Estimate of Q'' - . In this step our goal is to estimate fi" — ^2 
for 1 < k < n — 1. The point is that the error should be small as k —>■ 00, uniformly 
in n, so we can't just write — < — Id || + — Id ||. Instead, we start 
again from the differential equation satisfied by Z'^ and Z": 

= F" (r, x-v{t-T) + Zl^x, tO) - (^r, a; - - r) + 
[t, X - - r) + - F" (r, x - v{t - r) ^ Zl 
+ [F--F')(T,x-v{t-T) + Zl). 

This, together with the boundary conditions Z^^^ — Z^^^ = 0, dr{Z^.^ — Z^^^)\r=t = 0, 
implies 



I ^) VxF'^(5,a; - v{t - s) + {6 Z^, + {1-6) Zl, 

+ (s - r) (F" - F^) [s, x-v{t-s) + Z^^ix, t;)) ds. 



[Zl-Z')dsde 



We fix t and define the norm 



(^t,r)o<T 



<t 



sup 



bt 

k,n I 



„ o<r<i 7^.2'" (r,t) 
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where 7^.2'"' is defined in (110.141) . Using the bounds on Z^^Z^ in 

nil ■ nil < nil 



since 



(10.33) [Zl, - Zl) 



^ by using the fact that Ti\ < TZ^) and proceeding as before, we get 

1 



+ 



0<r<t 
t 



< 

k,n 2 



0<T<t 



k,n 



(^j\s - t) (F" - F') X - v{t -s) + Zl) ds^ 



0<T<t 



k,n 



Next we estimate 

{F''-F'')(^s,x-v{t-s) + Zl 



bt 
1 + 6 



iF--F')is,Xl) 

{F- -F'){s, nt) 



■'t bt 



^A*(1+6),M* 

bL. 

1+6 



< 



(F" -F^)(., ■) 



jpi/(s,t) + e(s,t) ' 



where the last inequahty follows from Proposition 14.251 u is again given by (110.231) . 
and 



e{s,t) = \\n''Xt,s-ld 



.A*,i+6),< <2 7^^(s,^)<2 7^^(s,^). 

^'-i^b 



The same reasoning as in Step 1 yields, under assumptions (Ci)-(C2), for k+1 < 
3 < n: 



(z/ + e)(s,t) 



<XjS + fij - {Xj - X*J s 
<//tt-(A,-A:)s 



for s > ht/{l + h) 
for s < bt/{l + b), 



and so 



< 



1 1 rpn rpk 

W''^ S ''^ S \\j:vJr£ — / J "3 

j=k+l 

For any r > 0, by integrating in time we find 



-2t(Aj-A*)s 



J^{s-r) (F" - F') (s, x-v{t-s) + Zl) ds 
(s-r) 5Je-2-(^^-^")^rfs<7^2•"(r,^). 



^A*(l+6),M?; 
■'i__M_ 

1+6 



j=fc4-l 
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Therefore 



(^l\s - t) (F" - F') (s, x-v{t-s) + Zl)i ds^ 
and by ffTOSS]) 



0<r<t 



< 1 



k'n 



0<T<t 



< 2. 



k,n 



Recalling the Sobolev correction, we conclude that 
(10.34) - 1^'X,,JL < 27^2''"(r,^). 



bt 
' 1+6 



For the velocity component, say U, we write 
d 



dr 



= F''{t, X - v{t - r) + Zl^{x, tO) - (t, x - v{t - t) + Z^^ix, t;)) 

= F-(t,x - v{t - t) + Z^^^) - F'^{t,x - v{t - r) + Zl) 

+ {F--F'){T,x-vit-r) + Zl,y 

where Z", Z^ were estimated above, and the boundary conditions are U^^^ — U^^ = 0. 
Thus 

^11 ^^'^^ ^ - ^'(^ - ^) + ^Is + (1 - ^) zi) ) ■ (z;:, - zl^) ds de 

+ j^{F''- F'') (^s, x-v{t-s) + Zl^ix, tO) ds, 
and from this one easily derives the similar estimates 



bt 
" 1+6 
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10.3.3. Step 3: Estimate ofWil"". Now we establish a control on the derivative of 
the scattering. Of course, we could deduce such a control from the bound on f2" — Id 
and Proposition I4.32( vi) : for instance, if A** < A*, /i** < /i*, then 

(10.35) - < 



■r M_ 

1 + 6 



min{A*-A;*; /i* - /x**} 



But this bound involves very large constants, and is useless in our argument. 
Better estimates can be obtained by using again the equation (110.201) . Writing 

(ni^ -ld){x,v)= [Zl^ {x + t;(t - r) , v) , Z^^ [x + v{t ~ t),v)), 

where the dot stands for 8/ Or, we get by differentiation 

V.ni^ - (/, 0) = (V.Z,", {x + v{t-T),v), VJl, {x + v{t-T),v)), 

V,QI,-{0,I) = (^{VMt-r)V^)Zll^{x+v{t-T),v), {VMt-r)'s/.)Zl[,{x+v{t-r),v) 

Let us estimate for instance V^f^— (/, 0), or equivalently VxZ^^. By differentiating 
ffnOnH . we obtain 



dr 



VxZl^{x, v) = V,F" (r, X - v{t - r) + Z^^ix, v)) ■ (Id + V.,Zl^). 



So V xZ^^ is a fixed point oi^ -.W i — > Q, where W and Q are functions of r G [0,t] 
satisfying 

^ = V.F" {r,x-v{t-T) + Zi:) (/ + ) , 
Qit) = 0, drQ{t) = 0. 

We treat this in the same way as in Steps 1 and 2, and find on Qx (the x component 
of Q) the same estimates as we had previously on the x component of Q. For 
the velocity component, a direct estimate from the integral equation expressing the 
velocity in terms of F yields a control by TZ"^ + • Finally for V„f2 this is similar, 
noting that (V^ + {t — 'T')'Vx){x — v{t — r)) = 0, the differential equation being for 
instance: 



dr' 



■(V,+(t-r)V,.)^;:.(x,tO = V,.F"(r,x-i;(t-r)+Z[:,(x,i;))-((V,+(t-r)V,.)^t':. 
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In the end we obtain 



(10.36) 



sup 

t>T>0 

sup 

t>T>0 



-(/,o) 
(0,/) 



J.A* (l + b),/j* 

\--bt_ 
1+6 



<nUT,t) + nUT,t). 



10.3.4. Step 4- Estimate of {Vt^)'^ o We do this by applying Proposition 14.281 
with F = Vt^ , G = . (Note: we cannot exchange the roles of VL^ and 1)" in this 
step, because we have a better information on the regularity of Vt^ .) Let e = e{d) 
be the small constant appearing in Proposition 14.281 If 



(10.37) 

then \\S/VlI^ - 
(10.38) (C4) 



then 



(Cs 



VA;>1, 3 7^^(r,^) +7^^(r,t) < £, 



/|| \l(l+b).,^,l < e; if in addition 



-bt/{l + b) 

VA; G {l,...,n-l}, Vt>r, 
2(1 + r) (1 + B) (3 7^2''" + 7^^") (r, t) < max{A^ - A; ; /i^ - /i;}, 



■ A:(1 +b)+2\\Q^- (^1L < A*(l + b) 



/i: + 2 1 





bt 






" 1 + 6 



bt 
1+6 

— fi*"'!! A*(i+b),n* < 

^~T+T 



(Once again, short times should be treated separately. Further note that the need for 
the factor (1 + t) in (C4) ultimately comes from the fact that we are composing also 
in the v variable, see the coefficient a in the last norm of (14.301) .) Then Proposition 
081 (ii) yields 



(f]t)"'°^i:.-M 



< 2 Wo, 

^A*(1+6),M^ — II t,T 



bt 
■ 1+6 



1+6 



<4(7^l'" + 7^2'")(r,^). 



10.3.5. Partial conclusion. At this point we have established (Ef^) + (E^) 



fE 



10.4. Estimates on the density and distribution along characteristics. In 



this subsection we establish (E""^-*- 
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10.4.1. Step 5: Estimate of o and {Vh^) o 1]" (k <n). Let A; G {1, . . . , n}. 
Since 

the control on o f2" will follow from the control on o VL^~^ in (E|J), together 
with the control on {Vt^-^Y^ o fi" in (E[?j). If 

(10.39) (l + r)||(f],^;i)-iof];;^-Id||^.*a+,)„* <min{(Afc-A:); 

then we can apply Proposition 14.251 and get, for any p G [l,p], and t > r > 0, 



(10.40) 



T t.T 



\--kL. 

1+6 



< 



6t 
1+6 



In turn, (110.391) is satisfied if 

(10.41) (Cs) VA;G {l,...,n}, VrG[0,t], 

4 (1 + r) (7^^"(r, t) + 7^^'"(r, t)) < min{ - A^ ; - 

we shall check later the feasibility of this condition. 
Then, by the same argument, we also have 



VA;G{l,...,n}, VpG[l,p] 



sup 

t>r>0 



J.A* (l + i)),^i* ;p 
bt 

^-t+t; 



^ bt 



10.4.2. Step 6: estimate on p[h"'~^^]. This step is the first where we shall use the 
Vlasov equation. Starting from (18.41) . we apply the method of characteristics to get, 
as in Section [HI 

(10.42) h''+'{t,X^/x,v),V,-/x,v)) = rS"+^(r,Xo",(x,^),yj:,(x,^))rfr, 

Jo 

where 

We compose this with (X"q, V^"q) and apply (15. 2p to get 

h^+\t,x,v) = f E^+'{t,XI^^^{x,v),V,1{x,v)) dr, 



140 C. MOUHOT AND C. VILLANI 

and so, by integration in the v variable, 

(10.43) p\h^+\t,x)= ! [ J:"+^{T,X]^^^{x,v),V^';^{x,v))dvdT 

Jo Jw 

* / {K+'-Gl,){x-v{t-T),v)dvdT 
-'R'' 

* / {Rl,-H:^^,){x-v{t-r),v)dvdT, 

-'R'' 

where (with a shght inconsistency in the notation) 

PI, = (V./") o ni^, = (v./i") o ni^. 

Since the free transport semigroup and fi"^ are measure-preserving, 
VO<r<t, [ [ {R'^f ■Gl,){x-v{t-T),v)dvdx 

R^^t^ ■G':^tdvdx 



;i0.44) 




j J F[K'+^]-V^rdvdx 

j y ■ (F[/i"+^] P) dvdx = 0, 



and similarly 

VO<r<t, [ I {R^,- H'^,){x-v{t-T),v)dvdx = 0. 

This will allow us to apply the inequalities from Section [61 
Substep a. Let us first deal with the source term 

(10.45) a"'"(t,x) := y j {R^, ■ H^,){x - v{t - t),v) dv dr. 
By Proposition 16. 2[ 

(10.46) ||cr"'"(^, ■ )||jpA*t+M;^ < / ||-Rr,J2^"(i+'')'^" ll-^T,J2^"''+''>'''"'' 

/n ^ bt ^ bt 
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On the one hand, we have from Step 5 



\Hr,t\\ ^yi,{l+b),nl,;l < 2 (1 + r) 5n- 



On the other hand, under condition (Ci), we may apply Proposition 14.251 (choos- 
ing (T = in that proposition) to get 



where 



^/„(t,r) =/i; + A;(l + 6) 



</i: + A:(l + 6) 



bt 



1 + 6 
bt 



' bt 



1+6 



1 + 6 



2 7^^(r,^). 



Proceeding as in Step 1 (treating small times separately), we deduce 

^ bt 



with 
(10.47) 



' Un{T,t) := fi^ when < r < 6t/(l + 6) 

UniT.t) := XnT + fin whcU T > 6t/(l + 6) . 



(We have used the gradient structure of the force to convert (gliding) regularity into 
decay.) Thus 



(10.48) 



^ 1+6 1+6 



<2Cf61 f e-2"(^"-^")"(l + r)rfr 



< 



(vr(A„-A*))2 

(Note: This is the power 2 which is responsible for the very fast convergence of the 
Newton scheme.) 
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Substep b. Now let us handle the term 

(10.49) a"'"+'(t,x) ■■= j [Kf ■ Glt){x - v{t - T),v)dvdT. 

This is the focal point of all our analysis, because it is in this term that the self- 
consistent nature of the Vlasov equation appears. In particular, we will make crucial 
use of the time-cheating trick to overcome the loss of regularity implied by composi- 
tion; and also the other bilinear estimates (regularity extortion) from Section [6], as 
well as the time-response study from Section [71 Particular care should be given to the 
zero spatial mode of G", which is associated with instantaneous response (no echo). 
In the linearized equation we did not see this problem because the contribution of 
the zero mode was vanishing! 

We start by introducing 

n 

(10.50) = V./° + V. {h'; o , 

k=l 

and we decompose a^'^+^ as 

(10.51) cx".«+i = + ^ + ^, 
where 

(10.52) a"'"+^(t,x) = ^ j F[h^+^]-Glt{x-v{t-T),v)dvdT 
and the error terms £ and £ are defined by 

(10.53) £{t,x) = J ((f[/i^+1]o1][;^-F[/i;^+1]) •G")(r,a;-t;(t-r),i;)di;dr, 

(10.54) £{t,x) = j {f[K'+^]- {G'' -G")^{T,x-v{t-T),v)dvdT. 
We shall first estimate £ and £. 
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Control of S: This is based on the time-cheating trick from Section El and the 
regularity of the force. By Proposition 16.21 



(10.55) < / F[K+']oni,-F[K'^'] 

Jo 



""^ bt 



7A*(i+b),,,* ;i dr. 



bt 
1+b 



From fllU.ll) and Step 5, 

n 

(10.56) ||G'"||^A*{i+6),M*;i < ||V.u/° O r2"^||^A*(i+b),M*;i + V.t,/l^ O fi"^||^A*(l+(,),M*;i 



bt 



k=l 



<Co+ E^M ^^ + ^) 



where C'q comes from the contribution of f^. 

Next, by Propositions S^l and i^S] (with y = 0, r = cr, 6 = 0), 

(10.57) 



n+ll 



,A* (l + b),n* 

\--bL. 
1+b 



< 



Jo 



7A*(1+6),M* 
bt 



dO I W^t^^ — Id ||2^R(i+b).M?i 



bt 

' 1 + b 



< WF[hr']\\^.. ||fii:,-Id||,A, 
where 



6t 



1 + 6 



^^Xt,T — X\ \ „A* • 
1+6 



Small times are taken care of, as usual, by the initial regularity layer, so we only 
focus on the case t > bt/{l + b); then 

< {Kr + /^:) - a: bit - r) + 2 7^"(r, t) 

min{t — r ; 1} 



1 + r 
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To make sure that Un < + yU* , we assume that 
and we note that 

minjt — r:l| ft — t 



1 + r - \l + t 

(This is easily seen by separating four cases: (a) t < 2, (h) t > 2 and t — r < 1, (c) 
t > 2 and t - r > 1 and r < t/2, (d) t > 2 and t - r > 1 and r > t/2.) 
Then, since 7 > 1, we have 



[10.59) \\VF[K+']\\^^^ < \\VF[K+' 



(Note: Applying Proposition 14. 101 instead of the regularity coming from the interac- 
tion would consume more regularity than we can afford to.) 
Plugging this back into (110.571) . we get 



7^n(l + i').Mn 
1 + 6 



<2 7^'^(r,t)C^||p[/^^l]||^.,.+,, 



Recalling (110.531) and (110.561) . applying Proposition 14. 24^ we conclude that 
(10.60) \\S{t,.)\\^,,,,^,<2ClCF (c', + J2A (E 



;^(27r(A.-A^))5^ 

\\p[K^ 



k=l / \fc= 

dr 



(We could be a bit more precise; anyway we cannot go further since we do not yet 
have an estimate on p[h"'~^^].) 
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b2. Control of S: This will use the control on the derivatives of h^. We start again 
from Proposition I6.2t 



(10.61^ 



J.A* (l + 6),^i* ;1 
bt 



\F[K 



n+ll 



where /?„ = A* (1 + 6) |r — + We focus again on the case r > + 

so that (with crude estimates) 

and the problem is to control — C": 



(10.62) IIG^-G' 



1+6 



1 < 



(V,J°) o - v./^ 



^A*(l+6),M?j;l 

1+b 



fc=i 



bt 
-T+B 



^A*(l+6),M?j;l' 



fc=l 



bt 

-T+B 



By induction hypothesis (E|J) , and since the Z^'^' norms are increasing as a func- 
tion of A, /i, 

It remains to treat the first and second terms in the right-hand side of (110.621) . 
This is done by inversion/composition as in Step 5; let us consider for instance the 
contribution oi , k > 1: 

r- M 



_A*(l+(,),M,*i;l 

^ bt 



„A*(l+b),M* '^^ 
^ bt 



< 2 



^A*(l+f,).M*;l 



it 
1+6 



bt 
' 1 + b 



<4 5fc(l + r)'7^'^-^'"(r,^) 



1 + r 
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where in the but-to-last step we used (EJ^), (E"), Propositions 14.24 1 and 14.281 Con- 
dition (C5) and the same reasoning as in Step 5. 

Summing up all contributions and inserting in (110. 6ip yields 



(10.63) \\S{t 



< 4Cf 



k=l 



1 (2ir(A, - A;.))< 



k=l 



* n+1 (^'^ 

\\P[K ]lljC-A*r + M* —— 77 

[1 + Ty 



b3. Main contribution: Now we consider 0""'"+^^^ which we decompose as 



— n,n+l — n,n+l i \ ^ — n,n+l 



k=l 



where 



— n,n+l/ \ 




(r, x-v{t-T),v)- VJ%v) dv dT, 



j {F[h^+^]-V,{h';on\-^)^{T,X-v{t-T),v)dvdT. 



Note that their zero mode vanishes. For any A; > 1, we apply Theorem 16.41 (with 

M = 1) to get 



-n,n+l 1 1 

t,k ||jF^n*+M' 



,< f K-''\t^T)\\F[h 
Jo 



+ fK-^'\t,T)\\F[K^']\\^,^^^dT, 
Jo 



where 



bt 



T — 



K^'^ {t, t) = sup 

0<T<t 



K^'^{t,T) = (I+t) sup e 

i^O, m^O 

\0<T<t 



1 + b 



■"T-bt/a + b) 



1 + r 



T^n,k 
^1 ) 



-2lT 



1 + \i-m\i 

n,k 



\i{t-T)+mT\ 



K, 



! 
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• Kq' {t, t) = e 

We assume 

(10.64) 



t - T 

TTt 



7] > small, 



and check that u'^ < A*r + /i*. Leaving apart the small-time case, we assume 
r > bt/{l + b), so that 



J^n = [Kt + ^nj ^ 

which is indeed bounded above by A* r + yU* as soon as 

(10.65) V<BX*^. 
Then, with the notation (17.101) . 

(10.66) K'^'\t,T) < Kf"''^'^(t,r), 
with 



(10.67) 



an,k = 27r min 



From the controls on (assumption (EjJ)) we have 



^Afe(l+b),Mfc;l 
1+6 



< 



^A(.(i+6),Mfe;i 

1+6 



and 



CAfc{l + 6);: 



< 



<5fc(l + r); 

<(V. + rV.)(/i^fi,V)> 
(V. + rV.)(/i^ofif;i) 



cAfc(i+6)a 

Afe(l+6)a 



^ 6t 
^-T+7^ 



<Sk. 
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After controlling ] by p[/i"^ ], we end up with 

(10.68) \\K;ri^K^.., <CF 



n+l] II 



dr 



t / " 



,fe=i 



fee 



\p[h 



n+l-| II 

T J lljf A*T + M* ) 



(ir, 



with an,k defined by (110. 67p . 



Substep c. Gathering all previous controls, we obtain the following integral inequality 
for 9? = 



(10.69) 



j {VW *ip){T,x -v{t-r)) ■VJ%v)dvdT 



<An + 



where, by mM^ . (fTO^OD and (fTO^SD . 
(10.70) A„ = sup ||a"''^(t, •)|U.,.+.,* 



< 



(l + r)2 



• )ll:F^?i-+M?i, 



fc=i 



n 

Ko"(t,r) = Cp5^4e 



q;„ = q;„ „ = 27r min 



vr(A.-A*))2' 



-27r 



^)(t-r) 



k=l 



4 



fe=i 



^ (2,r(Ai - AJ))* 



(We are cheating a bit when writing (110.691) . because in fact one should take into 
account small times separately; but this does not cause any difficulty.) 
We easily estimate K^: 



/ K^it,T)dT<CFj2 

Jo 



Sk 
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/27r(A, - A*) 

Let us assume that a.„ is smaller than 0(7) appearing in Theorem 17.71 and that 

(10.71) (c.) 3c,ci fc; + x: 4 + 1) (± < ' 



k=l 



-(27r(A,-A*))f 



4- 



(10.72) (Cg) CfY^= ^'' <-, 

(10.73) (C) C-.X:-^^<max|-;x|. 

(note that in these conditions we have strenghtened the inequalities by replacing 
Afc — A* by \k — \l where x > is also defined by Theorem 17. 7p . Applying that 
theorem with Aq = A, A* = Ai, we deduce that for any e G (0, a„) and t > 0, 

\2 



(10.74) ||pj^+i||^,,,+,, <CA„^i±2^e^^o (1 + ^) e^^-e^^"(^+^-)e^* 



where 



and 



Tp ri — Cry max 



,fc=i 



Cl ( C„ (c^)2/3 



Pick up Ajj < A* such that 27r(A* — AjJ < q;„, and choose e = 27r(A* — Aj^); recalling 
that p""''^(t, 0) = 0, and that our conditions imply an upper bound on c„ and Cq, we 
deduce the uniform control 

(10.75) llpM^4...?i < e-^-(^'-^'^^)* IIpM^.,.... 

< CA„ ( 1 + ^ ^ e^^", 

^.(A^-Al)^^; 
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where 
(10.76) 



7-1 



10.4.3. Step 7: estimate on F[h"'~^^]. As an immediate consequence of (110. 2p and 
fll0.75p . we have 



(10.77) 



sup ||F[p"+^]|| ,t,, . < CAn ( 1 , 



10.4.4. Step 8: estimate of h^~^^ o In this step we shall use again the Vlasov 
equation. We rewrite (110.421) as 



Jo 

but now we compose with (X"q, V^"q), where t > t is arbitrary. This gives 
h-+'{T,X^^^{x,v),V-,{x,v)) = [\-+'{s,Xl{x,v),V-{x,v))ds. 



Then for any p E [l,p] and < X^, using Propositions 14.19] and ITMt and the 
notation (110. 44p . we get 



y(l + 6)Atj,,.* ;p 



^(l + b)Ajj,(i* ;p 
't bt 



< 



s7 11 ^(l+6)A»,,M?i;p 

-'a bt 



(is 



<C I llR^^-'-ll llr^" II <7e 

— / II ^(l+6)AtMj ll'-^S,tll 7(l+b)A?i,M?;;p "-^ 



bt 



6t 



l|-^S,J2(l + b)A^,Mn \\^S,t\\ ^{l + b)>.i,^*^;p . 



bt 
' 1+6 



bt 
' 1+6 



Then (proceeding as in Step 6 to check that the exponents lie in the appropriate 
range) 



1+6 
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and 



with 



' Un{s,t) := fi^ when s <bt/{l + b) 

Pn{s, t) := Ajj s + /i* when s > bt/ {1 + b). 

On the other hand, from the induction assumption (EJ^)-(EJ^) (and again control 
of composition via Proposition 14.251 . . ), 



S,t\\ ^(1+6)A?,,mS. 
bt 

"1+5 



< 2 (1 + .) 5„ 



and 



We deduce that 



satisfies 



\\G[ 



s,t\\ _(i+6)a! 

1+6 



e-2-(^"-^")^ ||p:+i^.„(.) (1 + s) ds 



2CF6i / e 
/o 



-27r(At-A^)s 



(1 + s) ds; 



so 



(10.78) Vt>r>0, 

ll^r"^ ° ^t,r|L(i+6)A^,Mj;P — 



4Cf max{(^^^^4) ;!} / .2 , n n+i 



6f 
1+6 



(2vr(At-A^))' 



5„ + SUp||p^ 

s>0 
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10.4.5. Step 9: Crude estimates on the derivatives of h^'^^ . Again we choose p G 
From the previous step and Proposition 14.271 we deduce, for any A]^ such that 
< >t< ^i, and any fil < fi^, 



(10.79) v^K^'oni^) 



^A* (l + 6),Mn;P 
1+6 



< 



+ 



;.A* (l + i)),nJi;p 
1 + 6 



C{d) 



mm 



t,T\\ „A?,(l + 6),Mn;P' 



and 
(10.80) 

Similarly, 
(10.81) 



^A* (1+6),m|i;p 
\--bL. 
1+6 



< 



C{d){l + r) ||,n+1 on II 

TTlinlA^ A* ■ //* //H °"t,r||^A^,(l+6),M?i;p- 

mm |A„ - An , /i„ - /in| 



^a^(i+6),mLp — vnlnXX^ \^ ■ //* //H2 
^^_i>t mm |A„ - An , /i„ - /inj 



C(c/)(l + r)^ 



I T ° "t,r |UA,^(l+6),M?^;p• 



"T+5 



6f 
■ 1+6 



10.4.6. Step 10: Chain-rule and refined estimates on derivatives of /i""*"^. From 
Step 3 we have 



(10.82) 



bt 
1+6 



bt 
' 1+6 



and {via Proposition 14.271) 



(10.83) 



C(d)(l + r) 
^^"'r min {A* -Xh^^l- ^A} 



< 



" 1+6 



bt 
' 1+6 



< 



C{d) (1 



min {A* - Aj^ ; /i* - /i;^} 
Combining these bounds with Step 9, Proposition 14.241 and the identities 
f (V/i) ofi = (Vfi)"^ V(/i o fi) 

(vv/i) o = (vi])'2 vv(/i oQ)- (yn)-^ V^n {vn)~\vh o n), 



(10.84) <^ 
U 
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we get 
(10.85) 



C(c/)(l + r) 



< 



min{A^ - ; /x* - /xl} 



^A* (l + b),M^;l 



jA?j(l+6),M,*i;l 



and 
(10.86) 



;^A^{l + 6),nt,;l 



v\K^' o f^ig + ||v^i^i:.|L ||(v/.r^) o ^vIL A*,,,,,t, 



6t 



< 



C(d) (1 



in {A^ - A^ ; /i* - 



mm 



?-*'n(l + '>).Mn;P ■ 
1 + 6 



This gives US the bounds ||(V/i"+^)ofi" || = 0(l+r), || (V2/i"+i)ofi"|| = 0((l+r)2), 
which are optimal if one does not distinguish between the x and v variables. We 
shall now refine these estimates. First we write 

v{K^' o ni;) - o = v{ni, - id ) ■ [{vk^') o , 

and we deduce {via Propositions 14.24 1 and 14.271) 
(10.87) 



r-^n(H-'').Mn;P 



1 + 6 



< ||v(i]" -id)|L t \\iVK^')on- 



< C{d) 



(1+6), MS 

1 + r 



7-^n(l + 6),Mn;P 
1 + 6 



< 



min {A^ - AI ; /i* - /x^} 



t-, / -^"-^ Il2^n(l+'')'''n 11""^ " "t,r||2A^j{l+6),M?i;p 



6t 

^^t+t; 



T+E 



mm 



E 



4 



(i+r)-2 l|/^r^°^^i:.ll^A^(i+6),.,.. 
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(Note: r2" — Id brings the time-decay, while /i""*"^ brings the smallness.) 

This shows that (V/i"+^) o fi" ~ V(/i'"+^ o fi") as r oo. In view of Step 9, this 
also implies the refined gradient estimates 



(lOi 



T+E 



((V. + rV.)/i;'+^) o^i^ 



<(7||/l^ ° ^t,r |LA^(l+6),M?j;p; 



with 

C = C{d) 



6k 



min {A^ - Aii ; /i* - /il}^ (27r(Afc - A*))6 

10.4.7. Conclusion. Given A„+i < A*, fin+i < l^n^ we define 

An+l = -^n) Ain+l 



min{A^ - A^^ ; /i^ - /in} 



and we impose 



\* _ \t _ \t _ \t' =),''_ \t 

' '1-1 ' ' '1-1 ' *n ' ' * T) 



,1 



A* - A 



71+1 



Then from ffTirfKD . ffTTTTfL ffTirfSl) . ffTim . ffUM . ffTMTl) and ffHrssD we see that 
(E;J+^), (E;J+^), (Eh+^), (Eh+^) have all been established in the present subsection, 
with 
(10.89) 



5n 



+1 



min {A;; - A„+i ; - Hn+iY 



max 




> (27r(A, - A*))^ 



10.5. Convergence of the scheme. For any n > 1, we set 



(10.90) 



A,, — A! 



A! - A, 



+1 — /^n — /i„ — /^n " /^n+1 



A 



for some A > 0. By choosing A small enough, we can make sure that the conditions 
27r(Afe — A^) < 1, 27r(/ifc — /i^) < 1 are satisfied for all fc, as well as the other smallness 
assumptions made throughout this section. Moreover, we have A^ — A^ > A/Zc^, so 
conditions (Ci) to (Cg) will be satisfied if 

x:*"4<A«-. i:/*,<A^-(^-;^). 

fc=l 7=fc+l ^ ^ 
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for some small explicit constant lu > 0, depending on the other constants appearing 
in the problem. Both conditions are satisfied if 

CXj 

(10.91) ^ A;i2 < cu. 



k=i 



7+7 



Then from (110.761) we have T„ < C-y (n^/A)T-i, go the induction relation on 6„ 
allows 

2 \ 9 14+27 



2 



(10.92) 6i<C6, 5n+i = cij^\ e^("/^)^5: 

To establish this relation we also assumed that 6n is bounded below hj CpCn, the 
error coming from the short-time iteration; but this follows easily by construction, 
since the constraints imposed on 5„ are much worse than those on (n- 

Having fixed A, we will check that for 6 small enough, (110.921) implies both the 
fast convergence of {6k)kef>s, and the condition (110. 9ip . which will justify a posteriori 
the derivation of (110.921) . (An easy induction is enough to turn this into a rigorous 
reasoning.) 

For this we fix a G (1, aj), < z < zj < 1, and we check by induction 

(10.93) Vn>l, 5„<A^"". 

If A is given, (110.931) holds for n = 1 as soon as 5 < (A/C) z"". Then, to go from 
stage n to stage n + 1, we should check that 



/^„18 28+47 14+27 

"' ^Cn'^F^ /A~7^ ^2 ^2 a" ^ ^ ^a"+l 



this is true if 



\ (J / 28+47 14+27 



19 ^Cn^FT~/A^^ (2-a)a" 



Since a < 2, the supremum on the right-hand side is finite, and we just have to choose 
A small enough. Then, reducing A further if necessary, we can ensure (110.911) . This 
concludes the proof. 

Remark 10.1. This argument almost fully exploits the bi-exponential convergence 
of the Netwon scheme: a convergence like, say, 0(e~" ), would not be enough to 
treat values of 7 which are close to 1. In Subsection 111.21 we shall present a more 
cumbersome approach which is less greedy in the convergence rate, but still needs 
convergence like 0(e~"") for a large enough. 
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11. Coulomb/Newton interaction 

In this section we modify the scheme of Section [10] to treat the case 7 = 1. We 
provide two different strategies. The first one is quite simple and will only come 
close to treat this case, since it will hold on (nearly) exponentially large times in the 
inverse of the perturbation size. The second one, somewhat more involved, will hold 
up to infinite times. 

11.1. Estimates on exponentially large times. In this subsection we adapt the 
estimates of Section [10] to the case 7 = 1, under the additional restriction that 
< t < for some constant A>1. 

In the iterative scheme, the only place where we used 7 > 1 (and not just 7 > 1) 
is in Step 6, when it comes to the echo response via Theorem 17.71 Now, in the case 
7 = 1, the formula for K]^ should be 

n 

K^{t,T) = Y,6,K[''-'^'\t,T), 

k=l 

with = 2n min{(/ifc — fin)/2; (A^ — A*)/2; 77}. By Theorem 17.71 (ii) this in- 
duces, in addition to other well-behaved factors, an uncontrolled exponential growth 
0(e'"*), with 

n ^ 

— r \^ ^fc 

k=l 

in particular e„ will remain bounded and 0{6) throughout the scheme. 
Let us replace fll0.90p by 

\ \ * \* \ * * 

n (log(e -|- njj^ 

where A > is very small. (This is allowed since the series ^l/(n(log(e -|- n))^) 
converges — the power 2 could of course be replaced by any r > 1.) Then during the 
first stages of the iteration we can absorb the 0{e'^^^) factor by the loss of regularity 
if, say, 

^ A 
" ~ 2n (log(e -I- n)y 
Recalling that e„, = 0(5), this is satisfied as soon as 

K 
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where > is a positive constant depending on the other parameters of the problem 
but of course not on 6. So during these first stages we get the same long-time 
estimates as in Section [TOl 

For n > N we cannot rely on the loss of regularity any longer; at this stage the 
error is about 

N 

Sn C 6"" , 

where 1 < a < 2. To get the bounds for larger values of n, we use impose a restriction 
on the time-interval, say < t < Tmax- Allowing a degradation of the rate 5"" into 
6-" with a < a, we see that the new factor e'^"-^™'''' can be eaten up by the scheme if 

This is satisfied if 

max "-^ 1 it ^ 

Recalling f lll.ip . we see that the latter condition holds true if 

f ' log 1 

T^ax = O 

for some well-chosen constant A> 1. Then we can complete the iteration, and end 
up with a bound like 

Wft-m^.',,' <C6 Vt G [0,T^ax], 

where C is another constant independent of b. The conclusion follows easily. 

11.2. Mode-by-mode estimates. Now we shall change the estimates of Section 
[TOl a bit more in depth to treat arbitrarily large times for 7 = 1. The main idea is to 
work mode by mode in the estimate of the spatial density, instead of looking directly 
for norm estimates. 

Steps 1 to 5 remain the same, and the changes mainly occur in Step 6. 

Substep 6(a) is unchanged, but we only retain from that substep 

(11.2) VfeZ^ e2-(^^*+^")l^l |(ar'")lf)| < 



(7r(A„-A;))2- 



Substep b is more deeply changed. Let /i„ < /i* . 
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• First, for each £ e Z'*, we have, by Propositions 14.341 14.351 and the last part of 
Proposition 16.21 



jA* {1+6),M„ 
1+6 



;i dr 



< 



f Yl ' /' o (Id + e{Qi^ - Id )) 



^A*(l+6),M„ 



■ 1+6 



+ h'^^ — m'-J'r.t -7^n{l + 6),Mn 

6t 6t 



1 rt 



< 




^0 



ill 2^»*(l+'')>f'n'l ll^"r ^ Id 



\\Gl 

,-27r(At* -p:„)|m'| 



^"T+T 



E 



27r(/^* -p;„)|£-m| 



P™_^, (vF[/i^i] o (Id + e{ni^ - Id )) 



^^T+^ 



dr dO 



^ 6t 



E 



g~27r«-/^n)K-m| g-27r(At* -/i„)|m'[ 



^ 1+6 m, jr 
g-27r(Mn-A'n)|m-m'-q| (y ^ 



c/r. 



where 



/2„ + A:(1 + b) 



ht 



-,A* (l + 6),;i* . 
6t 



For a < 1 we have 



E 



^—2TTa\l—m\ ^—2-na\m'\ ^—2-na\m—m'—q\ ^ 



C{d) 



-2nf\e-q\ 



a" 
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we can argue as in Substep 6(b) of Section [10] to get 

(11-3) e-(^"*+'^")i^i \8{tA)\ < j-Srrr^ f^o + E^) (± 



5 



* dr 



327r(A*T+/x„)|q| I 



(l + r) 

Next, we use again Proposition 14.341 and simple estimates to bound S: 



J2 e-2.{<-An)H \\p,_^{F[h-^+^])\\^^„dT, 

where (3n = A* (1 + 6)|r — bt/ {1 + b) \ + Jin. Reasoning as in Substep 6(b) of Section 
[TOl we arrive at 

^ g-27r(M* -Mn)|m| / g27r(A*r+Mn)K-m||^^^n.+lj'|^£ _ ^^|| 



(l + r)2- 



• Then we consider the "main contribution" which we decompose as in 

Section [TO} 

n 

— n.,n+l — n,n+l i \ ^ — n,n+l 
k=l 

and we write for k > 1: 
e'^('"*+^")|(a;;f+f(£)|< V f Klt\t,T)\\P,.^{F[h^,+'])\\^.,, .,dT 



fK'\t,T)\\PeiF[hr' 
Jo 
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where 



I 



ht 



I 



= sup sup 

0<r<i 0<T<t 



1 + b 

v4h'^onlrr')-{v,{h'^on^-')) 



^Afc(l+6),Mfc 

-T-bt/il + b) 



)\m\ 



1 + r 

-27r(Ai^-Mn)K-m| 

1 + li-m]"' 



n,k 



I 



\l{t-T)+rm\ 



and the formula for Kq' is unchanged with respect to Section [TOl 

Assuming /i^ = /i„+?7 (t—r)/{l+t) and reasoning as in Substep 6(b) of Section [TOl 
we end up with the following estimate on the "main term" : 



(11.5) e2-(^"*+^")l'l |(a;^;"+')(£)| 

mGZ'* ° \k=l J 



+Mn)|^-m| \p[h'^+^]{£-m)\dT 

]{i-m)\dT. 



Then Substep 6(c) becomes, with $(£,t) = p[/i;^+^f(£), 



^2lT{\*t+Jln)\e\ 



$(£,t)- /" K^{i,t-T)^{i,T)dT 

Jo 



E 



(l + r)^J 



with 



a„ = 27r min 



k=l 



n 



g27r(A- r+M„)K-m| _ t)\ dr 

+ [ K^{t,T)e^'''-^>+^"^\^^mi,T)\dT, 



k=l 
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Note that 

Then we can apply Theorem 17.121 and deduce (taking already into account, for 
the sake of hsibility of the formula, that '^dk and '^5k/{\k — A^) are uniformly 
bounded) 

where 



T^r, = C max 



1 j. 

2\ 3 



7-2 



If AJj < A* and < /!„ are chosen as before and e = 27r(A,* — A^), this implies a 
uniform bound on 



obtained with the formula above with 



Te,n = Tn = C max 



1 1 1 



A* — An A„ A* /i„ /i* /i* fin 



max 



r 5+7+2d d+7+1/2 4d+l\ 
1 7 ' 7-1/2 ' 3 J 



Then Steps 7 to 10 of the iteration can be repeated with the only modification 
that /i* is replaced by /ij^. 

The convergence (Subsection 110.51) works just the same, except that now we need 
more intermediate regularity indices 

fin+l = fii < < fin < fJ'n, 

the obvious choice being to let ^1}^ — fi^^ = fin — fin = ~ P'n- 
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Choosing A„ — A„+i and /i„ — /i„+i of the order of A/n^, we arrive in the end at 
the induction 

with 

■5 + 7 + 2rf d + 7 + 1/2 4d+l 



:il.6) eK7) :=2ci + 2max , , . /o ' q 

'7 7 — 1/2 3 

Then the convergence of the scheme (and a posteriori justification of all the as- 
sumptions) is done exactly as in Section [101 

12. Convergence in large time 

In this section we prove Theorem 12.61 as a simple consequence of the uniform 
bounds established in Sections [10] and [TTl 

So let f'^, L, W satisfy the assumptions of Theorem 12.61 To simplify notation we 
assume L = 1. 

The second part of Assumption (12.141) precisely means that G C^'^. We shall 
actually assume a slightly more precise condition, namely that for some p E [1, 00], 

(12.1) Yl 7^ W^vfllLn^') <Co< +00, Vp G [i,p]. 

(It is sufficient to take p = 1 to get Theorem 12. 6j but if this bound is available for 
some p > 1 then it will be propagated by the iteration scheme, and will result in 
more precise bounds.) Then we pick up A G (0, A), ;U G (0, fi), P > 0, j3' E (0, /5). By 
symmetry, we only consider nonnegative times. 

If fi is an initial datum satisfying the smallness condition (12.151) . then by Theorem 
14.201 we have a smallness estimate on \\fi — f^\\z>^',M.';p for allp G [l,p], A' < A, /i' < /i. 
Then, as in Subsection 14. 121 we can estimate the solution to the linearized equation 



:i2.2) 



'dth^ + V ■ V^h^ + F[h^] ■ V„/° = 



and we recover uniform bounds in Z^'™ spaces, for any AG (A, A), /i G (/i, /i). More 
precisely, 

(12.3) SUp||p[/l,l]||^3;,+, + SUp||/ll(t, OLa.m;^ <C'5, 

t>0 t>0 t 
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with C = C{d, A', A, /i', /i, W, /°) (this is of course assuming e in Theorem 12.61 to be 
small enough). 

We now set Ai = A', and we run the iterative scheme of Sections [U HTUHTT] for all 
n >2. If e is small enough, up to slightly lowering Ai, we may choose all parameters 
in such a way that 

■^fc; K ~ ^ -^oo > A, /ifc, /ifc ^ /ioo > /i; 

fe— ►oo fc^oo — 

then we pick up 5 > such that 

Aoo(l + B)B > fi'^ > fi, 

and we let b{t) = B/{1 + t). 

As a result of the scheme, we have, for all /c > 2, 

(12.4) sup ofiJ^;M| < 4, 

t>r>0 ' .-^ 

where Y1T=2 5 and VL^ is the scattering associated to the force field generated 

hy + . . . + h^. Choosing t = r in (112. 4p yields 

sup \\h^ II 

< 4. 

^-^ i+B+t 

By Proposition 14.171 this implies 

sup ll/l^ll „Aoo(l + S),Moo-Aoo(l+S)S;l < 

In particular, we have a uniform estimate on h'l in z^°°'^°°'^ . Summing up over k 
yields for / = /° + X]a:>i the estimate 

(12.5) sup||/(t, O-rL^oo,.'^;! <C(5. 

t>0 t 

Passing to the limit in the Newton scheme, one shows that / solves the nonlinear 
Vlasov equation with initial datum /j. (Once again we do not check details; to 
be rigorous one would need to establish moment estimates, locally in time, before 
passing to the limit.) This implies in particular that / stays nonnegative at all times. 

Applying Theorem 14.201 again, we deduce from (112.51) 

sup||/(t, ■)- f\lA...<C5- 



or equivalently, with the notation used in Theorem 12.61 



12.6) sup \\f{t,x-vt,v)-f{v)\\.<C5. 
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Moreover, p = j f dv satisfies similarly 

sup \\p{t, ■ )|| < C6. 

t>o 

It follows that \p{t,k)\ < (7(5e-27rAoo|fc|t g-27rMoo|fc|^ fQj, g^j^y k ^ 0. On the one hand, 
by Sobolev embedding, we deduce that for any r G N, 

on the other hand, multiplying p by the Fourier transform of VW, we see that the 
force F = F[f] satisfies 

(12.7) Vt>0, \/keZ\ < C(5e-2"^'l^l*e-2"^'l'=l, 

for some \' > X, p' > p. 

Now, from (112. 6p we have, for any {k, 77) G Z*^ x W^, and any t > 0, 



12.8) f{t,k,7] + kt)-f{r]) 



so 



(12.9) \f{t,k,r]) \ < |/°(r/ + H)| +C5e-2"'^'l'=le"2-^''^+'=*'. 
In particular, for any k ^ 0, and any 77 G W^, 

(12.10) 7(t,A;,r/) = 0(e-2-^'*). 

Thus / is asymptotically close (in the weak topology) to its spatial average g 
(/) = / / dx. Taking /c = in (112.81) shows that, for any 77 G M"^, 

(12.11) |5^(t,r/)-70(r/)|<C5e"2-^'l''l. 

Also, from the nonlinear Vlasov equation, for any 77 G M*^ we have 

g{t, r]) = 7(0, r/) - / V / F(r, x) ■ V,/(r, x, ^) e-^-"'^' ciw c^x dr 



^Tf^ 



/,(0,7/)-2z7r5^ /" F{r,i)-7^f{T,-i,7^)dT. 



Using the bounds (112.71) and (112. lOp . it is easily shown that the above time-integral 
converges exponentially fast as t ^ 00, with rate 0(e~^ *) for any A" < A', to its 
limit 



T,i)-vf{T, -i,7^)dT. 
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By passing to the limit in (112.111) we see that 

|fi^oo(r/)-7°(r/)|<C5e-2-^'l^l 
and this concludes the proof of Theorem 



13. NON-ANALYTIC PERTURBATIONS 

Although the vast majority of studies of Landau damping assume that the per- 
turbation is analytic, it is natural to ask whether this condition can be relaxed. As 
we noticed in Remark 13. 5[ this is the case for the linear problem. As for nonlinear 
Landau damping, once the analogy with KAM theory has been identified, it is any- 
body's guess that the answer might come from a Moser-type argument. However, 
this is not so simple, because the "loss of convergence" in our argument is much 
more severe than the "loss of regularity" which Moser's scheme allows to overcome. 

For instance, the second-order correction h"^ satisfies 

dth^ + V ■ V,h^ + F[f'] ■ + F[h^] ■ VJ' = -F[h'] ■ VX- 

The action of F[f^] is to curve trajectories, which does not help in our estimates. 
Discarding this effect and solving by Duhamel's formula and Fourier transform, we 
obtain, with S = -F[h^] ■ V^, = J h"^ dv, 

(13.1) f{t,k)c^ f K'^{t-T,k)f{T,k)dT 

Jo 

+2i7r [ Y^{k-i)W{k-i)p^{T,k-ti^'^[T4,k{t-T))dT+ f S{r,k,k{t-r)) dr. 

Jo Jo 



(The term with K'^ includes the contribution of Vt,/°.) 

Our regularity/decay estimates on will never be better than those on the so- 
lution of the free transport equation, i.e., hi{x — vt,v), where hi = fi — /°. Let us 
forget about the effect of in (113.11) . replace the contribution of S* by a decaying 
term A{kt). Let us choose d = 1 and assume hi{i, ■ ) = if £ 7^ ±1. For k > let 
us use the long-time approximation 

k{t -t) -t) l[o,t]{T) (ir ~ 6_]a_, c= hi{-l, s) ds = hi{-l,0); 



166 



C. MOUHOT AND C. VILLANI 



note that c 7^ in general. Plugging all these simplifications in (113.11) and choosing 
W{k) = l/\k\^^'^ suggests the a priori simpler equation 

(13.2) ^(t, k) = Aikt) + ^^^-^ ^ k + 1 

Replacing f(t,k) by ip(t,k)/A{kt) reduces to A = 1, and then we can solve this 
equation by power series as in Subsection 17.1. 3[ obtaining 

(13.3) (^(t,A;) ~ A(H)e(^'=*)^. 

With a polynomial deterioration of the rate we could use a regularization argument, 
but the fractional exponential is much worse. 

However, our bounds are still good enough to establish decay for Gevrey pertur- 
bations. Let us agree that a function / = f{x,v) lies in the Gevrey class Q'^, f >1, 
if \f{k^ri)\ = O (exp{—c\{k, 7])]^^'^)) for some c > 0; in particular means analytic. 
(An alternative convention would be to require the nth derivative to be 0{n\'^).) As 
we shall explain, we can still get nonlinear Landau damping if the initial datum fi 
lies in Q'^ for u close enough to 1. Although this is still quite demanding, it already 
shows that nonlinear Landau damping is not tied to analyticity or quasi-analyticity, 
and holds for a large class of compactly supported perturbations. 

Theorem 13.1. Let A > 0. Let = f^{v) > be an analytic homogeneous profile 
such that 

E ^ iiv:;riiLi(M^) < +00, 

and let W = W{x) satisfy \ W{k) \ = 0{l/\k\), such that Condition (L) from Suhsec- 
tion \2.2\ holds. Let v & (1,1 + 6') with 9 = l/^((i, 7), where ^ was defined in flll.6p . 
Let (3 > and let a < \jv. Then there is e> ^ such that if 

(^-r)(A;,r/)|e^l''l^'''e^l'=l''^)+ f f \{f, - f'){x,v)\ e^^^^ dv dx <e, 



:= sup 

k,r] 

then as t ^ +00 the solution f = f{t,x,v) of the nonlinear Vlasov equation on 
j^d ^ yjj^if^ interaction potential W and initial datum fi satisfies, for all r G N, 

V(fc,r/), \f{tXri) - foo{ri)\ = 0(56"^*"); 

||F(t, ■)||c^(T^) = 0(5e-^*") 
for some c > and some homogeneous Gevrey profile foo, where F stands for the 
self- consistent force. 
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Remark 13.2. In view of (113.31) . one may hope that the result remains true for 
6 = 2. Proving this would require much more precise estimates, including among 
other things a qualitative improvement of the constants in Theorem 14.201 (recall 
Remark 1123]). 

Remark 13.3. One could also relax the analyticity of f'^, but there is little incentive 
to do so. 

Sketch of proof of Theorem \13.1[ We first decompose K = fi — f^ , using truncation 
by a smooth partition of unity in Fourier space: 

n>0 n>0 

where is the Fourier transform. Each bump function Xn should be localized around 
the domain (in Fourier space) 

Dn = {n'' <\{k,ri)\<{n + lfY 

for some exponent K > 1; but at the same time J^~^{Xn) should be exponentially 
decreasing in v. To achieve this, we let 

Then J-'~^{Xn) = ^~^{^d„) 1 has Gaussian decay, independently of n; so there is a 
uniform bound on J J \ h^{x, v) \ e^'"' dv dx, for some /3 > 0. 

On the other hand, if {k,r]) G Dn and {k',r]') ^ (L)„_i UD„ UD„+i), then \k-k'\ + 
\v ~ v'\ ^ cn^~^ for some c > 0; from this one obtains, after simple computations. 

So (with constants C and c changing from line to line) 

\hnk,v)\ < Ce^'\'\^ e-^l"!* l(n-i)-<|(M)l<(n+2)- + C e-^"^'""-'' e-^d'^l+l^'D 
< C max{e-t("-i)^, e"^"^^""^^ } e-^-^^l+I^D, 

where 

A„~^(n + 2)-M)^. 

If K > 2 then 2{K - 1) > K/i^; so Wh'^Wyj,,.-^ < Ce-t"'^^\ Then we may apply 
Theorem 14.201 to get a bound on /i" in the space Z^"'^"'^ with A.„ = A„/2, at the 
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price of a constant exp(C {n + 2)*^^ ^/'^)^^_ Assuming Kv > [1 — 1/ u)K ^ i.e., u < 2, 
we end up with 

(13.4) l|/^ril^x„,.„,, = 0(e-"^^), A„ = ^. 

Then we run the iteration scheme of Section [H] with the following modifications: 
(1) instead of /;."(0, ■ ) = 0, choose h^{0, ■ ) = h^, and (2) choose regularity indices 
An ~ A„ which go to zero as n goes to infinity. This generates an additional error 

term of size 0(exp(— cn^)), and imposes that A„ — A„+i be of order n' 

When we apply the bilinear estimates from Section [6l we can take A — A to be of 

the same order; so a = a„ and be chosen proportional to n' 

Then the large constants coming from the time-response will be, as in Section [TTl 
of order n'^ e'^^'' , with g G N and r = [(1 — l/h')K + 1]^, and the scheme will still 
converge like 0{e~'^"'") for any s < K/u, provided that K/u > r, i.e., 

V 1 

The rest of the argument is similar to what we did in Sections [10] to [121 In the 
end the decay rate of any nonzero mode of the spatial density p is controlled by 

^g-cn=g-A„t< j^e-'="'j sup [exp(-cn") exp(-cn"(^"^) t) 

< C exp(-ct'/^), 

and the result follows since s/K is arbitrarily close to □ 

Remark 13.4. An alternative approach to Gevrey regularity consists in rewriting 
the whole proof with the help of Gevrey norms such as 

wfWct = E ^-^^^^ wfWr^ = Y^e^^m^'" urn. 

which satisfy the algebra property for any u > 1. Then one can hybridize these 
norms, rewrite the time-response in this setting, estimate fractional exponential 
moments of the kernel, etc. 

Remark 13.5. In a more general C^' context, we do not know whether decay holds 
for the nonlinear Vlasov-Poisson equation. Speculations about this issue can be 
found in [50] where it is shown that (unlike in the linearized case) one needs more 
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than one derivative on the perturbation. As a first step in this direction, we mention 
that our methods imply a bound hke 0{6/{l + tY~^) for times t = 0(1/6), where f is 
a constant and r > r, as soon as the initial perturbation has norm 5 in a functional 
space W" involving r derivatives in a certain sense. The reason why this is nontrivial 
is that the natural time scale for nonlinear effects in the Vlasov-Poisson equation 
is not 0(1/6), but 0(1/ V6), as predicted by O'Neil [70] and very well checked in 
numerical simulations [5B]o Let us sketch the argument in a few lines. Assume 
that (for some positive constants c, C) 

(13.5) hi = '^h'l, ||/i"||2A„,A„;i < — , A„ = — . 

n 

Then we may run the Newton scheme again choosing a„ ~ cn/2", c„ = 0{^62~^^~'^^^^) 
and En = c' 6. Over a time-interval of length 0(1/5), Theorem I7.7( ii) only yields a 
multiplicative constant 0(e'^''*/^n) — 0(2^°'"). In the end, after Sobolev injection 
again, we recover a time-decay on the force F like 

5 ^ 2™'2 2""" e-^"* <C6 sup(2-"('^-"^) e"^"*) < - 



t)- 



as desired. Equation (113. 5p means that hi is of size 0(6) in a functional space 
whose definition is close to the Littlewood-Paley characterization of a Sobolev space 
with r derivatives. In fact, if the conjecture formulated in Remark 14.231 holds true, 
then it can be shown that W contains all functions in the Sobolev space ly+^O'^ 
satisfying an adequate moment condition, for some constant tq. 

14. Expansions and counterexamples 

A most important consequence of the proof of Theorem 12.61 is that the asymptotic 
behavior of the solution of the nonlinear Vlasov equation can in principle be deter- 
mined at arbitrary precision as the size of the perturbation goes to 0. Indeed, if we 
define g^{v) as the large-time limit of h'' (say in positive time), then \\g''\\ = 0{6k), 
so P + g\o + ■ ■ ■ + 9^ converges very fast to f^o- In other words, to investigate 
the properties of the time-asymptotics of the system, we may freely exchange the 
limits t — >■ oo and 5 — > 0, perform expansions, etc. This at once puts on rigor- 
ous grounds many asymptotic expansions used by various authors — who so far 
implicitly postulated the possibility of this exchange. 



1^ Passing from 0(l/^/^) to 0{l/5) is arguably an infinite-dimensional counterpart of Laplace's 
averaging principle, which yields stability for certain Hamiltonian systems over time intervals 
0( 1/(52) rather than 0(1/(5). 
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With this in mind, let us estimate the first corrections to the hnearized theory, in 
the regime of a very small perturbation and small interaction strength (which can 
be achieved by a proper scaling of physical quantities). We shall work in dimension 
d = 1 and in a periodic box of length L = 1. 

14.1. Simple excitation. For a start, let us consider the case where the perturba- 
tion affects only the first spatial frequency. We let 

• /°(^0 = e"'^"^: the homogeneous (Maxwellian) distribution; 

• e pi{x) = e cos(27rx): the initial space density perturbation; 

• e pi{x) 6{v): the initial perturbation of the distribution function; we denote by 
if the Fourier transform of 6; 

• aW: the interaction potential, with W{—x) = W{x). We do not specify its 
form, but it should satisfy the assumptions in Theorem 12.61 

We work in the asymptotic regime e — 0, a ^ 0. The parameter e measures 
the size of the perturbation, while a measures the strength of the interaction; after 
dimensional change, if W is an inverse power, a can be thought of as an inverse 
power of the ratio (Debye length) /(perturbation wavelength). We will not write 
norms explicitly, but all our computations can be made in the norms introduced in 
Section HI with small losses in the regularity indices — as we have done in all this 
paper. 

The first-order correction = 0{e) to /° is provided by the solution of the 
linearized equation (13. 3p . here taking the form 

dth^ + V + F[h^] ■ = 0, 

with initial datum h^{0, ■ ) = hi := fi — /°. As in SectionOwe get a closed equation 
for the associated density p[h^]: 

p[h^]{t,k) = hi{k,kt) -Aw^aWik) [ p{h\T,k)e-''^^^'-^^^\t-T)edT. 

Jo 

It follows that p|/i^](t, /c) = for A; 7^ ±1, so the behavior of pj/?-^] is entirely deter- 
mined by uiit) = p[h^](t, 1) and u_i(t) = f){h^](t, —1), which satisfy 

(14.1) ui{t) = -ip{t)-47r^aW{l) f ui{t) e-^^'-^^\t - t) dr 

2 Jo 

'^{t) + 0{a) . 
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(This equation can be solved explicitly [HI eq. 6], but we only need the expansion. 
Similarly, 



:i4.2) 



-7r(t-r)2 



(t - t) dr 



-2i7iaW{l)ui{t) 



The corresponding force, in Fourier transform, is given by F^{t, 1) 
and F\t,-1) = 2i7raW{l) U-i{t). 

From this we also deduce the Fourier transform of itself: 
(14.3) 

h\t,k,ri) = hiik,ri + kt)-A7r^aWik) I p[h^]{T,k) e-^^'^^^^'-^'^^^ {ri + k{t~T)) -kdr; 

Jo 

this is if /c 7^ ±1, while 

(14.4) h\t,l,r]) = -ip{7] + t) -Air^aWil) [ Ui(r) e-"(''+(*-"»' (r/ + (t - r)) dr 

2 In 



ip{ri + t) + 0{a) 



(14.5) 



h\t,-l,r]) = -^{r] -t) + 4n'aW{l) 



U- 



_^(^)e-(ry-(t-r))- (77 - (t - r)) c/r 



To get the next order correction, we solve, as in Section [TOl 

dth'^ + V ■ V^h^ + F[h^] ■ + F[h^] ■ (V^/° + VX) = -F[h^] ■ VX, 

with zero initial datum. Since h"^ = 0(5^), we may neglect the terms F[h^ ■ V„/i^ 
and ■ V^h^ which are both 0{ae^). So it is sufficient to solve 



:i4.6) 



dth'^ + V ■ + Flh'^] ■ V„f = -F[h^] ■ VX 



with vanishing initial datum. As t ^ 00, we know that the solution h'2{t,x,v) is 
asymptotically close to its spatial average (/jg) = j h'2 dx. Taking the integral over 
T'^ in ffT^ yields 

dt{K) = -{F[h^]-Vy). 
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Since converges to (hi), the deviation of / to (/,,) is given, at order e^, by 



+00 



{F[hi]-V^h^){t,v)dt 
p+00 

-^0 fcez 



Applying the Fourier transform and using (ll4.ip - (nA2l) - (ll4.41) - (ll4.5p . we deduce 



F[h^]{t,-1) {2iTTri) h\t, 1,7]) dt 

X) 

F[h^]{t,l){2iTT7])h\t,-l,r]) dt 





J ip{-t)(p{r] + t)dt- j (p{t)(p{r]-t)dt + 0{a) 

/ r+oo 

-7T^e^aW{l)r] i / <^{t) <^{r] - t) sign (t) dt + O {a) 



Summarizing: 

:i4.7) ^ 

\im f{t,k,r]) =0 if k 



t— ►oo 



hm fit, 0, v) = hit, 0, r/) - ^2 « (vr^ 1^(1)) r/ ( / (/.(t) (/.(r/ - t) sign (t) c/t + 0(« 



Since is an arbitrary analytic profile, this simple calculation already shows that 
the asymptotic profile is not necessarily the spatial mean of the initial datum. 
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Assuming e -C a, higher order expansions in a can be obtained by bootstrap on 
the equations (fm]) - (fm]) - (fmD - (11131) : for instance, 

/ + CO 
^{t)^{j]-t) sign(t) dt 
-CO 

- a' {2n'Wil)Yv ^ {^iv + t) ^(-r) - ^(r/ - t) ^(r)) e"-(*-^)^ (t - r) 

+ ifi-r) ipit) e-'^(''-(*-^))' {r,-{t- r)) | dr + 0{e'a''). 

What about the hmit in negative time? Reversing time is equivalent to changing 
f(t,x,v) into f(t,x,—v) and letting time go forward. So we define S{v) := —v, 
T{^){7]) := e^a7r^Wil)r]J^^^{tMv - t) sign (t) dt; then T(<^ o S) = T{ip) o 5, 
which means that the solutions constructed above are always homoclinic at order 
O(e^a). The same is true for the more precise expansions at order 0{e'^a'^), and in 
fact it can be checked that the whole distribution is homoclinic; in other words, 
/ is homoclinic up to possible corrections of order 0{e^). To exhibit heteroclinic 
deviations, we shall consider more general perturbations. 



14.2. General perturbation. Let us now consider a "general" initial datum fi{x, v) 
close to and expand the solution /. We write Sipkiji) = (/« ~ /") {k^v) ^'^^ 

= p[h^]. The interaction potential is assumed to be of the form a W with a <^ 1 
and W{x) = W{—x). The first equations of the Newton scheme are 

(14.8) p\t,k)=eipkikt)-A7r^aW{k) [ p\t, k) f {k{t - t)) {k^t - t) dr, 

Jo 

(14.9) 

h\t,k,ri) = eipk{v + kt)-^Ti^aW{k) I p^r, k) f {rj + k{t - r)) k- {r] + k{t-T)) dr, 

Jo 
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(14.10) 

h\t,k,7]) = -47i^aW{k) I p\r,k)p{r] + k{t-r))k-{r] + k{t-T))dT 



-An^a [ ^iy(£)pi(r,^)/ii(r,A;-£,r/ + A;(t-r))£- (r/ + A;(t-r))rfr 
Jo ^ 

- in^a [ J2 ^) {r,k-i,ri + k{t - t)) i ■ [r] + k{t - r)) dr 

Jo I 

- ATi'^a I V t?(£) f (r, l)h^{T,k-i,r] + k{t-T))£- {ri + k{t-T)) dr, 

Jo n 



(14.11) p\t,k) = -An^aWik) p\t, k) f {k{t - r)) {k^t - t)) dr 

Jo 

-An^a I Y^W{t)p\Tj)h^{T,k - l,k{t - t)) I ■ k{t - t) dr 
Jo ^ 

-Aii'^a [ Vl?(£)p2(r,£)/ii(r,A:-£,A:(t-r))£- A;(t-r)dr 
Jo ^ 

- Air^a [ J2 ^(^) ^) k-i,kit~T))£-k{t-T) dr. 
Jo ^ 

Here k and £ run over Z'^. 

From fll4.8l) - fll4.9l) we see that p^ and depend linearly on e, and 

(14.12) p\t,k) =e[Mkt) + 0{a)], h\t,k,r]) = e [ipk{r] + kt) + 0{a 



Then from fll4.10l) - fll4.1ip . and /i^ are 0(e^ a); so by plugging (114. 12p in these 
equations we obtain 
(14.13) 

f(t, k) = -47r V a /" V W{i) ipeiir) ipk-iikt-ir) i-k (t-r) dr +0{e^ a^)+0{e^ a), 
Jo ^ 

(14.14) 

h'^it, k, 7]) = -Air'^e^ " / MM ^k-iiv+kt-ir) £■ {ri+kit-r)) dr +0(6^ a^)+Oie^ a). 

Jo 
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We plug these bounds again in the right-hand side of (114.101) to find 
(14.15) h\t, 0, r/) = (E),(t, r/) + (n),(t, r/) + 0{e^ a^), 

where 

(E),(t,r/) = -47r2a / ■ ^)W{()p\T,l)h\T, -l,,^) dr 

Jo ^ 

is quadratic in e, and (JK)i,(t,T]) is a third-order correction: 
(14.16) 

(M), {t, 7]) =1671^6^ J2 W{i)W{m) 



t 

^m{ms) \ (pi-rni^T - uis) (p-e{r] - It) {i ■ m) {r - s) 



JO 



+ (feiir) if-e-miv — — ms) m ■ (rj — £{t — s)) > {£ ■ rj) ds dr. 



If is even, changing ifk into (pk{— ■ ) and t] into —r] amounts to change k into 
—k at the level of (I14.8p - (I14.9I) : but then (I)^ is invariant under this operation. 
We conclude that / is always homoclinic at second order in e, and we consider the 
infiuence of the third-order term (114.160 . Let 



C[vM := lim(n),(t,r/) 



t^oo 



After some relabelling, we find 



(14.17) C[^]{7]) = 1671^ ^ W{k)W{i) 

r>00 /-t 

M^r) lifk-iikt - ir) ip.kiv - kt) [k ■ £) {t - r) 



JO 



+ ipk{kt) ip^k-iiv - kt - £t)1 ■ {rj - k{t - r)) \ {k ■ r]) dr dt. 



Now assume that ip-k = a ipk with a = ±1. (a = 1 means that the perturbation is 
even in x; a = —1 that it is odd.) Using the symmetry {k,i) {—k, —i) one can 
check that 

C[ip o S] o S = a C[ipl 

where S{z) = —z. In particular, if the perturbation is odd in x, then the third-order 
correction imposes a heteroclinic behavior for the solution, as soon as C[!f] ^ 0. 
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To construct an example where C[ip\ 7^ 0, we set c? = 1, /° = Gaussian, fi — P = 
sin(27rx) 6'i(w) + sin(47rx) 62{v), ipi = —ip^i = 0i/2, = —^^2 = ^'2/2. The six pairs 



contributing to ([mT]) are (-1,1), (1,-1), (1,2), (2,1), (-1,-2), (-2,-1), 



By playing on the respective sizes of W{\) and 14^(2) (which amounts in fact to 
changing the size of the box), it is sufficient to consider the terms with coefficient 



If we let (yji and 1^2 vary in such a way that they become positive and almost con- 
centrated on M_|_, the only remaining term is the one in v'i(t) </22('7 + ^ ~t) </2i(^), and 
its contribution is negative for 77 > 0. So, at least for certain values of WiV) and 
VF(2) there is a choice of analytic functions (/Ji and 792, such that C\^\ 7^ 0. This 
demonstrates the existence of heteroclinic trajectories. 

To summarize: At first order in e, the convergence is to the spatial average; at 
second order there is a homoclinic correction; at third order, if at least three modes 
with zero sum are excited, there is possibility of heteroclinic behavior. 

Remark 14.1. As pointed out to us by Bouchet, the existence of heteroclinic trajec- 
tories implies that the asymptotic behavior cannot be predicted on the basis of the 
invariants of the equation and the interaction; indeed, the latter do not distinguish 
between the forward and backward solutions. 



We conclude this paper with some general comments about the physical implica- 
tions of Landau damping. 

Remark 114.11 show in particular that there is no "universal" large-time behavior 
of the solution of the nonlinear Vlasov equation in terms of just, say, conservation 
laws and the initial datum; the dynamics also have to enter explicitly. One can 
also interpret this as a lack of ergodicity: the nonlinearity is not sufficient to make 
the system explore the space of all "possible" distributions and to choose the most 
favorable one, whatever this means. Failure of ergodicity for a system of finitely 
many particles was already known to occur, in relation to the KAM theorem; this 



W{1)'^, i.e., the pairs (—1, 1) and (1, —1). Then the corresponding bit of C[ip]{ri) is 




+ <-Pi{t) (flit) (p2{r] + t - t) {r] + t - r) 
+ v?i(-r) V5i(?7 - t) Lp2{t + T){t- t) 

+ v?i(-'r) V9i(t) '^2{jl -t + T){t-T-ri) drdt. 



15. Beyond Landau damping 
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is mentioned e.g. in [571 p. 257] for the vortex system. There it is hoped that such 
behavior disappears as the dimension goes to infinity; but now we see that it also 
exists even in the infinite-dimensional setting of the Vlasov equation. 

At first, this seems to be bad news for the statistical theory of the Vlasov equation, 
pioneered by Lynden-Bell |53] and explored by various authors [16], [601 ES IMl [88] , 
since even the sophisticated variants of this theory try to predict the likely final 
states in terms of just the characteristics of the initial data. In this sense, our 
results provide support for an objection raised by Isichenko [321 P- 2372] against the 
statistical theory. 

However, looking more closely at our proofs and results, proponents of the statis- 
tical theory will have a lot to rejoice about. 

To start with, our results are the first to rigorously establish that the nonlinear 
Vlasov equation does enjoy some asymptotic "stabilization" property in large time, 
without the help of any extra diffusion or ensemble averaging. 

Next, the whole analysis is perturbative: each stable spatially homogeneous distri- 
bution will have its small "basin of damping" , and it may be that some distributions 
are "much more stable" than others, say in the sense of having a larger basin. 

Even more importantly, in Section [7] we have crucially used the smoothness to 
overcome the potentially destabilizing nonlinear effects. So any theory based on non- 
smooth functions might not be constrained by Landau damping. This certainly ap- 
plies to a statistical theory, for which smooth functions should be a zero-probability 
set. 

Finally, to overcome the nonlinearity, we had to cope with huge constants (even 
qualitatively larger than those appearing in classical KAM theory). If one believes 
in the explanatory virtues of proofs, these large constants might be the indication 
that Landau damping is a thin effect, which might be neglected when it comes to 
predict the "final" state in a "turbulent" situation. 

Further work needs to be done to understand whether these considerations apply 
equally to the electrostatic and gravitational cases, or whether the electrostatic case 
is favored in these respects; and what happens in "low" regularity. 

Although the underlying mathematical and physical mechanisms differ, nonlinear 
Landau damping (as defined by Theorem 12. 6p may arguably be to the theory of 
Vlasov equation what the KAM theorem is to the theory of Hamiltonian systems. 
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Like the KAM theorem, it might be conceptually important in theory and practice, 
and still be severely limitedo 

Beyond the range of application of KAM theory lies the softer, more robust weak 
KAM theory developed by Fathi [25] in relation to Aubry-Mather theory. By a nice 
coincidence, a Vlasov version of the weak KAM theory has just been developed by 
Gangbo and Tudorascu [28], although with no relation to Landau damping. Making 
the connection is just one of the many developments which may be explored in the 
future. 

Appendix 

In this appendix we gather some elementary tools, our conventions, and some 
reminders about calculus. We write Nq = {0, 1,2,.. .}. 

A.l. Calculus in dimension d. If n G Nq we define 

n\ = rill . . . nj 

and 

( ] = ( "^^ ) ( "''^ ] 
^ m ' ^ mi ' ' ' ' \ rrid ' ' 

If 2 G and n G Z^, we let 

\\z\\ = \zi\ + ... + \zd\] z"" = z'^'. . . z^" E C; kr=k"|. 

In particular, if 2; G we have 

II ^\\n 

^ n! 

We may write e'^' instead of e"^". 

A. 2. Multi-dimensional differential calculus. The Leibniz formula for func- 
tions /, ^ : M ^ M is 



m<n 



"'^^It is a well-known scientific paradox that the KAM theorem was at the same time tremendously 
influential in the science of the twentieth century, and so restrictive that its assumptions are 
essentially never satisfied in practice. 
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where of course /W = d'^f/dx''. The expression of derivatives of composed functions 
is given by the Faa di Bruno formula: 

J2jmj=n j=l 

These formulas remain valid in several dimensions, provided that one defines, for 
a multi-index n = (ni, . . . , Ud), 

Qni Qna 
fin) ^ f 

They also remain true ii [di, . . . , da) is replaced by a d-tuple of commuting deriva- 
tion operators. 

As a consequence, we shall establish the following Leibniz-type formula for oper- 
ators that are combinations of gradients and multiplications. 

Lemma A.l. Let f and g be functions of v E M.'^, and a,b E C^. Then for any 
n e N'^, 

(V. + (a + b)nfg) = E ( m ) + °)"^ ^"^^ + ^^^'"^3- 
Proof. The right-hand side is equal to 

E ( ) ( 7 ) ( "r ) 

m,q,r 

After changing indices p = q + r, s = m — q, this becomes 



s,p,r p 

= (V, + {a + b)nfg). 



□ 



A. 3. Fourier transform. If / is a function M'^ — > M, we define 
(A.l) m= [ e-'^-^-^f{v)dv- 

then we have the usual formulas 

f{v) = [ fiv) e"^-^-^ dTi; Vfiv) = 2z7rr/ /(r/). 
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Let Ti = M"'/(LZ'^). If / is a function M, we define 

(A.2) F\k)= I e-^^^i-^f{x)dx- 



then we have 

If / is a function TfxM.'^ ^ M, we define 
(A.3) p)^k,r])= [ [ e-^'^'i-^'e-^''''' '' fix,v)dxdv; 

SO that the reconstruction formula reads 

When L = 1 we do not specify it: so we just write 

(There is no risk of confusion since in that case, flA.Sp and flA.ip coincide.) 

A. 4. Fixed point theorem. The following theorem is one of the many variants 
of the Picard fixed point theorem. We write -8(0, R) for the ball of center and 
radius R. 

Theorem A.2 (Fixed point theorem). Let E be a Banach space, F : E E, and 
R = 2||F(0)||. If F is {l/2)-Lipschitz B{0, R) —>■ E, then it has a unique fixed point 
m B{0,R). 

Proof. Uniqueness is obvious. To prove existence, run the classical Picard iterative 
scheme initialized at 0: xq = 0, xi = F{0), X2 = F(F(0)), etc. It is clear that (x.„) is 
a Cauchy sequence and ||x„|| < ||F(0)||(1 + . . . + 1/2") < 2||F(0)||, so x„ converges 
in -8(0, R) to a fixed point of F. □ 
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